Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,280,527 A→G 24.2% 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,5270AG24.2% 57.6 / 15.8 33S332G (AGC→GGC) narKnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base A (11/14);  new base G (5/3);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCT  >  W3110S.gb/1280461‑1280596
                                                                  |                                                                     
gCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCGGCg                                                                    >  1:1269011/1‑70 (MQ=255)
 cATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGc                                                                  >  1:2688909/1‑71 (MQ=255)
     aTCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCGGCGGCctgc                                                              >  1:1482278/1‑71 (MQ=255)
     aTCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGGctg                                                               >  1:271992/1‑70 (MQ=255)
        tGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCctgct                                                             >  1:553454/1‑69 (MQ=255)
        tGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCctgct                                                             >  1:1837076/1‑69 (MQ=255)
        tGACCGTCTGGGCGAAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCctgct                                                             >  1:2285455/1‑69 (MQ=255)
           ccGTCTGGGCGGAACTCGTGTCACGCTGGTGAATTTTATTCTGATGGCGATTTTCGGCGGCCTGCTGTTcc                                                        >  1:1591526/1‑71 (MQ=255)
                  ggCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCt                                                  >  1:1561320/1‑70 (MQ=255)
                          tCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTc                                                         >  1:1534920/1‑55 (MQ=255)
                          tCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTc                                                         >  1:547466/1‑55 (MQ=255)
                              gtCACGCTGGTGAACTTTATTCTGATGGCGATTTTAGGCGGCCTGCTGTTCCt                                                       >  1:335469/1‑53 (MQ=37)
                               tCACGCTGGTGAACTTTATTCTGATGGCGATTTTCGGCGGCCTGCTGTTCCTGACCTTACCGACTGACggg                                    <  1:1032719/71‑1 (MQ=255)
                                 aCGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACggg                                    >  1:3993880/1‑69 (MQ=255)
                                       gTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTcc                                                        <  1:1714911/43‑1 (MQ=255)
                                       gTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTcc                                                        <  1:2819079/43‑1 (MQ=255)
                                       gTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCt                                                  <  1:2078911/49‑1 (MQ=255)
                                       gTGAACTTTATTCTGATGGCGACTTTCAGCGGCCTGCTGTTcc                                                        <  1:3318687/43‑1 (MQ=255)
                                        tGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACggg                                    <  1:3683086/62‑1 (MQ=255)
                                            cTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACt                                          <  1:1128766/52‑1 (MQ=255)
                                               tATTCTGATGGCGATTTTCGGCGGCCTGCTGTTGCTGACCTTACCGACTGACGGGCAGGGCGGAAGCtt                      <  1:1627149/69‑1 (MQ=255)
                                               tATTCTGATGGCGATTTTCGGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCtt                      <  1:42540/69‑1 (MQ=255)
                                               tATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCtt                      <  1:2046886/69‑1 (MQ=255)
                                               tATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCa                    >  1:2878330/1‑71 (MQ=255)
                                                 ttCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAg                         <  1:1467680/64‑1 (MQ=255)
                                                 ttCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAg                         <  1:3394000/64‑1 (MQ=255)
                                                     gATGGCGATTTTCGGCGGCCTGCTGTTCCTGACCTTAc                                               >  1:3973757/1‑38 (MQ=255)
                                                            aTTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGc                                   <  1:892610/43‑1 (MQ=255)
                                                            aTTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGc                                   <  1:263317/43‑1 (MQ=255)
                                                            aTTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGc                                   <  1:2611507/43‑1 (MQ=255)
                                                            aTTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGc                                   <  1:1506428/43‑1 (MQ=255)
                                                            aTTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGc                                   <  1:1440277/43‑1 (MQ=255)
                                                              attCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCa                                  >  1:25227/2‑42 (MQ=255)
                                                                  aGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCt  <  1:3582709/70‑1 (MQ=255)
                                                                  aGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCt  <  1:1527928/70‑1 (MQ=255)
                                                                  |                                                                     
GCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCT  >  W3110S.gb/1280461‑1280596

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