Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,503,822	A→G	14.8%	intergenic (‑391/‑236)	insG ← / → yjhB	IS4 transposase/putative MFS transporter, membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,503,822	0	A	G	14.8%	63.5 / 5.6	27	intergenic (‑391/‑236)	insG/yjhB	IS4 transposase/putative MFS transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base A (10/13);  new base G (2/2);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

GCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATA  >  NC_000913/4503692‑4503957                                                                                                                                   |                                                                                                                                        gCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGttt                                                                                                                      >  1:217952/1‑150 (MQ=255)    aaaGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGtttttt                                                                                                                   <  2:115676/150‑1 (MQ=255)    aaaGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGtttttt                                                                                                                   <  1:37816/150‑1 (MQ=255)                tctcAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCa                                                                                                       >  1:20233/1‑150 (MQ=255)                tctcAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCa                                                                                                       >  1:40354/1‑150 (MQ=255)                    aGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGTCCCaaatc                                                                                                   >  2:248752/1‑148 (MQ=255)                    aGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAAc                                                                                                   >  2:95794/1‑150 (MQ=255)                    aGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAAc                                                                                                   >  2:95791/1‑150 (MQ=255)                    aGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAAc                                                                                                   >  2:7142/1‑150 (MQ=255)                     gTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAGTATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACt                                                                                                  >  1:109927/1‑150 (MQ=255)                     gTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAGTATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACt                                                                                                  >  1:109929/1‑150 (MQ=255)                                                                     cGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTaagaa                                                  <  2:235332/150‑1 (MQ=255)                                                                         aaTGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAGTATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGaaa                                              <  2:109927/150‑1 (MQ=255)                                                                         aaTGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAGTATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGaaa                                              <  2:109929/150‑1 (MQ=255)                                                                                  aaTGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGaaa                                     >  1:98754/1‑150 (MQ=255)                                                                                  aaTGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGaaa                                     <  2:156348/150‑1 (MQ=255)                                                                                  aaTGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGaaa                                     >  1:190444/1‑150 (MQ=255)                                                                                     gAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGTTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGa                                  <  2:170327/150‑1 (MQ=255)                                                                                     gAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGa                                  <  1:76475/150‑1 (MQ=255)                                                                                             cAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTg                          <  2:7212/150‑1 (MQ=255)                                                                                             cAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTg                          <  2:133598/150‑1 (MQ=255)                                                                                                                   cTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAAcaatcaa    <  1:101097/150‑1 (MQ=255)                                                                                                                   cTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAAcaatcaa    <  1:142402/150‑1 (MQ=255)                                                                                                                    tgtgCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTCTTGATAAcaatcaat   >  1:9526/1‑150 (MQ=255)                                                                                                                     gtgCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATa  <  1:95794/150‑1 (MQ=255)                                                                                                                     gtgCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATa  <  1:95791/150‑1 (MQ=255)                                                                                                                     gtgCCAGTGGTGGAATATCTGGCGAACATTGTATTATGTTGGTGACCCAAAACTGTAGGTGGGTATCAGATTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATa  <  1:248752/150‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        GCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATA  >  NC_000913/4503692‑4503957

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑