Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,592,188 C→T H178Y (CAT→TAT)  tsr → methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,1880CT100.0% 62.7 / NA 21H178Y (CAT→TAT) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/10);  total (11/10)

CGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAA‑TGAATCGCCTGATT‑GACAGCATTCG‑TCATATTGCAGGCGGCGATCT‑GGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCC  >  NC_000913/4591946‑4592428
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
ggCTGAAACAGGCGGAAAAAAACTGGGCAATTTACGAAGCGTTGCCCCGTTACCCTCGGCAGAGCCCCGCCGCAGCGGCAGAGATCAAACGTAAATACGATATTTATCACAATGCGCTGGCGGATCTGATCCAACTATTAAGTGCAGGCAAAATCAACGAGTTCTTTTATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATc                                                                                                                                                                                                                                               <  1:123198/250‑1 (MQ=255)
    gAAACAGGCGGAAAAAAACTGGGCAGATTACGACGCGTTGCCGCGTGACCCGCGTCAGCGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGCTATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGatctctccg                                                                                                                                                                                                                                           >  2:86401/1‑244 (MQ=255)
      aaCAGGCGGAAAAAAACTGGGCAGATTCCGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTc                                                                                                                                                                                                                                         <  1:126957/251‑1 (MQ=255)
                  aaaccggggCAGATTACGAAGCGTTGCCGCGTGACCCGCGTAGGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAAt                                                                                                                                                                                                                             <  2:121145/245‑1 (MQ=255)
                         ggcagATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGcctcct                                                                                                                                                                                                                      >  1:91444/5‑251 (MQ=255)
                                                                           cGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTGTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCGGGTGGGCGTGATGATCGGCGTACTGGCg                                                                                                                                                                    <  2:257111/251‑1 (MQ=255)
                                                                            ggCAGAGATCAAACGTAATTACGCTATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAACCGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCgg                                                                                                                                                                   >  2:90011/1‑251 (MQ=255)
                                                                                              ttACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGGCTggg                                                                                                                                                >  1:233477/1‑250 (MQ=255)
                                                                                                         aTCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCATCGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTCTGGTTCGGTATTaa                                                                                                                                      <  1:75145/250‑1 (MQ=255)
                                                                                                              aaTGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTCTGGTTCGGTATTAAAGCCTc                                                                                                                                <  1:60450/251‑1 (MQ=255)
                                                                                                                 gcgcTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTCTGGTTCGGTATTAAAGCCTCGc                                                                                                                              <  2:233477/250‑1 (MQ=255)
                                                                                                                                       tattAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTAc                                                                                                                                                                         >  2:221557/3‑186 (MQ=255)
                                                                                                                                       tattAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAGCGTCGTAc                                                                                                                                                                         <  1:221557/184‑1 (MQ=255)
                                                                                                                                        attAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTTTGGTGTGGTTTTTAAGCCTCGCTGGTAGAGCCCA‑TGAAtagactg                                                                                                      >  1:89239/2‑244 (MQ=255)
                                                                                                                                           aGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTCTGGTTCGGTTTTTAAGCCCCGCTGGGTGCCCCAA‑TGAATCGCCTGATT‑g                                                                                                 >  1:134053/1‑251 (MQ=255)
                                                                                                                                            ggTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTCTGGTTTGGTATTAAAGCCTCCCTGGTAGCGCCAA‑TGAATCGCCTGATT‑g                                                                                                 >  1:3653/1‑250 (MQ=255)
                                                                                                                                                                gTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGTCGATGTGGATTCTGGTGGGCGTTATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAA‑TGAATCGCCTGATT‑GACAGCATTCG‑TCATATTGCAg                                                                           <  1:231897/251‑1 (MQ=255)
                                                                                                                                                                                             aGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGCCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCCCC‑TTTGCCGTCTGGGTCGGTATTAAAGCCTCGCTGGTTGCGCCAA‑AGAATAGCCTGATT‑GACAGCATTCG‑TCAAATTTCTGGCGGCGGTCT‑GGGGAAAACGATTGAGGTg                                             >  1:181931/1‑251 (MQ=255)
                                                                                                                                                                                                     ttCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAA‑TGAATCGCCTGATT‑GACAGCATTCG‑TCATATTGCAGGCGGCGATCT‑GGTGAAACCGATTGAGGTGGATGGct                                      <  2:91444/250‑1 (MQ=255)
                                                                                                                                                                                                                             aTGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGTATC‑TTCGCCCGCTGGGTCGGGATTAAAAGCTCGCTGGGAGAGACCA‑TGAATCGCCTGATT‑TGCAGCATTCG‑GCATATTGGAGGCGGCGGTCT‑GGTGGAACCGATTGAGGGGGATGGGTCTAATGGGATTGGGCaaaaggg                >  1:213446/1‑243 (MQ=255)
                                                                                                                                                                                                                                         gATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCCTCCTTC‑CCGGCTTGTTTGGTATTAAAGCCTCGCTGGTAGCACCAAATAAAT‑GCCTGATTTGA‑AGCCTTCGGTC‑TTTTGGAGGCGGCGGTCTTGGTGAAACGGAT‑GGGGTGGATGGGTCTAATGAGAGGGGGCAACTgggcgagaggtttcgcc  >  2:204423/1‑235 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
CGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATC‑TTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAA‑TGAATCGCCTGATT‑GACAGCATTCG‑TCATATTGCAGGCGGCGATCT‑GGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCC  >  NC_000913/4591946‑4592428

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