Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,077,384	G→A	T91T (ACG→ACA)	cbeA →	cytoskeleton bundling‑enhancing antitoxin CbeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,077,384	0	G	A	97.1%	93.9 / ‑5.0	34	T91T (ACG→ACA)	cbeA	cytoskeleton bundling‑enhancing antitoxin CbeA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/23);  minor base C (0/1);  total (10/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACGCTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCTACGGCCTGGTGTTTTCACCACGCCACTTTTCCATTTTTATATCTGCATATCAGGAAAATCTTCAGTATGAAAACATTACCTGTATTACCCGG  >  NC_000913/2077141‑2077594                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                   tCCCGCCGACAATCACGCCCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTa                                                                                                                                                                                                               <  1:531737/249‑1 (MQ=255)     gacgacACTCACGCCCGCCCCTGGTGGGGGCTGCCCTGCACCTTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAAGCCCCGCCATCAGCATACCGTCACACTGTATGCaaa                                                                                                                                                                                                         <  2:229526/251‑1 (MQ=255)     gacgacAATCACGACCGCCCCTGGGGGGGGCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAATCAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCaaa                                                                                                                                                                                                         <  2:483634/251‑1 (MQ=255)                 gACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTg                                                                                                                                                                                             >  1:457843/1‑251 (MQ=255)                      ccccTGGGGGGGGCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAg                                                                                                                                                                                        <  1:380424/251‑1 (MQ=255)                       cccTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGc                                                                                                                                                                                       <  1:541684/251‑1 (MQ=255)                           ggtggGGGCTGCCCTGCACCGTGCCGCCCTGTGTCGGGGCACGTATGGTGCAGGAGGTTAACCGGTGTCAGTCCCGTGCCGCCCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGTCCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAAGCCCCGCCATCAGCATACCGTCACCCTGTATGCAAAAGGGCTGACCTGCAAAGCCg                                                                                                                                                                                     <  1:305084/249‑1 (MQ=255)                                 ggCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTc                                                                                                                                                                             >  1:260085/1‑251 (MQ=255)                                        cTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGtt                                                                                                                                                                      >  1:115917/1‑251 (MQ=255)                                           cACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTtgtg                                                                                                                                                                   <  1:119962/251‑1 (MQ=255)                                            aCCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTtgtg                                                                                                                                                                   >  1:466644/1‑250 (MQ=255)                                               gtgccGCCCTGTTGCGGGGCAGGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGGCCGCGCCGGTATCAGAGGCCTGTTCAGCGATGAAGATGCGTACCACCTAGGCCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCTTTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGGATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTa                                                                                                                                                               <  2:462087/247‑1 (MQ=255)                                                         gTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTCCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCt                                                                                                                                                      <  2:598668/250‑1 (MQ=255)                                                             cGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCTCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAACCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTGTATCTGGCTg                                                                                                                                                 <  2:260085/251‑1 (MQ=255)                                                                                      aaCCGGTTGCATTACCGTGCCGACCGCGCCGGTAGCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCGGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCTCCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGaaaaa                                                                                                                        <  1:104505/251‑1 (MQ=255)                                                                                      aaCCGGTTGCAGTCCCTTGCCGCCCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCCCCTGGACAAGGCCTTTCGGCTGCGGATTAAACAACTGGAACCCCTGCTCACCAGCGGTGAACTGAATCCCCGACATCAGCATACCGTCTCACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATCCCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTCTCCGACGCCCGATATGaaaa                                                                                                                         <  2:534317/250‑1 (MQ=255)                                                                                          ggtagctttcccTTGCCGCCCGCGCCGGGATCAGGGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGAAGAAACAACTGAAACTCATGCTCGCCATCGGTTAACTGAATCCCCGCCATCAGCATACCGTCACACTGGATGCAAAAGGTCTGACCTGCAAAGCCGATACCCGCAGCAGTTGTGATGACGATTATCGGGCTGTTTATCCGACGCCCGAAAGGAAAAATTaa                                                                                                                    <  2:296490/241‑1 (MQ=255)                                                                                                   ccctggCCGACCGCGCCGGTATCAGAGGCGTGTTCAGCGATGCAGAGGCGTTCCACCTGGACTAGGCCGTTCCGCTGCTGAGGAAACAACTGGATCTCATGCTCACCAGCGGTGAACTGAAGCCCCGCCAGCAGCATACCGTCACACTGTAGGCAAAGGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGGTTATCCGCCGCCCGACATGAAAAATTAACTGTCCAg                                                                                                            <  2:629630/245‑1 (MQ=255)                                                                                                         ccGACCGCGCCGGTATCAGAGGCCTGTTCATCTATGCAGATGCGTACCCCCTGGACCAGGCCTTCCCGCTGCTGATGAAACAACTGGACCTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACGGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGcc                                                                                                     <  2:494270/251‑1 (MQ=255)                                                                                                                cgcaggTAGCAGAGGCCTGGTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCACGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCt                                                                                              <  2:49313/247‑1 (MQ=255)                                                                                                                    ggTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCTACgg                                                                                          <  1:594912/251‑1 (MQ=255)                                                                                                                    ggTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCGTCTGCTACgg                                                                                          <  2:231236/251‑1 (MQ=255)                                                                                                                     gtaccAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTAATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGGAGGCAAAAGTGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGAGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGaaaaa                                                                                                                        <  2:589313/216‑1 (MQ=255)                                                                                                                     gTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGaaaaa                                                                                                                        >  1:589313/1‑220 (MQ=255)                                                                                                                                             cAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCTACGGCCTGGTGTTTTCACCACGCCACttt                                                                 >  2:570884/1‑251 (MQ=255)                                                                                                                                                       ccaccTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGcc                                                                                                                                                                                      >  1:206723/1‑124 (MQ=255)                                                                                                                                                       ccaccTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACGGTATGCAAAAGGGCTGACCTGCAAAGcc                                                                                                                                                                                      <  2:206723/124‑1 (MQ=255)                                                                                                                                                          ccTGGACCCGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCTACGGCCTGGTGTTTTCACCACg                                                                        >  1:165927/1‑231 (MQ=255)                                                                                                                                                          ccTGGACCAGGCCTGTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAACCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCTACGGCCTGGTGTTTTCACCACg                                                                        <  2:165927/231‑1 (MQ=255)                                                                                                                                                             ggACCAGGCCTTTCCGCTGCTGATGAAACACCTGGAACCCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCTACGGCCTGGTGTTTGCACCACGCCACTTTTCCATTTTTATATc                                                   <  2:131331/249‑1 (MQ=255)                                                                                                                                                                       tttCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGAACTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGATGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATCGACTTCTGCTACGGCCTGGTGTTTTCACCACGCCACTTTTCCATTTTtatctcatcatttcagg                                        >  2:405722/1‑236 (MQ=255)                                                                                                                                                                               gctaatgaAACACCTGGAACTCATGCTCACCAGCGGGGAACTGACGCCCCGCCATCAGCATACCGTCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGGTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCTACGGCCTGGTGTTTTCACCACGCCACTTTTCCATTTTTATATCGGCATATCAGGAAAATCt                                 <  2:216082/245‑1 (MQ=255)                                                                                                                                                                                                accTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGGCACACTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCATTTGTGAGTACGTTTAGCTGGCTGTTTATCGGACGCCCGAAATGAAAAAGTAACTCTCCAGAATAGCCTTCTGCTACGGCCTGGTGTTTTCACCACGCCACTTTTCCATTTTTATATCTGCGTATCAGGCAAGTCTGCAGTATGAAAACAtttcc              <  2:120692/249‑4 (MQ=255)                                                                                                                                                                                                            ccAGCGGTGAACTGAATCCCCGCCATCAGCATACCGACACACTGTATGCAAACGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCTACGGCCTGGTGTTTTAACCACGCCACTTTTCCATTTTTATATCTGCATATCAGGAAAATCTTCAGTATGAAAACATTACCTGTATTACCCgg  >  2:602704/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                   TCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGGCCTTTCCGCTGCTGATGAAACAACTGGAACTCATGCTCACCAGCGGTGAACTGAATCCCCGCCATCAGCATACCGTCACGCTGTATGCAAAAGGGCTGACCTGCAAAGCCGATACCCTCAGCAGTTGTGATTACGTTTATCTGGCTGTTTATCCGACGCCCGAAATGAAAAATTAACTCTCCAGAATAGCCTTCTGCTACGGCCTGGTGTTTTCACCACGCCACTTTTCCATTTTTATATCTGCATATCAGGAAAATCTTCAGTATGAAAACATTACCTGTATTACCCGG  >  NC_000913/2077141‑2077594

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