Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	416,770	C→A	T61T (ACG→ACT)	sbcD ←	ATP dependent, structure specific DNA nuclease ‑ SbcD subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	416,770	0	C	A	100.0%	175.1 / NA	52	T61T (ACG→ACT)	sbcD	ATP dependent, structure specific DNA nuclease ‑ SbcD subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (36/16);  total (36/16)

GGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGT  >  NC_000913/416635‑416906                                                                                                                                        |                                                                                                                                         ggCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGt                                                                                                                                        >  2:769110/1‑138 (MQ=255) ggCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGt                                                                                                                                        <  1:769110/138‑1 (MQ=255)  gCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCg                                                                                                                                     >  2:1045411/1‑140 (MQ=255)  gCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCg                                                                                                                                     >  1:926165/1‑140 (MQ=255)                 ttGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGggcgg                                                                                                                      >  1:890304/1‑140 (MQ=255)                    aGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGa                                                                                                                   >  1:639344/1‑140 (MQ=255)                    aGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGa                                                                                                                   >  2:218840/1‑140 (MQ=255)                    aGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGa                                                                                                                   >  2:294575/1‑140 (MQ=255)                    aGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGa                                                                                                                   >  2:583809/1‑140 (MQ=255)                      gAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGc                                                                                                                 >  1:1244623/1‑140 (MQ=255)                      gAACGCCAGGATCGCGCCGGATTCCTTCAGCGTGGCGCCGGGGTCAGGTTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTATCTGGGCGGCGAGc                                                                                                                 <  2:653453/140‑1 (MQ=255)                         cGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCgg                                                                                                              >  2:61972/1‑140 (MQ=255)                         cGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCgg                                                                                                              >  1:23949/1‑140 (MQ=255)                         cGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCgg                                                                                                              >  2:918915/1‑140 (MQ=255)                               gATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGa                                                                                                        >  2:1100128/1‑140 (MQ=255)                                      cgcgATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATc                                                                                                          >  1:41925/1‑131 (MQ=255)                                      cgcgATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATc                                                                                                          <  2:41925/131‑1 (MQ=255)                                      cgcgATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATc                                                                                                 >  2:798905/1‑140 (MQ=255)                                      cgcgATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATc                                                                                                 >  2:154397/1‑140 (MQ=255)                                        cgATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATTAc                                                                                               >  1:668195/1‑140 (MQ=255)                                            tcattcaGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGc                                                                                           <  1:294575/140‑1 (MQ=255)                                            tcattcaGCGTGGCGACCAAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGc                                                                                           <  1:218840/140‑1 (MQ=255)                                             cattcaGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGca                                                                                          <  1:203116/140‑1 (MQ=255)                                                 gaGCGTGGCGACCGGGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGAGAACAAAACGGTTGTATAAAGTGCGGGCTTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGcaacaa                                                                                      <  2:117151/139‑1 (MQ=255)                                                   gCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACaata                                                                                    >  1:673934/1‑140 (MQ=255)                                                    cGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACaataa                                                                                   >  1:894959/1‑140 (MQ=255)                                                    cGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGGTGTATAAAGGGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACaataa                                                                                   >  2:623569/1‑140 (MQ=255)                                                     gTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAAc                                                                                                    <  1:20816/122‑1 (MQ=255)                                                     gTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAAc                                                                                                    >  2:20816/1‑122 (MQ=255)                                                       ggCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCg                                                                                >  2:398844/1‑140 (MQ=255)                                                           aCCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATc                                                                            <  1:955684/140‑1 (MQ=255)                                                               aGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCAcc                                                                        >  2:1058245/1‑140 (MQ=255)                                                                          ccTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTtg                                                             >  2:618563/1‑140 (MQ=255)                                                                            cgCCAGTACCACCAGATGACAGACTGTTTGCCGTAAATTAGCAACAAACCGGTGGTAAAAAGTGCGGGCTTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTtgtg                                                           <  2:668195/139‑1 (MQ=255)                                                                                            tGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCa                                           >  2:702388/1‑140 (MQ=255)                                                                                            tGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCa                                           >  2:494675/1‑140 (MQ=255)                                                                                            tGACAGCCAGTTTGCTGTAAATTCAGAACAAAACGTTTGTATAAAGTGCGGGCGTAACTGGGGGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCa                                           <  2:926165/140‑1 (MQ=255)                                                                                             gACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGGATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAg                                          >  1:1162841/1‑140 (MQ=255)                                                                                              aCAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGt                                         <  1:918915/140‑1 (MQ=255)                                                                                              aCAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGt                                         >  1:647889/1‑140 (MQ=255)                                                                                                     gTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGaaa                                  >  2:749814/1‑140 (MQ=255)                                                                                                     gTTTGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGaaa                                  >  1:510287/1‑140 (MQ=255)                                                                                                        tGCTGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGc                               >  1:1160428/1‑140 (MQ=255)                                                                                                           tGTAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTg                            <  2:23949/140‑1 (MQ=255)                                                                                                             tAAATTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGAt                          >  2:98621/1‑140 (MQ=255)                                                                                                                aTTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGtt                       <  2:375701/140‑1 (MQ=255)                                                                                                                aTTGACAACAAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGtt                       <  1:798905/140‑1 (MQ=255)                                                                                                                    acaacaAAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGc                   >  2:885929/1‑140 (MQ=255)                                                                                                                       acaaAACGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTc                >  2:422814/1‑140 (MQ=255)                                                                                                                          aaaCGGTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTcgcg             <  1:154397/140‑1 (MQ=255)                                                                                                                               gTTGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCttt        >  1:591881/1‑140 (MQ=255)                                                                                                                                 tGTATAAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTa      <  2:1300599/140‑1 (MQ=255)                                                                                                                                  ttataAAGTGCGGGCGTAACGGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTAc     <  2:639344/139‑1 (MQ=255)                                                                                                                                   tataAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACt    >  2:822487/1‑140 (MQ=255)                                                                                                                                     taAAGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGt  >  2:1224007/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                         GGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGT  >  NC_000913/416635‑416906

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