Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,282,409 A→G 81.3% P323P (CCA→CCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,282,4090AG81.3% 71.8 / 10.6 32P323P (CCA→CCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (12/14);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCT  >  W3110S.gb/1282342‑1282459
                                                                   |                                                  
gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGc                                                 >  1:2530441/1‑71 (MQ=255)
   ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAg                                              <  1:257469/71‑1 (MQ=255)
     gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTa                                            >  1:1578355/1‑71 (MQ=255)
     gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTa                                            >  1:1727316/1‑71 (MQ=255)
         cGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTAttt                                         >  1:388678/1‑70 (MQ=255)
         cGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTAttt                                         >  1:2690660/1‑70 (MQ=255)
           gCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCAc                                      >  1:220700/1‑71 (MQ=255)
             aaTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCACCTCGACAACCCGAGCCAGTAt                                           >  1:963003/1‑63 (MQ=255)
              aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGa                                   >  1:2352732/1‑71 (MQ=255)
               tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                  >  1:1168155/1‑71 (MQ=255)
                ggCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                  >  1:1933777/1‑70 (MQ=255)
                 gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTa                                <  1:2502220/71‑1 (MQ=255)
                  cgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTa                                <  1:2521489/70‑1 (MQ=255)
                     tGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTAtgtg                            >  1:1008072/1‑71 (MQ=255)
                           tGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGc                      <  1:1970150/71‑1 (MQ=255)
                              gCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTac                   <  1:22926/71‑1 (MQ=255)
                              gCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTac                   <  1:1645469/71‑1 (MQ=255)
                                  cGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACt                                 <  1:1943381/53‑1 (MQ=255)
                                   gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGc                      <  1:18980/63‑1 (MQ=255)
                                    tAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                  >  1:2641207/1‑50 (MQ=255)
                                     aaTGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCAACGACTATGTGCGt                         >  1:2196635/1‑58 (MQ=255)
                                         ctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGcc        <  1:2749475/71‑1 (MQ=255)
                                         ctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGcc        <  1:2752755/71‑1 (MQ=255)
                                         ctgcGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGcc        >  1:2577528/1‑71 (MQ=255)
                                            cGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGAt     <  1:1653318/71‑1 (MQ=255)
                                                aaTTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGcc        <  1:1117055/64‑1 (MQ=255)
                                                aaTTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCt  <  1:1649666/70‑1 (MQ=255)
                                                  ttCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa              <  1:2733049/56‑1 (MQ=255)
                                                  ttCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCt  <  1:2241778/68‑1 (MQ=255)
                                                  ttCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCt  <  1:522669/68‑1 (MQ=255)
                                                      aCCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGcc        <  1:1193460/58‑1 (MQ=255)
                                                       ccTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCg       >  1:305649/1‑58 (MQ=255)
                                                                   |                                                  
GCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCT  >  W3110S.gb/1282342‑1282459

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: