Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,507,634	G→A	N328N (AAC→AAT)	insI‑3 ←	KpLE2 phage‑like element; IS30 transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,507,634	0	G	A	95.5%	69.8 / ‑3.4	22	N328N (AAC→AAT)	insI‑3	KpLE2 phage‑like element; IS30 transposase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (8/13);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGATCGCAGAAGTAAACTTTAACGCC  >  NC_000913/4507400‑4507693                                                                                                                                                                                                                                           |                                                            tAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATtct                                                         <  2:1138322/239‑1 (MQ=255)       gAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGt                                                  >  2:393538/1‑240 (MQ=255)        aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGtt                                                 >  2:626292/1‑240 (MQ=255)           aaaaaGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTccc                                              <  2:927555/240‑1 (MQ=255)           aaaaaGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTccc                                              <  1:1087164/240‑1 (MQ=255)           aaaaaGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTccc                                              <  1:1155629/240‑1 (MQ=255)             aaaGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCg                                            <  2:798140/240‑1 (MQ=255)                     aaCGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCaa                                     <  1:266849/239‑1 (MQ=255)                     aaCGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAg                                    <  1:463257/240‑1 (MQ=255)                       cGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAGGa                                  >  2:683717/1‑240 (MQ=255)                          aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCAAGGAct                                >  2:509028/1‑239 (MQ=255)                          aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAGGAct                                >  1:1015846/1‑239 (MQ=255)                            tGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAGGActct                              <  2:675870/239‑1 (MQ=255)                              aaaTGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAg                           <  1:929759/240‑1 (MQ=255)                                    agaATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGATCGc                     <  1:529047/240‑1 (MQ=255)                                    agaATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCAGAATTAGCCTATTTGTATTCTTATATg                                                   >  2:80775/1‑210 (MQ=255)                                    agaATAAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCAt                                                      <  2:670011/207‑1 (MQ=255)                                      aaTGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCg                                            <  1:769004/215‑1 (MQ=255)                                                tAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGATCGCAGAAGTAAACtt         >  2:81528/1‑240 (MQ=18)                                                   aaaGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGATCGCAGAAGTAAACTTTaa      <  1:683720/240‑1 (MQ=17)                                                       gCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGATCGCAGAAGTAAACTTTaa      >  1:1091843/1‑236 (MQ=11)                                                       gCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTATTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGATCGCAGAAGTAAACTTTAACGcc  >  1:579590/1‑240 (MQ=11)                                                                                                                                                                                                                                           |                                                            TAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGATCGCAGAAGTAAACTTTAACGCC  >  NC_000913/4507400‑4507693

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