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_000913	3,516,104	0	T	C	100.0%	141.1 / NA	40	K401K (AAA→AAG)	damX	cell division protein that binds to the septal ring
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/24);  total (16/24)

CAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGTCGTAGT  >  NC_000913/3515988‑3516237                                                                                                                     |                                                                                                                                      cAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAg                                                                                                                 >  1:442502/1‑139 (MQ=255)       tAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTaca                                                                                                           <  1:792571/139‑1 (MQ=255)            cATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCa                                                                                                      <  1:628776/139‑1 (MQ=255)            cATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCa                                                                                                      <  1:248192/139‑1 (MQ=255)                gCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGaaa                                                                                                  <  2:757404/139‑1 (MQ=255)                     agaTAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAgg                                                                                             <  1:390716/139‑1 (MQ=255)                                aTTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACgg                                                                                  >  1:540520/1‑139 (MQ=255)                                     ttCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGAc                                                                             >  1:193520/1‑139 (MQ=255)                                     ttCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGAc                                                                             <  2:193520/139‑1 (MQ=255)                                      tCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGAcc                                                                            >  1:391806/1‑139 (MQ=255)                                      tCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGACCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCa                                                                                                      <  2:271340/113‑1 (MQ=255)                                      tCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGACCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCa                                                                                                      >  1:271340/1‑113 (MQ=255)                                              gCCTGTACCTGACGCAGCGGTTTCGCCTACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTAcgcg                                                                    <  1:36047/139‑1 (MQ=255)                                                cTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGcgtc                                                                  <  2:181310/139‑1 (MQ=255)                                                cTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGcgtc                                                                  <  2:507880/139‑1 (MQ=255)                                                  gTACCTGACGCAGCGGTTTCGCCCACGGGCTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCCAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGcgtcgt                                                                <  2:540520/139‑1 (MQ=255)                                                        gaagCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATa                                                          <  1:583845/136‑1 (MQ=255)                                                        gACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATAc                                                                                       >  1:328370/1‑110 (MQ=255)                                                        gACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATAc                                                                                       <  2:328370/110‑1 (MQ=255)                                                         aCGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCa                                                                                     <  1:517700/111‑1 (MQ=255)                                                         aCGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCa                                                                                     >  2:517700/1‑111 (MQ=255)                                                         aCGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTAcgcg                                                                    <  2:756678/128‑1 (MQ=255)                                                         aCGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTAcgcg                                                                    >  1:756678/1‑128 (MQ=255)                                                             aGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACaa                                                     >  2:433589/1‑139 (MQ=255)                                                                  tttCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAg                                                >  1:72513/1‑139 (MQ=255)                                                                     cGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGa                                                                                            <  1:174459/92‑1 (MQ=255)                                                                     cGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGa                                                                                            >  2:174459/1‑92 (MQ=255)                                                                         cACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCa                                         <  2:274177/139‑1 (MQ=255)                                                                           cGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGa                                       >  2:195295/1‑139 (MQ=255)                                                                                        tGGACATCTGCTGGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCg                          >  1:331523/1‑139 (MQ=255)                                                                                               ctgctgGCAATGTAGATACCGCCTTTTTCGCCTCTCCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAAcc                   >  2:255686/1‑139 (MQ=255)                                                                                                  ctgGCAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGtt                <  1:195295/139‑1 (MQ=255)                                                                                                      cAATGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAgg            <  1:433589/139‑1 (MQ=255)                                                                                                         tGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCaa                     <  2:21097/127‑1 (MQ=255)                                                                                                         tGTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCaa                     >  1:21097/1‑127 (MQ=255)                                                                                                          gTAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGt        <  2:331523/139‑1 (MQ=255)                                                                                                           tAGATACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGTc       <  2:72513/139‑1 (MQ=255)                                                                                                              aTACCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGTCGTa    <  1:540588/139‑1 (MQ=255)                                                                                                              aTACCGCCTTTTGCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGTCGTa    <  1:339974/139‑1 (MQ=255)                                                                                                                aCCGCCTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGTCGTAGt  >  1:212625/1‑139 (MQ=255)                                                                                                                     |                                                                                                                                      CAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCCTCTTCTTTCGAAGCGTACACGCCAGAAACCAGGACATACCACGGCTGACCATTACGCGTCGTTTCATAGACAACGTAGTTTTTCAGATTCTCTTTCTTCGCCCAACCGTTCAGGTTGTCGTAGT  >  NC_000913/3515988‑3516237

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