Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,912,914	G→A	E43K (GAG→AAG)	spuI →	glutamylpolyamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,912,914	0	G	A	100.0%	107.8 / NA	34	E43K (GAG→AAG)	spuI	glutamylpolyamine synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/16);  total (18/16)

AACAACAAAGGTGTTGCATCATGTCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTATCCAAT  >  NC_002947/5912768‑5913056                                                                                                                                                   |                                                                                                                                               aacaacAAAGGGGTTGCATCATGTCGGTACCCCCGCGTCACCTTCATCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAgc                                                                                                                                             >  2:681208/1‑150 (MQ=255)     acaaAGGTGTTGCATCATGTCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCAc                                                                                                                                         >  1:589824/1‑150 (MQ=255)                 catcatGTCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACaa                                                                                                                             >  1:712359/1‑150 (MQ=255)                          ggTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACga                                                                                                                    <  1:649279/150‑1 (MQ=255)                              cccccGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGaaa                                                                                                                <  2:712359/150‑1 (MQ=255)                                 ccGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGc                                                                                                             <  2:330036/150‑1 (MQ=255)                                 ccGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGc                                                                                                             >  2:713633/1‑150 (MQ=255)                                      tGCCGTCCAGCTTAACGAAGCGTGCGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGGTGACCTGCTGATTGCCGATATCAATGGTGTAGTGCGTGGCAAGCGCATAAAGCACACCAGTCTCCACAAGGTTTACGAGAAATGCATCaa                                                                                                        <  2:268560/150‑1 (MQ=255)                                                   aaCGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGcctcc                                                                                           <  1:713633/150‑1 (MQ=255)                                                               gcgTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTg                                                                               <  2:299214/150‑1 (MQ=255)                                                                     cTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAGAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATc                                                                         >  2:386332/1‑150 (MQ=255)                                                                           gAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGt                                                                   >  1:689574/1‑150 (MQ=255)                                                                             aCATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTc                                                                 <  1:206330/150‑1 (MQ=255)                                                                                 ccTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGcctccct                                                                                         <  1:640762/122‑1 (MQ=255)                                                                                 ccTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGcctccct                                                                                         >  2:640762/1‑122 (MQ=255)                                                                                            ctACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCtt                                                                                      <  1:346764/114‑1 (MQ=255)                                                                                            ctACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCtt                                                                                      >  2:346764/1‑114 (MQ=255)                                                                                            ctACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGaatggca                                                  <  1:681208/150‑6 (MQ=255)                                                                                                gtggACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCggg                                              >  1:865269/4‑150 (MQ=255)                                                                                                            aTTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGctggacctggac                                  >  2:647540/1‑150 (MQ=255)                                                                                                            aTTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGAc                                  >  2:647677/1‑150 (MQ=255)                                                                                                             ttGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGc                                             >  1:608973/1‑138 (MQ=255)                                                                                                             ttGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGc                                             <  2:608973/138‑1 (MQ=255)                                                                                                                cAGATATGAATGGGGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCg                              >  2:364840/1‑150 (MQ=255)                                                                                                                  gATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAg                                                                                                               <  2:266436/67‑1 (MQ=255)                                                                                                                  gATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAg                                                                                                               >  1:266436/1‑67 (MQ=255)                                                                                                                  gATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGc                            <  2:865269/150‑1 (MQ=255)                                                                                                                     atGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCCTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGAt                         <  2:589824/150‑1 (MQ=255)                                                                                                                                  tgCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCt            >  1:257966/1‑150 (MQ=255)                                                                                                                                      tGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTa        >  1:420348/1‑150 (MQ=255)                                                                                                                                           aGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTATCCaa   <  2:257966/150‑1 (MQ=255)                                                                                                                                            gcgcATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTATCCAAt  >  2:157322/1‑150 (MQ=255)                                                                                                                                              gcATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGcc                                                                                  <  2:385166/68‑1 (MQ=255)                                                                                                                                              gcATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGcc                                                                                  >  1:385166/1‑68 (MQ=255)                                                                                                                                                   |                                                                                                                                               AACAACAAAGGTGTTGCATCATGTCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTATCCAAT  >  NC_002947/5912768‑5913056

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