Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	977,322	G→A	A333T (GCC→ACC)	mukB →	chromosome condensin MukBEF, ATPase and DNA‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	977,322	0	G	A	100.0%	129.5 / NA	33	A333T (GCC→ACC)	mukB	chromosome condensin MukBEF, ATPase and DNA‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/18);  total (15/18)

CCATCTCGACAAAGCCCTGGAGTTTCGTCGCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGGCCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGCAAAATGAAGTGGTGGCAGAAGCCATCGAACGCCAGCAAGAGAATG  >  NC_000913/977174‑977469                                                                                                                                                     |                                                                                                                                                    ccATCTCGACAAAGCCCTGGAGTTTCGTCGCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGAcc                                                                                                                                                   <  2:207390/151‑1 (MQ=255)           aaaGCCCTGGAGTTTCGTCGCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCAcc                                                                                                                                         >  2:196942/1‑151 (MQ=255)            aaGCCCTGGAGTTTCGTCGCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCt                                                                                                                                        <  2:184922/151‑1 (MQ=255)               cccTGGAGTTTCGTCGCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGa                                                                                                                                      >  2:337118/1‑150 (MQ=255)               cccTGGAGTTTCGTCGCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGa                                                                                                                                      >  2:250138/1‑150 (MQ=255)                  tGGAGTTTCGTCGCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCt                                                                                                                                  >  2:401029/1‑151 (MQ=255)                          cgtcgCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCaaa                                                                                                                          <  2:105944/151‑1 (MQ=255)                          cgtcgCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCaa                                                                                                                           <  2:330458/150‑1 (MQ=255)                          cgtcgCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCa                                                                                                                            >  1:79170/1‑149 (MQ=255)                             cgcgAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAAcc                                                                                                                        <  2:675436/150‑1 (MQ=255)                                   cTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGc                                                                                                                 >  2:684066/1‑151 (MQ=255)                                     aCATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGt                                                                                                               >  1:520649/1‑151 (MQ=255)                                       aTACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTc                                                                                                              <  1:141771/150‑1 (MQ=255)                                         aCTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTcagc                                                                                                           >  1:83465/1‑151 (MQ=255)                                            tCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAgg                                                                                                        >  1:261263/1‑151 (MQ=255)                                                  cAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGgaga                                                                                                     >  1:48090/1‑148 (MQ=255)                                                  cAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGagaag                                                                                                   <  2:347848/150‑1 (MQ=255)                                                      aaCTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCg                                                                                               >  2:139962/1‑150 (MQ=255)                                                      aaCTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCg                                                                                               >  1:330434/1‑150 (MQ=255)                                                       aCTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGaa                                                                                             <  2:64303/151‑1 (MQ=255)                                                                  aGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGc                                                                                  <  2:192648/151‑1 (MQ=255)                                                                            acacGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGAt                                                                        <  1:171524/151‑1 (MQ=255)                                                                              acGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGa                                                                      <  1:585920/151‑1 (MQ=255)                                                                                   gATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGc                                                                 <  2:486702/151‑1 (MQ=255)                                                                                     tatGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGc                                                                 >  2:600382/1‑149 (MQ=255)                                                                                        ggCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGAt                                                             <  1:291614/150‑1 (MQ=255)                                                                                                     gCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGc                                               <  1:688275/151‑1 (MQ=255)                                                                                                            acaACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGCAAAATg                                         >  2:317558/1‑150 (MQ=255)                                                                                                                     ccGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGCAAAATGAAGTGGTGGc                               >  2:253314/1‑151 (MQ=255)                                                                                                                     ccGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGCAAAATGAAGTGGTGGc                               <  2:517854/151‑1 (MQ=255)                                                                                                                      cGAAGGCGATCTGGAAGCGGATTATCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGCAAAATGAAGTGGTGGCa                              <  1:34518/151‑1 (MQ=255)                                                                                                                                             aTCAGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGCAAAATGAAGTGGTGGCAGAAGCCATCGAACGCCAGCAag        <  2:11540/150‑1 (MQ=255)                                                                                                                                                aGGCGACCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGCAAAATGAAGTGGTGGCAGAAGCCATCGAACGCCAGCAAGAGaa    <  1:636161/151‑1 (MQ=255)                                                                                                                                                  gcgacCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGCAAAATGAAGTGGTGGCAGAAGCCATCGAACGCCAGCAAGAGAATg  <  2:295672/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCATCTCGACAAAGCCCTGGAGTTTCGTCGCGAGCTACATACTTCGCGTCAGCAACTGGCGGCTGAGCAGTACAAACACGTCGATATGGCGCGTGAGCTGGCAGAGCACAACGGTGCCGAAGGCGATCTGGAAGCGGATTATCAGGCGGCCAGTGATCACCTGAACCTGGTGCAAACCGCACTGCGTCAGCAGGAGAAGATCGAACGCTACGAAGCGGATCTCGATGAGCTGCAGATCCGTCTGGAAGAGCAAAATGAAGTGGTGGCAGAAGCCATCGAACGCCAGCAAGAGAATG  >  NC_000913/977174‑977469

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