Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,281,830 T→C L656P (CTG→CCG)  narG → nitrate reductase A subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,281,8300TC100.0% 49.8 / NA 16L656P (CTG→CCG) narGnitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/5);  total (11/5)

TAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAA  >  NC_000913/1281702‑1281957
                                                                                                                                |                                                                                                                               
taACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGccgccg                                                                                                                             >  1:376295/1‑133 (MQ=255)
       tCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATAGGCTGGCCGCCGTCTGCAc                                                                                                                      >  2:83309/1‑133 (MQ=255)
       tCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTCCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCAc                                                                                                                      >  2:23605/1‑133 (MQ=255)
         cagccagTGGTGTTACGAAACCGTGACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCg                                                                                                                    <  2:182379/133‑1 (MQ=255)
             cagTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGACGGCGTCTGCACCGCAGt                                                                                                                >  1:206042/1‑133 (MQ=255)
                tGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGGCACTTGATCGACTTTAACGTCCGTGCGGAACGCCGGGGCTGGCCGCCGTCTGCACCGCAGCTAg                                                                                                             >  2:273427/1‑133 (MQ=255)
                tGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAg                                                                                                             >  1:276117/1‑133 (MQ=255)
                tGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAg                                                                                                             >  2:7660/1‑133 (MQ=255)
                                           gagTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCg                                                                                  <  1:318439/133‑1 (MQ=255)
                                                                     aTCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCg                                                        <  1:435188/133‑1 (MQ=255)
                                                                     aTCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCg                                                        <  2:376295/133‑1 (MQ=255)
                                                                          gCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGa                                                   <  1:23605/133‑1 (MQ=255)
                                                                                        ttGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAAt                                     >  2:484409/1‑133 (MQ=255)
                                                                                                    aaCGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATGCGGTGAAATCCCTGAAAGAgg                         >  2:277722/1‑133 (MQ=255)
                                                                                                    aaCGTCCGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGAGATGAATCCGGTCGACCATACGGTGTAATCCCTGAAAGAgg                         >  2:189470/1‑133 (MQ=255)
                                                                                                          cGTGCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGGCTACACGGTGAAATCCATGAACGAGGGTTCCa                   >  2:146691/1‑133 (MQ=255)
                                                                                                             gCGGAACGCATGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATTCGGTGAATTCCATGTAAGAGGCTTCCATcc                >  1:240562/1‑133 (MQ=255)
                                                                                                                      aTGGGCTGGCCGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAACCCCTGAACGAGAGTTCCATCCGTTTTGCgg       >  1:201546/1‑133 (MQ=255)
                                                                                                                           ctggccgCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGaa  >  2:110223/1‑133 (MQ=255)
                                                                                                                                |                                                                                                                               
TAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAATCCCGCTATACCGGACACTTGATCGACTTTAACGTCCGTGCGGAACGCATGGGCTGGCTGCCGTCTGCACCGCAGTTAGGCACTAACCCGCTGACTATCGCTGGCGAAGCGGAAAAAGCCGGGATGAATCCGGTGGACTATACGGTGAAATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAA  >  NC_000913/1281702‑1281957

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