Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,534,255	G→T	intergenic (‑202/+175)	topAI ← / ← yjhZ	KpLE2 phage‑like element; toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/KpLE2 phage‑like element; putative acetyltransferase YjhZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,534,255	0	G	T	96.2%	79.6 / ‑2.8	26	intergenic (‑202/+175)	topAI/yjhZ	KpLE2 phage‑like element; toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/KpLE2 phage‑like element; putative acetyltransferase YjhZ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (12/13);  total (13/13)
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.26e-01

CAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGGGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATGGCCTCATCGGCTTTTCCCTTCGGG  >  NC_000913/4534136‑4534389                                                                                                                        |                                                                                                                                       cAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGaaaaa                                                                                                                    <  2:18462/140‑1 (MQ=255)    tCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAgtg                                                                                                                 >  2:356717/1‑140 (MQ=255)    tCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAgtg                                                                                                                 >  1:373528/1‑140 (MQ=255)                                              aTTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGaaa                                                                       <  1:412386/140‑1 (MQ=255)                                                tctcCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAg                                                                     >  1:313108/1‑140 (MQ=255)                                                    cAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAg                                                                 >  2:277688/1‑140 (MQ=255)                                                    cAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAg                                                                 >  2:17857/1‑140 (MQ=255)                                                     aGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGa                                                                <  1:17857/140‑1 (MQ=255)                                                            aaGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGGGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACt                                                         >  1:299179/1‑140 (MQ=255)                                                                   gggTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAAt                                                  <  2:89555/140‑1 (MQ=255)                                                                     gttgagttgaAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCGGATGACACTGACGAATGc                                                <  2:373528/140‑1 (MQ=255)                                                                         agttgaAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAGGGCAGATGACACTGACGAATGCATGt                                            <  1:356717/140‑1 (MQ=255)                                                                              aaaGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAg                                       <  1:446763/140‑1 (MQ=255)                                                                              aaaGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAg                                       <  1:280189/140‑1 (MQ=255)                                                                                   tATTTACGGAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGc                                  <  2:171870/140‑1 (MQ=255)                                                                                          ggAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGt                                      >  2:179106/1‑129 (MQ=255)                                                                                          ggAGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGt                                      <  1:179106/129‑1 (MQ=255)                                                                                            aGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATgg                         >  1:157772/1‑140 (MQ=255)                                                                                            aGCGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATgg                         >  2:266406/1‑140 (MQ=255)                                                                                              cGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGc                                  <  2:18948/129‑1 (MQ=255)                                                                                              cGAATATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGc                                  >  1:18948/1‑129 (MQ=255)                                                                                                aaTATTAACACGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATGGCCTc                     <  2:273240/140‑1 (MQ=255)                                                                                                          cGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATgg                         <  2:21984/126‑1 (MQ=255)                                                                                                          cGTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATgg                         >  1:21984/1‑126 (MQ=255)                                                                                                           gTAGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATGGCCTCATCGGCTTTTc          >  2:275128/1‑140 (MQ=255)                                                                                                             aGTTATAGGAGTGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATGGCCTCATCGGCTTTTccc        >  1:328515/1‑140 (MQ=255)                                                                                                                   aggagtGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATGGCCTCATCGGCTTTTCCCTTCggg  >  2:236081/1‑140 (MQ=255)                                                                                                                        |                                                                                                                                       CAGTCCGGCCTGAAAAGGTCGGGGAACTGATTATCGGATGATTACATTCTCCAGCATCAAAGCCTCGGGTTGAGTTGAAAGGTATTTACGGAGCGAATATTAACACGTAGTTATAGGAGGGATATGTTTCAGTAGAAAAAGTGTGATGTACTTTGGAAATTGAATGCCAACCCATTTACTGGAAAAGGCAGATGACACTGACGAATGCATGTTGAAGTACGCCTGAAATGGCCTCATCGGCTTTTCCCTTCGGG  >  NC_000913/4534136‑4534389

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