Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 4,305,151 G→C P51P (CCG→CCC LY180_20645 → rhsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065844,305,1510GC100.0% 54.8 / NA 20P51P (CCG→CCCLY180_20645rhsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (8/12);  total (8/12)

CCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACA  >  CP006584/4304862‑4305312
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                 
ccgcACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCgg                                                                                                                                                           <  1:307367/298‑1 (MQ=255)
        cgcTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGCGCGGAAAACCGGCGGCGCGTCCGGGCGACATGACGCAGTATGGCGGTAGCAGTGTCCCGGTTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGGGGCCTGTTCGGTGTGCCCCGGCGTGGTGACTTCCGGCCCCCCGGCCCATCCCCTGCCCGTGGCaacgt                                                                                                                                                   >  2:428078/1‑295 (MQ=255)
              ttttGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTc                                                                                                                                             <  1:422931/298‑1 (MQ=255)
                                  ttttATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGcc                                                                                                                         >  1:39351/1‑298 (MQ=255)
                                  tttcttaAATTATGAAAATGTCCGGCGTTTTACTTCTGTTTTGTTTCCGGTTGGTGATTGTGGCTTGGCTTCTGGTTCTTATTTGTTTAAAAAGGATTTTACGATGAGCGGAACACCGGCGGCGCGTCAGGGCGCCATGACGACGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGTACGTATGGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGcc                                                                                                                         <  2:191140/293‑1 (MQ=14)
                                          aTTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccgg                                                                                                                 >  1:168970/1‑297 (MQ=255)
                                          aTTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGCTGTTGACTGGCGCTAGGCCTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGCAAAACCGGCGCCGCGTCAGGGCGACATGACCCCGTCTGGCGGTAGCATTGTCCAGGGATCACCCGGGGCCCCTCCTGGTCCCCCCACCGGCGCGGCCTCCTCGCTGTCCCCCGGCGGCGCGCACTCCGGCCCTCCGTCCCACCCCCCGCCCGGGCCAATCGGCCTCCCCGGTGCAACCCACCTCGCCCTgccctg                                                                                                                 >  2:424571/1‑296 (MQ=255)
                                             cTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGcccggccc                                                                                                              <  1:156258/298‑1 (MQ=255)
                                                                  aCTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCg                                                                                         <  1:82780/298‑1 (MQ=255)
                                                                   cTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTTCCAGGGTTCAGCCGGGGTACGTATTGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCCTCCCGTCAACCCCCTGCCCCGTGCAAAGGTCCTTCCCGGTGAAACCGACACCCCCCTGCCCGGCCCCCCGCCGTTCCCCCTCTccccc                                                                                        >  2:376683/1‑296 (MQ=18)
                                                                           cATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTAcag                                                                                <  1:184586/298‑1 (MQ=255)
                                                                                      gttgACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGa                                                                     <  1:428078/298‑1 (MQ=255)
                                                                                               gtgtatgCTTTTGGTGCTTAGTGATTTCAAGAGGTTTTTACGAGGGGGGGAAAACCGGCGGCGCGTCAGGGCGACAGGACGCAGTATGGCGGTAGCAGTGTCCAGGGTTCAGCCGGGGTACGGATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAgcgc                                                            <  2:131292/292‑1 (MQ=11)
                                                                                                          tGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTg                                                 >  1:365065/1‑298 (MQ=255)
                                                                                                              tCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGa                                             >  1:294169/1‑298 (MQ=255)
                                                                                                                ttatttatTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAACATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGCTCGGTGTGCCCCGGCGGAGTGACGCCCGGCCCTCCGGTCAATCCCCTGCGCGGTGCAACGGCCCTTCCCGGGGAAACCGGCACCCCCCTGCCCTGCCCGCCGCCGTTCACCCCCCCCCGCACCTACAGCAGTTCCCGGCACAACCCGCCCGCGCCGGTGGGgcgc                                           >  2:158706/1‑298 (MQ=14)
                                                                                                                    ttatttCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCg                                       <  1:216917/298‑1 (MQ=255)
                                                                                                                         tCAAAAGGATTTTACGATGAGCGGAAACCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcccc                                  <  1:188658/298‑1 (MQ=255)
                                                                                                                                          tGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCTCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGc                 <  1:423809/298‑1 (MQ=14)
                                                                                                                                                         cggcgCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACa  <  1:25328/298‑1 (MQ=11)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                 
CCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACA  >  CP006584/4304862‑4305312

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