Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,280,527 A→G 17.6% S332G (AGC→GGC)  narK → nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,280,5270AG17.6% 74.5 / 10.4 34S332G (AGC→GGC) narKnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base A (15/13);  new base G (2/4);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCG  >  W3110S.gb/1280460‑1280583
                                                                   |                                                        
tgCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCg                                                       <  1:1794123/71‑1 (MQ=255)
      aTCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCGGCGGCctgc                                                 >  1:4030064/1‑71 (MQ=255)
      aTCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCGGCGGCctgc                                                 >  1:2510106/1‑71 (MQ=255)
      aTCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCctgc                                                 >  1:2724990/1‑71 (MQ=255)
      aTCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCctgc                                                 >  1:4023347/1‑71 (MQ=255)
          gACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGtt                                             >  1:4066810/1‑71 (MQ=255)
           aCCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGt                                              <  1:1333968/69‑1 (MQ=255)
             cGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCGGCGGCCTGCTGTTCCt                                          <  1:2669483/71‑1 (MQ=255)
              gTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTg                                         <  1:3133584/71‑1 (MQ=255)
                 tGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTg                                         <  1:2391074/68‑1 (MQ=255)
                   ggCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCGGCGGCCTGCTGTTCCTGACCtt                                    <  1:2205856/71‑1 (MQ=255)
                   ggCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCg                                                       >  1:3208270/1‑52 (MQ=255)
                     cGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCctgctg                                               <  1:1882849/58‑1 (MQ=255)
                     cGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTAc                                  >  1:196762/1‑71 (MQ=255)
                       gAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGa                                        <  1:363641/63‑1 (MQ=255)
                           tCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGa                                        <  1:2396740/59‑1 (MQ=255)
                           tCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACt                             >  1:299056/1‑70 (MQ=255)
                           tCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACt                             >  1:1814012/1‑70 (MQ=255)
                              tgtCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACg                         <  1:1361466/71‑1 (MQ=255)
                               gtCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACgg                        >  1:3003538/1‑71 (MQ=255)
                                tCACGCTGGTGAACTTTATTCTGATGGCGATTTTCGGCGGCCTGCTGTTCCTg                                         <  1:553423/53‑1 (MQ=255)
                                tCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTa                                   >  1:1893778/1‑59 (MQ=255)
                                        gTGAACTTTATTCTGATGGCGATTTTCGGCGGCctgc                                                 <  1:2992899/37‑1 (MQ=255)
                                        gTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTAcc                                 <  1:595080/53‑1 (MQ=255)
                                        gTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACggg                       <  1:2531351/63‑1 (MQ=255)
                                          gAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCa                     >  1:3070920/1‑63 (MQ=255)
                                              tttATTCTGATGGCGATTTTCAGCGGCCTGCTGtt                                             <  1:945874/35‑1 (MQ=255)
                                              tttATTCTGATGGCGATTTTCAGCGGCCTGCTGtt                                             <  1:4490501/35‑1 (MQ=255)
                                               ttATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCa                     >  1:807178/1‑58 (MQ=255)
                                               ttATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCtt         >  1:1443538/1‑70 (MQ=255)
                                                    cTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTATCGACTGACGGGCAGGGCGGAAGCTTCATGGc   >  1:3129800/1‑71 (MQ=255)
                                                    cTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGc   >  1:2061618/1‑71 (MQ=255)
                                                     tGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTa                                   >  1:1207136/1‑38 (MQ=255)
                                                       aTGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCg  <  1:3943518/69‑1 (MQ=255)
                                                                   |                                                        
TGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCG  >  W3110S.gb/1280460‑1280583

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: