Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,541,037 A→G I518I (ATT→ATC)  narZ ← nitrate reductase Z subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,541,0370AG100.0% 27.2 / NA 11I518I (ATT→ATC) narZnitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (2/9);  total (2/9)

TGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGC  >  NC_000913/1540834‑1541260
                                                                                                                                                                                                           |                                                                                                                                                                                                                               
tgtgGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGATCGCCCATGCTTTTTATGGGCGGTATCGGCAAATTCACGGTCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTgggg                                                                                                                                                                                   <  1:345974/250‑1 (MQ=255)
                                                  ttGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCAGCGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTAACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGcc                                                                                                                                >  1:173051/1‑251 (MQ=255)
                                                                         gcAGCATAGTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACTGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGGGGCCGTAATTCGcc                                                                                                         <  1:45375/250‑1 (MQ=255)
                                                                          cAGCATATTGATCATCCCGCGGTAGTTCATGTACATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAAGTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGcc                                                                                                         <  1:27019/250‑1 (MQ=255)
                                                                              atatTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCGCTCCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTAc                                                                                                    <  1:174838/251‑1 (MQ=255)
                                                                                                        gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCTTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGan                                                                          <  1:92763/251‑2 (MQ=255)
                                                                                                                          nttAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACg                                                                                                                                                                                                               >  1:210092/2‑100 (MQ=255)
                                                                                                                          gTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACg                                                                                                                                                                                                               <  2:210092/100‑1 (MQ=255)
                                                                                                                                          agggtAATCATAGAGCGCCCATGCGGTGTATGGGCGGGAGCGGCAAATTCACGGGCGATGGTTTCGATAGAGGGGCGCGGAACGCCAGTAATTTGCGCACCTCAGGATGGGGGGTACGGTTGGATTTCAGCGTAAGCGTTCGCAAGGTTTTCATCTGCCAGCCCGCGATGGAGGCCGTAATGCGCCAGTGCGAAATCATAAACGCTGACCACCGGACAGGTGGTGCCATCGACTAACGTCAGGTTTTTAAc                                        <  1:207482/247‑1 (MQ=255)
                                                                                                                                                                          gggggTATCGGCAAATTCACGGGCGATGGTTTCGATATTCTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTGCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTc                                                                                                            <  1:166300/148‑1 (MQ=255)
                                                                                                                                                                                aTCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTTTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGTCGCACCAGCACCGGATTGTGTTTTACGc  <  1:301841/251‑1 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                               
TGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGC  >  NC_000913/1540834‑1541260

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