Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,541,247 C→A V1020F (GTC→TTC)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,541,2470CA76.6% 72.1 / 30.5 47V1020F (GTC→TTC) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base A (0/36);  total (6/41)
Fisher's exact test for biased strand distribution p-value = 4.30e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACG  >  W3110S.gb/1541198‑1541295
                                                 |                                                
ctgcACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGCCGTGATCGCTTTc                                      <  1:1695081/59‑1 (MQ=255)
ctgcACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGACGTGATCGCTTTc                                      <  1:1577931/59‑1 (MQ=255)
ctgcACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGACGTGATCGCTTTc                                      <  1:1693266/59‑1 (MQ=255)
ctgcACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGACGTGATCGCTTTc                                      <  1:404916/59‑1 (MQ=255)
ctgcACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGACGTGATCGCTTTc                                      <  1:1971746/59‑1 (MQ=255)
         aTTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTccagaccaga                             <  1:1349932/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTccagaccagag                            <  1:1821497/62‑2 (MQ=38)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:548553/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2675323/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2710347/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2783488/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2916075/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:327949/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:37602/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:399793/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:428276/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:500354/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2393338/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:584430/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:617878/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:619837/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:642291/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:708030/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:743950/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:858451/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:941071/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:1971397/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:1140920/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:1206225/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:1419763/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:1679953/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:1810098/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:207337/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2278437/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2287334/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2358011/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:1005977/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2515135/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:2547508/62‑1 (MQ=255)
          ttAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCTAGACCAGAc                            <  1:1732754/62‑1 (MQ=38)
            aaCTCATGAACGTTGGTGTATCCCGCGTTATAGGAGAAGTGATCGCTTTCcagaccagac                            <  1:1443191/60‑1 (MQ=37)
                                    gcgTTATAGGAGCCGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTAcg  >  1:2013392/1‑62 (MQ=255)
                                    gcgTTATAGGAGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTAcg  >  1:354328/1‑62 (MQ=255)
                                    gcgTTATAGGAGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTAcg  >  1:2614576/1‑62 (MQ=255)
                                    gcgTTATAGGAGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTAcg  >  1:1706222/1‑62 (MQ=255)
                                    gcgTTATAGGAGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTAcg  >  1:159226/1‑62 (MQ=255)
                                    gcgTTATAGGAGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTAcg  >  1:1143842/1‑62 (MQ=255)
                                                 |                                                
CGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACG  >  W3110S.gb/1541198‑1541295

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: