Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,095,379	C→T	R20H (CGC→CAC)	insE4 ←	IS3 element protein InsE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,095,379	0	C	T	96.8%	94.5 / ‑3.7	31	R20H (CGC→CAC)	insE4	IS3 element protein InsE
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (14/16);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGG  >  NC_000913/1095257‑1095622                                                                                                                           |                                                                                                                                                                                                                                                    cAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGcc                                                                                                                     >  2:513025/1‑251 (MQ=14)     agaCGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGc                                                                                                                 <  2:482536/251‑1 (MQ=18)             gctgATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTGACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgc                                                                                                           <  2:625635/249‑1 (MQ=255)                gATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTg                                                                                                      >  1:418668/1‑251 (MQ=38)                 aTTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGATTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTg                                                                                                      >  1:269824/1‑250 (MQ=255)                       cTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTAAGCAAGCTTCAGGGCTTAACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGGTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGaa                                                                                                >  2:405906/1‑250 (MQ=255)                        tGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGa                                                                                                 >  1:158298/1‑248 (MQ=255)                         gTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGa                                                                                              <  1:445939/250‑1 (MQ=255)                            tACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACgg                                                                                           <  2:489655/250‑1 (MQ=255)                              cTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGt                                                                                        >  1:590599/1‑251 (MQ=255)                                gcgcCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACGGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGt                                                                                        <  2:549884/249‑1 (MQ=255)                                   ccAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAgg                                                                                   >  2:491882/1‑251 (MQ=255)                                      gTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCa                                                                                 >  2:127880/1‑250 (MQ=255)                                        tGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCTCCTCTGACTTAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACTCGTTTTTGATGGAAGACGGCGTTCAGGCAAt                                                                               <  1:579311/250‑1 (MQ=255)                                        tGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAAt                                                                               <  2:635253/250‑1 (MQ=255)                                         gTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATgc                                                                             <  1:315224/251‑1 (MQ=255)                                          tAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACAACGATGCGTTCAGCAAGCTTCAGGGCTTAACTGTGAAATTCAGGCGAATGATGTTTACGGGGGTTTTTACTGGTTGATAATGTTTTTGACATGTGAGTCACATCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATgcg                                                                            >  2:120474/1‑251 (MQ=17)                                             agTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGACGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAAGCGTTTGATGGAAGACGGCGTTCAGGCAATGcgcg                                                                          >  1:301662/1‑250 (MQ=255)                                                tgtgATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAtt                                                                       <  2:456959/250‑1 (MQ=255)                                                tgtgATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAtt                                                                       <  2:304067/250‑1 (MQ=255)                                                   gATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                                    <  1:484224/250‑1 (MQ=255)                                                   gATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACGGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                                    <  1:260830/250‑1 (MQ=255)                                                      tCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTg                                                                 <  2:321113/250‑1 (MQ=255)                                                       cATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGc                                                                <  2:410565/250‑1 (MQ=255)                                                             ggCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTg                                                         <  1:127880/251‑1 (MQ=255)                                                                    ttCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCCCTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAgg                                                                                   >  2:226099/1‑218 (MQ=255)                                                                    ctCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAgg                                                                                   <  1:226099/217‑1 (MQ=255)                                                                              gCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGcc                                           >  2:252100/1‑248 (MQ=255)                                                                                      aaCACCGATGCGTTCAGCAAGCTTCAGGGCTTCCCTGTGAAATTCAGGCGCATGCTGTTTACGGGGTTTTTTACTGGTTAATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCCCTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGACGACGGCGTTCAGTCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAgg                                >  2:50490/1‑251 (MQ=21)                                                                                      aaCACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAgg                                <  1:550050/251‑1 (MQ=255)                                                                                                                    ttCACTGTGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTAATCACTTAGCCGCGTGTCCACTCTTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGACATCGCTTCGCCGGAGCATAAGGTTCTGATGAACGTAACAGCAAAAAGCGAgg  >  2:364085/1‑251 (MQ=255)                                                                                                                           |                                                                                                                                                                                                                                                    CAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGG  >  NC_000913/1095257‑1095622

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