Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,538,270 A→G 69.2% F21S (TTT→TCT)  narV ← nitrate reductase 2 (NRZ), gamma subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,538,2700AG69.2% 33.4 / 18.1 26F21S (TTT→TCT) narVnitrate reductase 2 (NRZ), gamma subunit
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base G (10/8);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCC  >  W3110S.gb/1538204‑1538336
                                                                  |                                                                  
aTTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAacac                                                                <  1:2575226/71‑1 (MQ=255)
         gaggagGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGcac                                                        >  1:1471176/1‑70 (MQ=255)
             aggCGCGCCAGGTGTACTGCCCTTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATg                                                   <  1:1782119/71‑1 (MQ=255)
             aggCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATg                                                   <  1:356417/71‑1 (MQ=255)
                    ccAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGac                                                                  <  1:1583198/49‑1 (MQ=255)
                       ggTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACggg                                             >  1:2013732/1‑67 (MQ=255)
                          gtACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGata                                       >  1:1543584/1‑70 (MQ=255)
                             cTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCg                                    >  1:1520383/1‑70 (MQ=255)
                             cTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCg                                    >  1:1108229/1‑70 (MQ=255)
                                cccGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaa                                 >  1:930342/1‑70 (MQ=255)
                                cccGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaa                                 >  1:1036686/1‑70 (MQ=255)
                                    tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGacac                                                                <  1:1725103/35‑1 (MQ=255)
                                    tAGTCGTAGCGCAGCCAGCTGACGAGGAAAAACACCGTCGGACaa                                                      <  1:1403901/45‑1 (MQ=38)
                                         gTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGt                       >  1:1158327/1‑71 (MQ=255)
                                         gTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGt                       >  1:1882470/1‑71 (MQ=255)
                                            gcgcAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCa                    <  1:3911/71‑1 (MQ=255)
                                              gcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGtt                      >  1:1309193/1‑67 (MQ=255)
                                                  ccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGa                          >  1:2472246/1‑59 (MQ=255)
                                                      cTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGa                          <  1:242209/55‑1 (MQ=255)
                                                      cTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGa                          <  1:2030361/55‑1 (MQ=255)
                                                      cTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAAt          <  1:213139/71‑1 (MQ=255)
                                                        gCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGCCGTTCAGATACTGAATCa        >  1:612595/1‑71 (MQ=255)
                                                        gCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCa        >  1:1740463/1‑71 (MQ=255)
                                                        gCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCa        >  1:1837782/1‑71 (MQ=255)
                                                             ggAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTcc   <  1:2616492/71‑1 (MQ=255)
                                                              gAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTccc  <  1:2503207/71‑1 (MQ=255)
                                                               aaaaaCACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCa        <  1:742354/64‑1 (MQ=255)
                                                                 aaaCACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGa                  >  1:1892331/1‑52 (MQ=255)
                                                                  |                                                                  
ATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTCAGATACTGAATCATTTCCC  >  W3110S.gb/1538204‑1538336

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: