Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,516,104 T→C K401K (AAA→AAG damX ← cell division protein that binds to the septal ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,516,1040TC100.0% 108.1 / NA 31K401K (AAA→AAGdamXcell division protein that binds to the septal ring
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/17);  total (14/17)

CAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGTCG  >  NC_000913/3515975‑3516233
                                                                                                                                 |                                                                                                                                 
cAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGcc                                                                                                                          <  1:134928/139‑1 (MQ=255)
  ggCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCctc                                                                                                                        <  2:373079/139‑1 (MQ=255)
                          ttGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGCCGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTCGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGa                                                                                                >  2:55926/1‑139 (MQ=255)
                             cGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAAc                                                                                             <  2:693237/139‑1 (MQ=255)
                             cGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAAc                                                                                             >  1:697839/1‑139 (MQ=255)
                                  gaTAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTa                                                                                                         >  1:193910/1‑122 (MQ=255)
                                  gaTAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTa                                                                                                         <  2:193910/122‑1 (MQ=255)
                                            aTTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACgg                                                                              >  2:744925/1‑139 (MQ=255)
                                            aTTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACgg                                                                              <  1:213523/139‑1 (MQ=255)
                                                 ttCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGAc                                                                         >  2:797236/1‑139 (MQ=255)
                                                      aTCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTa                                                                    >  2:319491/1‑139 (MQ=255)
                                                            cTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGcgtc                                                              <  1:797236/139‑1 (MQ=255)
                                                                    gACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCa                                                                                                  <  2:2212/95‑1 (MQ=255)
                                                                    gACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCa                                                                                                  >  1:2212/1‑95 (MQ=255)
                                                                     aCGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTc                                                         >  1:705630/1‑135 (MQ=255)
                                                                     aCGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTc                                                         <  2:705630/135‑1 (MQ=255)
                                                                      cGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTTCGCGTCGTTTCATAGa                                                    <  2:469034/139‑1 (MQ=255)
                                                                         aGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGaa                                                                                                              >  2:174749/1‑78 (MQ=255)
                                                                         aGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGaa                                                                                                              <  1:174749/78‑1 (MQ=255)
                                                                         aGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACaa                                                 >  1:600119/1‑139 (MQ=255)
                                                                                   cccACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGa                                                                                        <  2:440536/90‑1 (MQ=255)
                                                                                   cccACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGa                                                                                        >  1:440536/1‑90 (MQ=255)
                                                                                       cGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGa                                   >  1:194093/1‑139 (MQ=255)
                                                                                               tGGCATGAACATGTGCTGGCAATGTAGATACCGCCTGTTTCGCCTCTTCTTTCGAAGCGTACACGGCAGAAACCAGGACAGACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCtcttt                           <  1:94534/139‑1 (MQ=255)
                                                                                                  ccTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCtctttctt                        >  1:689878/1‑139 (MQ=255)
                                                                                                         aTCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCaa                 <  2:600119/139‑1 (MQ=255)
                                                                                                                    aTGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGtt      <  2:689878/139‑1 (MQ=255)
                                                                                                                      gTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGt    <  1:319491/139‑1 (MQ=255)
                                                                                                                        aGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGTCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGTCg  <  2:697839/139‑1 (MQ=255)
                                                                                                                          aTACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAg                                                     <  2:523401/86‑1 (MQ=255)
                                                                                                                          aTACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAg                                                     >  1:523401/1‑86 (MQ=255)
                                                                                                                                 |                                                                                                                                 
CAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGTCG  >  NC_000913/3515975‑3516233

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