Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,637,295	C→T	R21H (CGT→CAT)	der ←	GTPase, multicopy suppressor of ftsJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,637,295	0	C	T	100.0%	121.4 / NA	38	R21H (CGT→CAT)	der	GTPase, multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/18);  total (20/18)

TCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCA  >  NC_000913/2637166‑2637394                                                                                                                                  |                                                                                                    tctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATa                                                                                            >  1:794325/1‑139 (MQ=255)   tgtgCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAAc                                                                                          >  1:596193/1‑139 (MQ=255)     tgCCATCAATCCCGCCGGTATCAATAAAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGt                                                                                        >  2:494589/1‑139 (MQ=255)         aTCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                    <  2:444595/139‑1 (MQ=255)         aTCAATCCCGCCGGTATCAATAAAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                    <  1:494589/139‑1 (MQ=255)          tCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAtt                                                                                   >  1:423450/1‑139 (MQ=255)          tCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAtt                                                                                   >  2:323447/1‑139 (MQ=255)             aTCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                <  1:595146/139‑1 (MQ=255)                                aGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGAc                                                             >  2:447548/1‑139 (MQ=255)                                  aTAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaaa                                                                                              >  1:412692/1‑104 (MQ=255)                                  aTAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaaa                                                                                              <  2:412692/104‑1 (MQ=255)                                   tAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAg                                                          >  1:248806/1‑139 (MQ=255)                                       cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                >  1:390612/1‑113 (MQ=255)                                       cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                <  2:390612/113‑1 (MQ=255)                                              ccTTCAAGTTCCGCACGACCGTACTTACGGTCACGAGTAAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAg                                               <  2:669352/139‑1 (MQ=255)                                                ttCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGt                                             >  2:351036/1‑139 (MQ=255)                                                 tCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTa                                            >  2:380158/1‑139 (MQ=255)                                                                   tACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                     >  1:327536/1‑80 (MQ=255)                                                                   tACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                     <  2:327536/80‑1 (MQ=255)                                                                     cTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAttt                        <  2:423450/139‑1 (MQ=255)                                                                          ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                            >  2:162151/1‑98 (MQ=255)                                                                          ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                            <  1:162151/98‑1 (MQ=255)                                                                          ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                   <  2:546761/139‑1 (MQ=255)                                                                             cACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAat                >  1:784063/1‑139 (MQ=255)                                                                              aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                    >  2:766640/1‑102 (MQ=255)                                                                              aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                    <  1:766640/102‑1 (MQ=255)                                                                              aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                   >  2:389034/1‑135 (MQ=255)                                                                              aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                   <  1:389034/135‑1 (MQ=255)                                                                               cGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatca              >  1:569278/1‑139 (MQ=255)                                                                                         ccGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg    <  1:447548/139‑1 (MQ=255)                                                                                           gggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca  <  2:248806/139‑1 (MQ=255)                                                                                           gggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca  <  2:596193/139‑1 (MQ=255)                                                                                            ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTAc                                                                            >  2:371715/1‑64 (MQ=255)                                                                                            ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTAc                                                                            <  1:371715/64‑1 (MQ=255)                                                                                                 tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                               <  2:200500/40‑1 (MQ=255)                                                                                                 tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                               >  1:200500/1‑40 (MQ=255)                                                                                                 tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                      >  2:597914/1‑81 (MQ=255)                                                                                                 tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                      <  1:597914/81‑1 (MQ=255)                                                                                                                                  |                                                                                                    TCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCA  >  NC_000913/2637166‑2637394

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