Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	16,749	T→A	intergenic (+192/+2)	insL1 → / ← hokC	IS186/IS421 family transposase/protein HokC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	16,749	0	T	A	96.9%	91.5 / ‑6.3	32	intergenic (+192/+2)	insL1/hokC	IS186/IS421 family transposase/protein HokC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (25/6);  minor base C (0/1);  total (25/7)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGG‑TCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGT  >  NC_000913/16502‑16983                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           tCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGG‑TCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACtt                                                                                                                                                                                                                                         >  1:348247/1‑251 (MQ=25) cccagccagCGCTGGATTTACCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGG‑TCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑ACGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACtt                                                                                                                                                                                                                                         >  1:348245/2‑251 (MQ=255)             tGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAACTGG‑TCATCTGGAGCTTACCGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTCCTCGGATTCGt                                                                                                                                                                                                                             >  2:274362/1‑251 (MQ=255)              ggATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGG‑TCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGt                                                                                                                                                                                                                             >  2:456250/1‑250 (MQ=255)               gATTTCCCCCCAAGATGGCCCGGATCGAAAAAGAAGACCTACCT‑GGTGGGGAAAAAAAAAAAAATGGGTCATTTGGAGTTTAAAGGTTGCCA‑TCCGTGGGACAGTATCCCTGACAGCCTACAAAGCGCAATTGAAGAACGCGAGGCTCCGTCTTAAGGAGCCCCCGCGGCGGCGCATTCTACAGATGGTT‑AACCCTTGAGGTTAGCGTTTATGGGATTACCCCCCGCCGTTGCACTTACTCGGATTCGTaa                                                                                                                                                                                                                           <  2:570309/251‑1 (MQ=9)                           aGAAGTGCCGGAGCCGAAAAGAAGAACTAACGCGTTGTGGAGAATAACAAAAATGG‑TCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGACTCCCTGACAGCCTCCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTATCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAAACCTT‑AAGCTCGCGCTTATGGGATCACTCCCAGCCGTTGCACTTACTCGGCTTCGTAAGCCGTGaaa                                                                                                                                                                                                                  >  1:575624/1‑248 (MQ=14)                                    gtaaccGAAAAAAAAAACTAACTCGTTGGGGAAAATAACAAAAAGGG‑TCACCGGGAGCTGACGGGTGGCCA‑TGCGTGGGACAGTAGCCCGGACAGCCTCCAAAACGCAATGGAAGAACGGGAGGAAGCGCCTAAACGAGGCGCCGAGGCGTCGCTTTCTCCAGGTGTTTCACCCCTT‑AAGTGAGGCCTTATGGGACCCCGCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTcc                                                                                                                                                                                                      <  1:274362/247‑1 (MQ=255)                                          gAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGG‑TCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTgg                                                                                                                                                                                                >  2:145000/1‑251 (MQ=255)                                            aaaGAAGAACAAACTCGTGGGGGAGAATAACAAAAATGG‑CGATCTGGCGTTTACAGG‑GGCCAATTCGTGGAACAGTACCCCTGACAGCCTACAAACCGCAATTGAAGAACGCGAGGCAGCGTCATAACGAGGCACCGAGGCGTCGCATGCTTCAGATGGTTCACCCCTA‑AAGTTAGCGCTTAGGTGATCACTCCCCGCCGTTGCACTTACTCGGATGCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGcc                                                                                                                                                                                              <  2:259940/251‑1 (MQ=255)                                                          cctttggagaaaaTAACAAAAATGG‑TCATCTGGGCCTTACAGGTGGCCA‑TTCGGGGGACAGTTCCCCTGACAGCCAAAAAAACGAAATTGAAGAACGCGAGGCACCGTCTTACCGGGCCACCGGGGCGTCGCATTCTTCAGATGGTCCACCCTTT‑AAGTTAGCGTTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGaa                                                                                                                                                                                <  1:111476/240‑1 (MQ=255)                                                                          cAAAAATGG‑TCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCttt                                                                                                                                                                >  1:48922/1‑251 (MQ=255)                                                                                  g‑TCATCTGGAGCTTACAGGCGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGCGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAACCAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGtt                                                                                                                                                         >  1:325097/1‑250 (MQ=255)                                                                                            aGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGcgccg                                                                                                                                               >  1:431476/1‑251 (MQ=255)                                                                                                   aGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACACAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTGACCAGCGCCGCCACTAc                                                                                                                                        >  2:223300/1‑251 (MQ=255)                                                                                                                  gggACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATAc                                                                                                                             >  2:269246/1‑248 (MQ=255)                                                                                                                      cAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGAttac                                                                                                                      >  1:214361/1‑248 (MQ=255)                                                                                                                        gTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAACACAGCAACCTCCGTCTGGCCAGTTCGGCTGTGAACCTCACAGAGGTCttt                                                                                                                                                                >  1:305047/1‑207 (MQ=255)                                                                                                                              cTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTCGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACggtccggg                                                                                                              >  2:263560/1‑244 (MQ=255)                                                                                                                                            aaCGCAATTGACGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGcc                                                                                                >  1:480018/1‑251 (MQ=255)                                                                                                                                                 aaTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATAAAGATGACGATccagaacgggaaagcacccgtctg                                                                                           >  1:165702/1‑228 (MQ=37)                                                                                                                                                   ttGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGGGATACAGATGACGATCAGGGCGACAATCATCGCCTTAtgct                                                                                         >  1:532857/1‑251 (MQ=255)                                                                                                                                                                     cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGCCGATCCGGGCGACAATCATCGACTTATGCTGCTTCATTGCTCtcttct                                                                       >  1:445138/1‑251 (MQ=255)                                                                                                                                                                           aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGGTACAGATGACGATCAGGGAGACAATAATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTg                                                                  >  1:342178/1‑250 (MQ=255)                                                                                                                                                                                      gAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGCACTTACgggcc                                                      >  2:470877/1‑248 (MQ=255)                                                                                                                                                                                             cGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCCGGGCGCCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAag                                                 >  1:303084/1‑249 (MQ=255)                                                                                                                                                                                              gCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGCATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTAGGATGTGAACCTCACAGCGGTCTTTTCTCGTTACAAGCGCCGCCACTACGGCGGTGATACCGATGACGATCAGGGGGCAAAACATAGCCTTATACTGATTCATTGCTCTCTTTTCCTTGCACTTACGGtacgaaaaagg                                              >  1:564695/1‑241 (MQ=255)                                                                                                                                                                                                             ttCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGGGATACAGATGACGATCAGGGCGACAATCATCGCCTTATTCTTCCTCATTTCTCTCTTCTCCTTGACCTTTAGGTCCGTAAAAGGCACTCTAcctttgtt                               >  1:55128/1‑244 (MQ=255)                                                                                                                                                                                                                            tactttttgtggtccccccccccGCCGTTCCACTTACTCGGATTCGGAAGCCGGAAAAACAGCATCCTCCGTCGGGCCAGTGCGGATGTGACCCTCACAGAGGTCTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATagggg                  <  2:348247/231‑3 (MQ=255)                                                                                                                                                                                                                            tAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGGTAAAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACAt                                    >  1:100599/1‑232 (MQ=255)                                                                                                                                                                                                                                    agggaacccccccccGCCGTTGCCCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGt         <  1:171849/240‑1 (MQ=255)                                                                                                                                                                                                                                     gggtatccccccccGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGtt        <  1:456250/241‑1 (MQ=255)                                                                                                                                                                                                                                            aCTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTggtagt  >  1:15036/1‑245 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           TCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGG‑TCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTT‑AAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGT  >  NC_000913/16502‑16983

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