Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,540,902 G→A F563F (TTC→TTT)  narZ ← nitrate reductase Z subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,540,9020GA92.3% 26.8 / ‑3.5 13coding (1689/3741 nt)narZnitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (4/8);  minor base . (1/0);  total (5/8)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

TACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCG‑CAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGG  >  NC_000913/1540682‑1541145
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
tacgggggtCGGCGAGCGGTTAAATTATCCTCTGGGCAGAGACTTTTTCATAGCGGCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTAGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGGGGCAACCAGCCGGTTTGTGGGGGCAGTTTTTCCTGGCCGGCATAGGGTGCCCAGCCGCCACCTCTTTGCCCGACACAACCG‑CAAAAGATCAGCATATTTATCATCCCGCgg                                                                                                                                                                                                                          <  1:251037/240‑1 (MQ=255)
 actGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCCTTTGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACAAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACCACCG‑CAAAAGATCAGCATATTGATCATCCCGCGGTa                                                                                                                                                                                                                        >  1:114260/1‑249 (MQ=255)
                                 gAGCATAGTCTTTTTCATAGCGCCATTTGCTGGAATGATTGTAGAAAAACGAGTTGCTGTTCATTTGGCGCTGTGTTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCATTCAGCCTGTTTGTGGGCGCAGTTTTTCCGGGCCGACATAGTGGGCCCATCCTCCACTGCTTTGCCCGACACTACCT‑CAAAAGATCATCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAacac                                                                                                                                                                                        <  1:250909/249‑1 (MQ=255)
                                                                                                                  tCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGCCACAACCG‑CAAAAGATCAGCATATTGATCATCCCGCGGAAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTgg                                                                                                     >  1:5904/1‑251 (MQ=255)
                                                                                                                   ccAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCG‑CAAAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCTCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGgc                                                                                                    <  1:39868/251‑1 (MQ=255)
                                                                                                                                          gCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGCCATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCG‑CAACAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCg                                                                                                                                      >  1:45628/1‑194 (MQ=255)
                                                                                                                                           cAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCG‑CAAAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGcc                                                                                                                                                                                 <  1:118842/150‑1 (MQ=255)
                                                                                                                                                  ncGGTTTGTGGGCGCAGTTTTTCCTGGCCGAAATAGTGTCCTCAGCCGCCACCGCTTTGCACGCCACCACCGACA‑AAGATCAGTATATTGATCCTCCCGCGGTCGTTCCTGTCCATGTGATCCCAGTGGATAACACCTGCGCCGAGGCTAatca                                                                                                                                                                      >  1:27881/2‑154 (MQ=255)
                                                                                                                                                                  gTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCG‑CAAAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATC‑AGCGCCCATGCGTTTTATGGGCGGTATCgg                                                                                                                                     <  1:71610/170‑1 (MQ=255)
                                                                                                                                                                                          cccAGNCGCCANCNNTTTGCCCGACTCAACCG‑CANAAGATCAGCATATTGNTCATNCCGCGNTANtn                                                                                                                                                                                                                     >  2:96990/1‑66 (MQ=255)
                                                                                                                                                                                          cccAGCCGCCACCGCTTTGCCCGACACAACCG‑CAAAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCGGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTc                                                                                                               <  1:96990/169‑1 (MQ=255)
                                                                                                                                                                                                          nnccccGANNCAACCG‑CCAAAGATCAGCATATTGATCATCCCGCTGTAGTTCATGTCCATGTGATACCAGTTGATAACACCTGCGCCGAGGATAatca                                                                                                                                                                      <  2:307/95‑1 (MQ=255)
                                                                                                                                                                                                            gCCCGACACAACCG‑CAAAAGATCAGCATATTGATCATACCGCGGTAGTTCATGCCCATGAGATCCCAATGGCTAACACCTGCGCCGAGGATAATCATAGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTACTTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCcg           >  1:166946/1‑251 (MQ=255)
                                                                                                                                                                                                                     aaCCG‑CAAAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAng  <  1:198011/251‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
TACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCG‑CAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGG  >  NC_000913/1540682‑1541145

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑