Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,286,859 T→C intergenic (+333/+207) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,8590TC100.0% 47.0 / NA 18intergenic (+333/+207)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/8);  total (10/8)

GGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAA  >  NC_000913/1286805‑1286989
                                                      |                                                                                                                                  
aatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAg                                                                                                           >  2:490339/3‑80 (MQ=11)
aatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAg                                                                                                           <  1:490339/78‑1 (MQ=11)
                     cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggg                                                                                              >  1:105721/1‑72 (MQ=11)
                     cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggg                                                                                              <  2:105721/72‑1 (MQ=11)
                         catGGGCAGAATATTTGATTGCTGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggg                                                                                              <  2:553807/66‑1 (MQ=11)
                         catGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggg                                                                                              >  1:553807/3‑68 (MQ=255)
                          atGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                            >  2:1047910/2‑69 (MQ=255)
                          atGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                            <  1:1047910/68‑1 (MQ=255)
                          atGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                      >  1:466904/2‑139 (MQ=255)
                            gggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                                                     <  1:199220/106‑1 (MQ=25)
                            gggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                                                     >  2:199220/1‑106 (MQ=25)
                                        tGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                                                                 <  1:1454773/82‑1 (MQ=255)
                                        tGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                                                                 >  2:1454773/1‑82 (MQ=255)
                                            tGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                            >  1:1380103/1‑51 (MQ=255)
                                            tGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                            <  2:1380103/51‑1 (MQ=255)
                                              cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtg                                                                                                   <  2:112121/42‑1 (MQ=255)
                                              cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtg                                                                                                   >  1:112121/1‑42 (MQ=255)
                                              cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGaa  >  2:1144114/1‑139 (MQ=255)
                                                      |                                                                                                                                  
GGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAA  >  NC_000913/1286805‑1286989

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