Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	657,294	G→A	M1M (ATG→ATA) †	cspE →	transcription antiterminator and regulator of RNA stability CspE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	657,294	0	G	A	96.4%	91.1 / ‑5.4	28	M1M (ATG→ATA) †	cspE	transcription antiterminator and regulator of RNA stability CspE
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (15/12);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATC  >  NC_000913/657071‑657535                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  acaacaATGGCAATGTGTTCTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTg                                                                                                                                                                                                                         <  1:37056/250‑1 (MQ=255)       aTGGCAATTTGTACGTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAATTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTAGTTTGCATGTAAAGGTAATTTTGATATCTAAGAt                                                                                                                                                                                                                                           <  2:50511/226‑1 (MQ=255)       aTGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGAt                                                                                                                                                                                                                                           >  1:50511/1‑226 (MQ=255)             aTGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTcc                                                                                                                                                                                                            >  1:275299/1‑251 (MQ=255)                         ccTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGt                                                                                                                                                                                                >  1:69520/1‑251 (MQ=255)                                gcgcTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCAtt                                                                                                                                                                                          <  2:69520/250‑1 (MQ=255)                                           tttGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTATATTTTTCATGTAAAGGAAATTTTGATATCTAAGATTAAAGGTAACGTTATGTGGTTTATTGAGTCCAAAGGAATCGGTTTCATTACTCCGGCAGAc                                                                                                                                                                              >  1:311547/1‑251 (MQ=255)                                              gagaCAAATGAAGTTTTATAAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTATCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACgg                                                                                                                                                                            <  1:1422/250‑1 (MQ=255)                                                      tGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCaa                                                                                                                                                                      >  1:307878/1‑248 (MQ=255)                                                        aaGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACg                                                                                                                                                                 >  2:4625/1‑251 (MQ=255)                                                             ttAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTc                                                                                                                                                            >  1:71039/1‑251 (MQ=255)                                                                atAACATCTTTAAAATCAATAACTAAAATCAGTAACAGCAAAAATAACGCGACTTTTATCACTTTGTAGTAACGTTACACTGGACAAAGCGTACCAATATTTGTGTACTGGTATCCGATGCAGCATTGGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAATGATGCGGTTTCATTACTCCTGAAGACGGCAGCACACTCGTGTTCGTa                                                                                                                                                         <  2:49064/250‑1 (MQ=255)                                                                      tcttTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtc                                                                                                                                                   >  2:283285/1‑251 (MQ=255)                                                                       cttTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct                                                                                                                                                  >  1:117114/1‑251 (MQ=255)                                                                                          aaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAACGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATgg                                                                                                                               >  2:188204/1‑251 (MQ=255)                                                                                                         aaaTACCGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCGCAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGc                                                                                                                <  2:319225/251‑1 (MQ=255)                                                                                                           aTAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                                                                                                              >  1:8213/1‑251 (MQ=255)                                                                                                                        tATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCATAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTACCGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGt                                                                                                 >  2:162332/1‑251 (MQ=255)                                                                                                                            aCTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTagag                                                                                             <  2:71039/251‑1 (MQ=255)                                                                                                                                      aaaGTTACTCTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATc                                                                                    <  1:177328/250‑1 (MQ=255)                                                                                                                                        aGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCCGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACt                                                                                 <  1:283285/251‑1 (MQ=255)                                                                                                                                         gTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACt                                                                                 >  1:253732/1‑250 (MQ=255)                                                                                                                                                               acaatagagtgaCTGGTAACCGACACAGCATTTGAGTCCATTTTTCATGTAAAGGTAAGTTTGATGTACAAGATTAAAGTTAACGTTAAGTGGTTTAATGAGTCCTAAGGATTAGGTTTCATTACTCCGGAAGACGGCAGCAAAGAGGTGTTCGTACACTACTCTGAACTCCAGCCTAATGGTTTTAAAACTCTTGCTGAAGTTCTGCGCGAAGAGTTCGAAATAACAAACGGTCCTACAGGCACTTctg                                                           >  2:55765/12‑250 (MQ=255)                                                                                                                                                                caATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTctgct                                                         <  2:117114/251‑1 (MQ=255)                                                                                                                                                                                    gACACAGCATTTGTGTCTATTTTTCATGGAATGGTAATTTTGATATATAAGAATAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACTGCAGCAAATACGTGTTCGTACACTTCGCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTCCGAAAGCACTAACGGTGCCAAAGGCCCTTCTTCTGCAATCGTAATCGCTCTGtt                                     <  2:215476/251‑2 (MQ=255)                                                                                                                                                                                                            cATGTAAAGGTAATTTTGCTATCTAAGATTAAAGGTAACGTTACGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCaaa                                                                     >  2:204397/1‑195 (MQ=255)                                                                                                                                                                                                            cATGTAAAGGTAATTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCaaa                                                                     <  1:204397/195‑1 (MQ=255)                                                                                                                                                                                                                        aTTTTGATATCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGATGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAACCCCGCTTAATc  >  2:177135/1‑250 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  ACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATC  >  NC_000913/657071‑657535

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