Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,684,783 A→G K175K (AAA→AAG tus → DNA replication terminus site‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,7830AG100.0% 105.8 / NA 31K175K (AAA→AAGtusDNA replication terminus site‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/14);  total (17/14)

GAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTG  >  NC_000913/1684649‑1684919
                                                                                                                                      |                                                                                                                                        
gAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCatat                                                                                                                                     >  1:48478/1‑140 (MQ=255)
     aaacccccccccCGCGGCATTTTTTGAATGGGTGCATCGTCTTTTGCGGGGGCTGATCCCCTTTAATGCTTACCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTa                                                                                                                                <  2:114213/134‑1 (MQ=255)
            cccACCGAGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGCTCACCCTTAATGCTTACAGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTCATACGCATATCATTAAGTAttt                                                                                                                         >  2:349015/1‑140 (MQ=255)
                                             cATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaaa                                                                                        >  1:220760/1‑140 (MQ=255)
                                                        gCTGATCCCCTTTAATGCTTCCCGCCCGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAAt                                                                             <  1:88440/140‑1 (MQ=255)
                                                          tGATCCCCTTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATca                                                                           <  2:561426/140‑1 (MQ=255)
                                                               ccccTTAATGTTCCCCGCCCTCTCCCCTTCCGGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACg                                                                      <  2:220760/139‑1 (MQ=255)
                                                                    tAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtg                                                                 >  2:277976/1‑140 (MQ=255)
                                                                           tACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCAcc                                                          >  1:9870/1‑140 (MQ=255)
                                                                                  cGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcg                                                   >  2:1005658/1‑140 (MQ=255)
                                                                                     tCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgc                                                >  2:1075738/1‑140 (MQ=255)
                                                                                           ttCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGATGCGCgaggag                                          >  2:620970/1‑140 (MQ=255)
                                                                                              cgCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTgg                                       <  2:574894/139‑1 (MQ=255)
                                                                                                 aCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCaa                                    >  2:105100/1‑140 (MQ=255)
                                                                                                 aCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCaa                                    >  2:1161852/1‑140 (MQ=255)
                                                                                                     ccccGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGaa                                >  1:434053/1‑140 (MQ=255)
                                                                                                     ccccGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGaa                                >  2:685791/1‑140 (MQ=255)
                                                                                                     ccccGCCACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGaa                                >  2:155244/1‑140 (MQ=255)
                                                                                                          ccACTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTg                           >  1:649901/1‑140 (MQ=255)
                                                                                                            aCTTTACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGa                         <  1:1161852/140‑1 (MQ=255)
                                                                                                               ttACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCg                      <  2:66096/140‑1 (MQ=255)
                                                                                                                tACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAg                        <  1:1063837/137‑1 (MQ=255)
                                                                                                                tACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAg                        >  2:1063837/1‑137 (MQ=255)
                                                                                                                tACGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCga                     <  2:434053/140‑1 (MQ=255)
                                                                                                                 aCGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCgag                    <  1:789166/140‑1 (MQ=255)
                                                                                                                 aCGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCgag                    <  2:649901/140‑1 (MQ=255)
                                                                                                                 aCGCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCgag                    <  1:1005658/140‑1 (MQ=255)
                                                                                                                   gCTTTGGTTGGGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGagag                  <  1:685791/140‑1 (MQ=255)
                                                                                                                        ggttggGCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCa             >  1:185082/1‑140 (MQ=255)
                                                                                                                              gCTAATAAGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGatat       >  2:180151/1‑140 (MQ=255)
                                                                                                                                   taaGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGctg  <  1:698302/140‑1 (MQ=255)
                                                                                                                                   taaGCATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGctg  >  2:893453/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                        
GAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTG  >  NC_000913/1684649‑1684919

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