Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,866,764 G→A 7.3% Q216* (CAG→TAG)  nlpD ← predicted outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,866,7640GA 7.3% 174.5 / 3.4 55Q216* (CAG→TAG) nlpDpredicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (22/29);  new base A (2/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

ACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATTTTGTGCA  >  W3110S.gb/2866697‑2866831
                                                                   |                                                                   
aCCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTg                                                                     >  1:1948883/1‑68 (MQ=255)
   gTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTc                                                                  <  1:1323385/68‑1 (MQ=255)
      aCAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCAcc                                                               >  1:448591/1‑68 (MQ=255)
          gcgcTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAg                                                          >  1:1013444/1‑69 (MQ=255)
          gcgcTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAg                                                          <  1:2764314/69‑1 (MQ=255)
          gcgcTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAg                                                          >  1:1784819/1‑69 (MQ=255)
               gtgACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTc                                                                  >  1:1103770/1‑56 (MQ=255)
               gtgACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTc                                                                  >  1:1128985/1‑56 (MQ=255)
               gtgACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTc                                                                  <  1:3405864/56‑1 (MQ=255)
               gtgACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTc                                                                  >  1:3276844/1‑56 (MQ=255)
               gtgACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTc                                                                  <  1:3227746/56‑1 (MQ=255)
                  aCCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAg                                                    <  1:89885/67‑1 (MQ=255)
                  aCCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAg                                                    <  1:2295499/67‑1 (MQ=255)
                  aCCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGa                                                   <  1:1657741/68‑1 (MQ=255)
                     gTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTATTCACCCGAAGACTCAGAATa                                                >  1:933373/1‑68 (MQ=255)
                       ggTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGt                                             <  1:1749308/69‑1 (MQ=255)
                       ggTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGt                                             <  1:1395470/69‑1 (MQ=255)
                       ggTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGt                                             <  1:645690/69‑1 (MQ=255)
                       ggTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGt                                             <  1:2258769/69‑1 (MQ=255)
                        gTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGt                                             <  1:147782/68‑1 (MQ=255)
                        gTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTa                                            >  1:2178826/1‑69 (MQ=255)
                        gTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTa                                            <  1:1840051/69‑1 (MQ=255)
                        gTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTa                                            >  1:462736/1‑69 (MQ=255)
                        gTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTa                                            <  1:2550914/69‑1 (MQ=255)
                        gTCGCAGTTGGCTTGTTGTTCGGCAACATTCTGTTAGCACTCTATTCACCCGAAGACTCAGAATACGTa                                            >  1:1344450/1‑69 (MQ=255)
                           gCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAg                                                    <  1:837680/58‑1 (MQ=255)
                           gCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAg                                                    <  1:1688916/58‑1 (MQ=255)
                           gCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAg                                                    <  1:3641410/58‑1 (MQ=255)
                           gCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAAtt                                         >  1:1717466/1‑69 (MQ=255)
                           gCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAAtt                                         >  1:3107611/1‑69 (MQ=255)
                                 ggCTTGTTGTTCGGCAACATTTTGTTAGCACTCTg                                                                     >  1:2196710/1‑35 (MQ=255)
                                 ggCTTGTTGTTCGGCAACATTTTGTTAGCACTCTg                                                                     >  1:1419550/1‑35 (MQ=255)
                                 ggCTTGTTGTTCGGCAACATTTTGTTAGCACTCTg                                                                     >  1:90853/1‑35 (MQ=255)
                                 ggCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGt                                   >  1:1291478/1‑69 (MQ=255)
                                    ttgttgTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGa                                                   <  1:1630042/50‑1 (MQ=255)
                                    ttgttgTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGc                                <  1:2655644/69‑1 (MQ=255)
                                    ttgttgTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGc                                <  1:249688/69‑1 (MQ=255)
                                    ttgttgTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGc                                <  1:2447684/69‑1 (MQ=255)
                                     tgttgtTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAg                                                    <  1:283647/48‑1 (MQ=255)
                                     tgttgtTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAg                                                    <  1:3089092/48‑1 (MQ=255)
                                         gttCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGc                           <  1:177053/69‑1 (MQ=255)
                                         gttCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGc                           <  1:322490/69‑1 (MQ=255)
                                                   aTTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGc                           <  1:2259803/59‑1 (MQ=255)
                                                   aTTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGc                           >  1:385934/1‑59 (MQ=255)
                                                   aTTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGc                           >  1:69782/1‑59 (MQ=255)
                                                     tttGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACgg                >  1:3111429/1‑68 (MQ=255)
                                                      ttGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGt               >  1:3556980/1‑68 (MQ=255)
                                                      ttGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGt               >  1:3004756/1‑68 (MQ=255)
                                                         ttAGCACTCTATTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGa            <  1:1105686/68‑1 (MQ=255)
                                                         ttAGCACTCTATTCACCCGAAGACTCAGAATACGTAATTGTCGGTGGCGACGCAACAGCAACGGTGGa            <  1:2396890/68‑1 (MQ=255)
                                                          tAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGaa           <  1:3436327/68‑1 (MQ=255)
                                                            gCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAAttt        >  1:789840/1‑69 (MQ=255)
                                                              aCTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAAttt        <  1:843238/67‑1 (MQ=255)
                                                                tctGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATTTtgtg    <  1:2264317/69‑1 (MQ=255)
                                                                 ctGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATTTTGTGc   >  1:2665694/1‑69 (MQ=255)
                                                                 ctGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATTTTGTGc   >  1:3684394/1‑69 (MQ=255)
                                                                 ctGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATTTTGTGc   >  1:1691274/1‑69 (MQ=255)
                                                                  tGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATTTTGTGCa  >  1:2483909/1‑69 (MQ=255)
                                                                  tGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATTTTGTGCa  >  1:734516/1‑69 (MQ=255)
                                                                   |                                                                   
ACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATTTTGTGCA  >  W3110S.gb/2866697‑2866831

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: