Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,535,060	G→A	intergenic (‑1007/‑1537)	topAI ← / → ythB	toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,535,060	0	G	A	92.3%	28.2 / ‑6.1	13	intergenic (‑1007/‑1537)	topAI/ythB	toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (9/3);  total (9/4)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGG  >  NC_000913/4534879‑4535235                                                                                                                                                                                      |                                                                                                                                                                                ccAAGTTCNANNNCNNAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGNCTACTGTAACTCACTNTGNTATANCGGNAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAAGGAGACGTCATANCCCNCTTGANNCNGAAACANNngcnnnnnnnngnnnnggnnann                                                                                                                        >  2:179918/1‑219 (MQ=255)             ncgGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAAGGAGc                                                                                                                                                                            >  1:76527/2‑175 (MQ=255)               ggAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGCTATTTACATATTGATGGTAGAGGCATGAAGCAAAATGGAAGGAGcc                                                                                                                                                                           >  1:183044/1‑174 (MQ=255)                    gTATCTGTATGAGCATCGTCGCTGTTACAATAACATCATCGACTACGGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACAGATTGATGGTAGAGGCATGAAGCCAAACGGAAGGAGCCGTCATAACCCGCTTGAAGCGGAATCAATAGCAGGATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAAcc                                                                                         <  1:179918/251‑1 (MQ=255)                                         cTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAACGAGGGATGGAAAAAGGACGATTTTCCCCGCACTGGGATATTTACATATTGATGGTAGAGGCATGAAGCCCAATGGAAGGCGCCGTCATAACCCTCTTGAAGCTGAAACAATCGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGgtgt                                                                    >  1:9370/1‑251 (MQ=255)                                                                                  aTCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAACGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCCTGAAGCCAAATGGAAGGAGCCGTCATAACACTCATGAAGCTGAAACAATAGCAGCATGGATGGTGGc                                                                                                                             >  1:220293/1‑153 (MQ=255)                                                                                            gTTGCAGCCTAAAAGAGGGAGGGAAAAAGGACGATTGTCCCCGCAATGGGATATTGACATATGGATGGTAGAGGCATGAAGCCAAATGGAAGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCtggtg                                                                                                                               <  1:231327/141‑1 (MQ=255)                                                                                                  gCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAAGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTgg            >  1:180707/1‑250 (MQ=255)                                                                                                                 ggAAAAAGNACGATTTTCCCCGAGATGGGAGATTTACATANTNACGGGCGAGGAATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGNAGCTGAAGCANTAGNNGCAGNGNTNGTGGCACATAAGGATGATAGAGGGCGCCANTGCTGTGAACCCCTGGACAAAGTNNGTNNNGTNNTGACg                                                             <  2:193152/186‑1 (MQ=255)                                                                                                                                                                    ngCATGAAGCCAAATGGAAGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATgg  >  1:32314/2‑194 (MQ=255)                                                                                                                                                                      cATGAAGCCAANTGGAAGGAGCCGTCATAAccc                                                                                                                                                                 >  2:36610/1‑33 (MQ=255)                                                                                                                                                                      cATGAAGCCAAATGGAAGGAGCCGTCATAAccc                                                                                                                                                                 <  1:36610/33‑1 (MQ=255)                                                                                                                                                                        tGAAGCCACATGGAAGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCgg                                                   >  1:213650/1‑141 (MQ=255)                                                                                                                                                                                      |                                                                                                                                                                                CCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGG  >  NC_000913/4534879‑4535235

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