Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,282,627 C→T Q922* (CAG→TAG)  narG → nitrate reductase A subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,282,6270CT80.0% 42.4 / ‑3.3 20Q922* (CAG→TAG) narGnitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base T (4/12);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGC  >  NC_000913/1282378‑1282873
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
gTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGATTGCCCGGGACAccc                                                                                                                                                                                                                                                        >  1:82960/1‑250 (MQ=255)
     ccagcTGGGTAAAGAAACCGACATCGTCCCTCTGCCTATCCAGCATGCCTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCg                                                                                                                                                                                                                                                  <  2:8476/247‑1 (MQ=255)
                     aCCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCg                                                                                                                                                                                                                                  <  2:191852/251‑1 (MQ=255)
                         aCATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAg                                                                                                                                                                                                                              <  1:164300/251‑1 (MQ=255)
                                                        tGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGGCTGGAAAACCCAGAGGGAGGAGGAGCTGCTGCGTAAGCTCAACTACAACAAAACGGAAGGGCCGGCgg                                                                                                                                                                                               >  1:196346/1‑250 (MQ=255)
                                                                       gcAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGAATGCCTGGGAAACCCAGAGCGAGGAGGGTCTGCTTCGGAAACTCAACTACACCAAAGCGGAAAGGCCGGCGAAAAGCCAGCCGacgc                                                                                                                                                                                >  2:293965/1‑248 (MQ=255)
                                                                                                 aaaaaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTAGCCGGCGACTTACGAACGCGTTACCACTAGCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAAGTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCagag                                                                                                                                                      <  2:82960/251‑1 (MQ=255)
                                                                                                             gCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGac                                                                                                                                           <  1:3908/250‑1 (MQ=255)
                                                                                                                    gATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCAc                                                                                                                                   <  1:292043/251‑1 (MQ=255)
                                                                                                                                aaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACg                                                                                                                       <  1:206477/251‑1 (MQ=255)
                                                                                                                                           acaTTCTGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGCTGGAGAAAATCGGTAATGGCGGTAAAGGGATTGACTGGGAAACCCAAAGCGAGATGGATCTGCTGCGTAAGCTCCACTACACCCAAAAGGGAGGGCCGGCGGAAGGCCAGCACCTGCTGAACACCCCACATGATGCGGGAGAGACGATACCGGCACCGGCCCCGGAAACCaaaaggccggtagcc                                                                                                            >  2:47995/1‑238 (MQ=255)
                                                                                                                                                        agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTGCCTGGGACACCTAGAGCGAGATGGATCTGCTTCGTAAGATCAACTACACCAAAAAGGAAGGGCCGGCGAAAGGCCAGCCGATGCTGAACACCCCAATTGATGCGGCAGAGATGATCCTGACCCTGGCAACGGAAACCAACGGTCAGGGAGCCGTGGAAGCCTggg                                                                                               >  1:269549/1‑251 (MQ=255)
                                                                                                                                                                               aCGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGt                                                                        <  2:217778/251‑1 (MQ=255)
                                                                                                                                                                                cGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCAGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTc                                                                       <  2:212464/251‑1 (MQ=255)
                                                                                                                                                                                    cGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGa                                                                   <  1:293965/251‑1 (MQ=255)
                                                                                                                                                                                          aCCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATAc                                                             <  2:262013/251‑1 (MQ=255)
                                                                                                                                                                                                        gCTGATGGCGACAATCGGTAATGGCGGTCAAAGGAGTTGGTGGAAAACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAACCGATGCCGAATACCGCAAGTGTTGCGGGAGAGATGATCCTGACACTTGCAACGGAAACCAAAAGTCAAGTAGCCGTGAAAGACTGGGGTGCCCTGAGATAATTTAACGGTCGGGAACATACgaattcggcggca                                                >  1:226414/1‑238 (MQ=255)
                                                                                                                                                                                                                         gTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATGTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGagaag                              <  1:47995/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGAGCTGAATAAAGAAGACGAGAAGATCCg                         >  2:262187/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aaGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGAAACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCcgc                  >  2:83032/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      cacCTAGACCGAGCTGGATCTGCTGCGTCAGCTAAACTACCCCAAAGCGGAAGGTCAGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGc  >  2:152564/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
GTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGC  >  NC_000913/1282378‑1282873

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