Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,134,811 Δ1 bp coding (681/693 nt) yghT → putative ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,134,8110T.100.0% 120.7 / NA 30coding (681/693 nt)yghTputative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/14);  total (16/14)

TGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACC  >  NC_000913/3134721‑3134941
                                                                                          |                                                                                                                                  
tgCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGcct                                                                                   >  2:369815/1‑137 (MQ=255)
       cTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                            <  2:506658/139‑1 (MQ=255)
        tGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGa                                                                           <  1:181513/139‑1 (MQ=255)
          gTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAg                                                                         >  1:1004979/1‑139 (MQ=255)
                    aTGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTaatc                                                                                                                      >  2:64905/1‑81 (MQ=255)
                    aTGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATc                                                                                                                      <  1:64905/84‑1 (MQ=255)
                       ccAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                            <  1:225926/139‑1 (MQ=255)
                             aGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGc                                                      >  2:952491/1‑139 (MQ=255)
                                               tggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                       >  2:49940/1‑120 (MQ=255)
                                               tggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                       <  1:49940/120‑1 (MQ=255)
                                                ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTaatc                                                                                                                      <  1:245414/56‑4 (MQ=255)
                                                ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATc                                                                                                                      >  2:245414/1‑56 (MQ=255)
                                                ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                       >  1:531627/1‑119 (MQ=255)
                                                ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                       <  2:531627/119‑1 (MQ=255)
                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                     >  1:898829/1‑68 (MQ=255)
                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                     <  1:1310730/68‑1 (MQ=255)
                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                     <  2:898829/68‑1 (MQ=255)
                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                     >  2:1310730/1‑68 (MQ=255)
                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCagaaag                              >  2:16671/1‑139 (MQ=255)
                                                      tGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCggg                                                                                      >  1:1548531/1‑82 (MQ=255)
                                                      tGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCggg                                                                                      <  2:1548531/82‑1 (MQ=255)
                                                      tGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCagaaaga                             >  2:12460/1‑139 (MQ=255)
                                                                gCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCggg                                                                                      <  2:692378/72‑1 (MQ=255)
                                                                gCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCggg                                                                                      >  1:692378/1‑72 (MQ=255)
                                                                gCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTgccg                   >  2:1192809/1‑139 (MQ=255)
                                                                   aTTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                     >  2:1159792/1‑54 (MQ=255)
                                                                   aTTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                     <  1:1159792/54‑1 (MQ=255)
                                                                                 gACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTTAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATAcc  <  2:1004979/139‑1 (MQ=255)
                                                                                    cgcgcgGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCggg                                                                                      >  2:234476/1‑52 (MQ=255)
                                                                                    cgcgcgGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCggg                                                                                      <  1:234476/52‑1 (MQ=255)
                                                                                          |                                                                                                                                  
TGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACC  >  NC_000913/3134721‑3134941

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