Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,337,480 A→G 76.0% S5S (AGT→AGC wcaF ← predicted acyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,337,4800AG76.0% 137.9 / 41.1 79S5S (AGT→AGCwcaFpredicted acyl transferase
Reads supporting (aligned to +/- strand):  ref base A (15/4);  new base G (30/30);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 3.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

TACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGT  >  Exported/1337414‑1337545
                                                                  |                                                                 
tACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCtt                                                                >  1:2957886/1‑70 (MQ=255)
  cttcccACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTaa                                                              <  1:3579935/67‑1 (MQ=255)
  cTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTaaa                                                             <  1:4114445/71‑1 (MQ=255)
  cTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTaa                                                              >  1:46917/1‑70 (MQ=255)
  cTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTaa                                                              <  1:4368718/70‑1 (MQ=255)
       caccaTAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTg                                                        >  1:156448/1‑71 (MQ=255)
       caccaTAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTg                                                        <  1:4302375/71‑1 (MQ=255)
           aTAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCtt                                                         >  1:426087/1‑66 (MQ=255)
           aTAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCAtt                                                    <  1:947865/71‑1 (MQ=255)
           aTAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCAtt                                                    >  1:2073851/1‑71 (MQ=255)
                 gCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTaaa                                                             <  1:5266441/56‑1 (MQ=255)
                   aCTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTcc                                            >  1:3385274/1‑71 (MQ=255)
                   aCTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTcc                                            >  1:3237177/1‑71 (MQ=255)
                   aCTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGttttac                                            >  1:3395679/1‑69 (MQ=255)
                   aCTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTcc                                            >  1:1739101/1‑71 (MQ=255)
                    cTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCt                                           <  1:4096020/71‑1 (MQ=255)
                    cTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCt                                           <  1:623693/71‑1 (MQ=255)
                      ttAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCtt                                          >  1:4441986/1‑70 (MQ=255)
                      ttAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCtt                                          >  1:4993869/1‑70 (MQ=255)
                      ttAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCtt                                          >  1:3273173/1‑70 (MQ=255)
                      ttAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCTta                                         >  1:2446581/1‑71 (MQ=255)
                       tAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTta                                         <  1:883486/70‑1 (MQ=255)
                       tAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTta                                         >  1:1927161/1‑70 (MQ=255)
                       tAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTta                                         >  1:1676737/1‑70 (MQ=255)
                       tAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCTta                                         >  1:1659535/1‑70 (MQ=255)
                          tAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCTTATAtt                                     >  1:60187/1‑71 (MQ=255)
                            gCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGtt                                                >  1:5695464/1‑58 (MQ=255)
                            gCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGtt                                                >  1:1621615/1‑58 (MQ=255)
                             cGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGttt                                               <  1:4005145/58‑1 (MQ=255)
                             cGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGttt                                               <  1:1717938/58‑1 (MQ=255)
                                  ccgccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACtt                              <  1:5000553/70‑1 (MQ=255)
                                   cgccCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTg                            <  1:1352867/71‑1 (MQ=255)
                                     ccccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATAtt                                     <  1:3722712/60‑1 (MQ=255)
                                     ccccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGt                           <  1:4801746/70‑1 (MQ=255)
                                     ccccGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGt                           >  1:1128681/1‑70 (MQ=255)
                                      cccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTta                         <  1:5070010/71‑1 (MQ=255)
                                       ccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCtt                                          <  1:2313482/53‑1 (MQ=255)
                                        cGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGtt                                                <  1:4883674/46‑1 (MQ=255)
                                          gAACCCTTTCGGCACCGAGAAACCGCTTAAATCtt                                                         >  1:306495/1‑35 (MQ=255)
                                          gAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGt                           >  1:2010454/1‑65 (MQ=255)
                                          gAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTtatat                      >  1:5034088/1‑70 (MQ=255)
                                               cTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTtat                        >  1:1604682/1‑63 (MQ=255)
                                               cTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTtat                        >  1:3186818/1‑63 (MQ=255)
                                               cTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTtat                        >  1:3392318/1‑63 (MQ=255)
                                               cTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCt                <  1:2199534/71‑1 (MQ=255)
                                                tttCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCtt               >  1:1147308/1‑71 (MQ=255)
                                                  tCGGCACCGAGAAACCGCTTAAATCTTGCATTGGtt                                                >  1:1327646/1‑36 (MQ=255)
                                                  tCGGCACCGAGAAACCGCTTAAATCTTGCATTGGtt                                                >  1:310229/1‑36 (MQ=255)
                                                  tCGGCACCGAGAAACCGCTTAAATCTTGCATTGGtt                                                >  1:2996442/1‑36 (MQ=255)
                                                  tCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCtt               <  1:59025/69‑1 (MQ=255)
                                                     gCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTc          <  1:865400/71‑1 (MQ=255)
                                                      cACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCt         <  1:199590/71‑1 (MQ=255)
                                                       aCCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCtt        >  1:606810/1‑71 (MQ=255)
                                                        ccGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATAtt                                     <  1:3264361/41‑1 (MQ=255)
                                                        ccGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACtt                              <  1:5527516/48‑1 (MQ=255)
                                                         cGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCttt       >  1:3898029/1‑70 (MQ=255)
                                                         cGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTg      <  1:1943437/71‑1 (MQ=255)
                                                         cGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTg      <  1:4554737/71‑1 (MQ=255)
                                                          gagaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGtt                          <  1:3256696/50‑1 (MQ=255)
                                                          gagaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGtt                          <  1:3249902/50‑1 (MQ=255)
                                                          gagaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGtt                          <  1:1602760/50‑1 (MQ=255)
                                                          gagaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTtatat                      <  1:4879406/54‑1 (MQ=255)
                                                          gagaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCttt       >  1:2971226/1‑69 (MQ=255)
                                                          gagaAACCACTTAAATCCTGCATTGGTTTTCCTTATATTCAGACttt                             <  1:3724934/47‑1 (MQ=255)
                                                           agaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGa                                 >  1:2743975/1‑42 (MQ=255)
                                                           agaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGa                                 >  1:2737287/1‑42 (MQ=255)
                                                           agaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTtat                        >  1:2873613/1‑51 (MQ=255)
                                                           agaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTtat                        >  1:587810/1‑51 (MQ=255)
                                                            gaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAgg   <  1:2334129/71‑1 (MQ=255)
                                                            gaAACCACTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGGTTTCGTCTTTGAgg   >  1:2792955/1‑71 (MQ=255)
                                                            gaAACCACTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAgg   >  1:5684829/1‑71 (MQ=255)
                                                            gaAACCACTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAgg   >  1:2511516/1‑71 (MQ=255)
                                                             aaaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACttt                             <  1:3478123/44‑1 (MQ=255)
                                                             aaaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCttt              >  1:4011231/1‑59 (MQ=255)
                                                             aaaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAgg   >  1:2647933/1‑70 (MQ=255)
                                                             aaaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGgt  <  1:5440015/71‑1 (MQ=255)
                                                              aaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAgg   >  1:447298/1‑69 (MQ=255)
                                                              aaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAgg   >  1:1878728/1‑69 (MQ=255)
                                                              aaCCGCTTAAATCTTGCATTGGTTATCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAg    >  1:24230/1‑68 (MQ=255)
                                                                  |                                                                 
TACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCACTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGT  >  Exported/1337414‑1337545

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: