Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,282,409 A→G 15.2% 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,4090AG15.2% 107.3 / 13.4 46P323P (CCA→CCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (28/11);  new base G (5/2);  total (33/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

ACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGG  >  W3110S.gb/1282344‑1282470
                                                                 |                                                             
aCCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                                                            <  1:3290069/69‑1 (MQ=255)
      gCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt                                                      <  1:961881/69‑1 (MQ=255)
      gCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt                                                      <  1:686438/69‑1 (MQ=255)
         gCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTc                                                   >  1:2658983/1‑69 (MQ=255)
         gCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTc                                                   >  1:2362842/1‑69 (MQ=255)
              ggCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTc                                                   >  1:716389/1‑64 (MQ=255)
               gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGGCAACCCAAGCCAGTATTTCACCGAc                                             >  1:1716144/1‑69 (MQ=255)
               gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGAc                                             >  1:1081175/1‑69 (MQ=255)
               gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGAc                                             >  1:358848/1‑69 (MQ=255)
                           ggCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCg                                  >  1:1659001/1‑68 (MQ=255)
                           ggCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGc                                 >  1:1230434/1‑69 (MQ=255)
                           ggCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGc                                 >  1:2280599/1‑69 (MQ=255)
                           ggCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                                                            >  1:1630083/1‑42 (MQ=255)
                           ggCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                                                            >  1:3532629/1‑42 (MQ=255)
                           ggCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                                                            >  1:2997328/1‑42 (MQ=255)
                           ggCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                                                            >  1:1173750/1‑42 (MQ=255)
                             ccACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTa                               >  1:39311/1‑69 (MQ=255)
                             ccACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTa                               >  1:707095/1‑69 (MQ=255)
                             ccACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTa                               >  1:1273188/1‑69 (MQ=255)
                                  tAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACAcc                           >  1:674754/1‑68 (MQ=255)
                                  tAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACAcc                           >  1:2546879/1‑68 (MQ=255)
                                   aaTGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                         <  1:2399500/69‑1 (MQ=255)
                                   aaTGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                         <  1:92219/69‑1 (MQ=255)
                                      gctgcGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAGCGACTATGTGCGTCGCTACACCGACa                       <  1:1105164/68‑1 (MQ=255)
                                      gctgcGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACa                       <  1:631374/68‑1 (MQ=255)
                                      gctgcGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACa                       <  1:2700384/68‑1 (MQ=255)
                                       ctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATg                     >  1:1862649/1‑69 (MQ=255)
                                        tgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                         <  1:1794888/64‑1 (MQ=255)
                                         gcGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGcc                   <  1:2559911/69‑1 (MQ=255)
                                             gtaTTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                         >  1:380395/3‑59 (MQ=255)
                                             gAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                         >  1:1276778/1‑59 (MQ=255)
                                             gAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                         >  1:1553264/1‑59 (MQ=255)
                                              aaTTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTa                                           >  1:1466791/1‑40 (MQ=255)
                                              aaTTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAtg                                         >  1:2797995/1‑42 (MQ=255)
                                              aaTTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAtg                                         >  1:1931463/1‑42 (MQ=255)
                                              aaTTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAtg                                         >  1:2435311/1‑42 (MQ=255)
                                              aaTTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAtg                                         >  1:2577482/1‑42 (MQ=255)
                                                  ccACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGt          >  1:1506408/1‑69 (MQ=255)
                                                  ccACCTCGACAACCCAAGCCAGAATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGt          >  1:713416/1‑69 (MQ=255)
                                                   cACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCt                                >  1:2545060/1‑46 (MQ=255)
                                                   cACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTac                              >  1:1631674/1‑48 (MQ=255)
                                                   cACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTac                              <  1:2088276/48‑1 (MQ=255)
                                                   cACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTac                              <  1:3045252/48‑1 (MQ=255)
                                                   cACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTac                              >  1:2813609/1‑48 (MQ=255)
                                                         gACAACCCAAGCCAGTATTTCACCGACTATTTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCt    >  1:867587/1‑68 (MQ=255)
                                                           cAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTgg  <  1:1668880/68‑1 (MQ=255)
                                                                 |                                                             
ACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGG  >  W3110S.gb/1282344‑1282470

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: