Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	316,519	G→A	intergenic (+66/+1207)	insF1 → / ← rclC	IS3 element protein InsF/DUF417 domain‑containing inner membrane protein RclC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	316,519	0	G	A	97.0%	108.2 / NA	33	intergenic (+66/+1207)	insF1/rclC	IS3 element protein InsF/DUF417 domain‑containing inner membrane protein RclC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/16);  minor base C (1/0);  total (17/16)
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.83e-01

AGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTGTCTCAGGAAAGTTGATTAATCTACACTCCGTCTT  >  NC_000913/316280‑316764                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       aGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGGCCATATTACGTGGGGAGGATCATAAAGTCCTCGGGTCTTTGGCC‑ACTCTTTTGGGATATCAGCATAtt                                                                                                                                                                                                                                              >  2:365962/1‑250 (MQ=255)                   aaaGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCaaa                                                                                                                                                                                                                           <  1:128805/251‑1 (MQ=255)                             gTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAAcc                                                                                                                                                                                                                 <  2:411942/251‑1 (MQ=255)                              tATCCAAGGAGAACCATTTATCCGCCGGGAAATAATGCGGGCAACGGTGTTTCATTCTATAGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGACCGGAACAATTTGCAAACAAGACACTCGCTTAGGGCTGTGTCCATACTACGTGGGTAGGAACATAAAACCCTCGGGTCGATGTCC‑ACTCTTTTGGCCTATCAGAATATTCACCTACCGTACCAGCAAACGtttaaaacc                                                                                                                                                                                                                >  2:162353/1‑245 (MQ=14)                                        gaACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATc                                                                                                                                                                                                      >  1:350785/1‑251 (MQ=255)                                                       cGGGAAATAATGCGGGCACCGGTGTTTAATTATATCGACTGTGATTACAATCGGTGGCGGCGGCACGGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTCGGCCTTTTCCCCTATTACGAGGGGAGGATCATAAAGACCTCGGGCCGTTGTCC‑ACTCTTTTGCCATATCTGCCTATTCCCCTACCGTCCCAGCCAACGTTTCAACCGTGAGAAAATCGCCCGATTttgttt                                                                                                                                                                                        >  2:374113/1‑249 (MQ=17)                                                       cGGGAAATAATACGGGCAACGGAGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCATCATTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCATCGGGTCGTTGTCCCAC‑CTTTTGGCATATAAGCATATTCACATACCGGACCCGCAAAAGTTTCAACCGTGAGCAAATCGACCGATAttgttgc                                                                                                                                                                                       >  2:374141/1‑250 (MQ=18)                                                            aataatGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑tt                                                                                                                                                                                 <  2:301894/251‑1 (MQ=255)                                                                    gggCAGCGGCGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGG‑TCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTgg                                                                                                                                                                        >  2:239937/1‑251 (MQ=255)                                                                        aaCGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGaaa                                                                                                                                                                     <  1:300026/251‑1 (MQ=255)                                                                         aCGGTGTTTAATTATATCGAATGGGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCGCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGTCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAAc                                                                                                                                                                    <  1:356901/251‑1 (MQ=255)                                                                                      atatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATcc                                                                                                                                                       <  1:159344/251‑1 (MQ=255)                                                                                      atatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATcc                                                                                                                                                       <  2:217411/251‑1 (MQ=255)                                                                                      atatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATcc                                                                                                                                                       <  2:145316/251‑1 (MQ=255)                                                                                                    tACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTGCAACCGTTAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTTGAAACGAAATATATATCCTTCTTAAAAttct                                                                                                                                          <  2:239322/250‑1 (MQ=255)                                                                                                           ggTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACAt                                                                                                                                  <  2:93941/251‑1 (MQ=255)                                                                                                                   cggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAAACTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTa                                                                                                                           >  1:52100/1‑250 (MQ=255)                                                                                                                                                ggAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCAt                                                                                             >  1:359393/1‑251 (MQ=255)                                                                                                                                                 gAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGAt                                                                                                                                                                                               <  2:193768/153‑1 (MQ=255)                                                                                                                                                 gAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGAt                                                                                                                                                                                               >  1:193768/1‑153 (MQ=255)                                                                                                                                                       tttGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTT‑GTTTACGGGTTCAAACGTGGAA‑CGGAATAGATATCCGTCTTAAAATTTTTCCTTCCTCTTGCTTAATGAATACGGCGTTGTTTTActtctttcttcggtggt                                                                                     >  1:19473/1‑240 (MQ=255)                                                                                                                                                           aaaaCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGAttt                                                                                   >  2:14254/1‑250 (MQ=255)                                                                                                                                                            aaaCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATCTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTAt                                                                                 >  1:190674/1‑251 (MQ=255)                                                                                                                                                               cAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCg                                                                              >  1:30010/1‑251 (MQ=255)                                                                                                                                                                     ccTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACtt                                                                        <  1:85311/251‑1 (MQ=255)                                                                                                                                                                      cTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTc                                                                       >  2:88878/1‑251 (MQ=255)                                                                                                                                                                                      tCCATATTACGTGGATAGGATCATAAAGTCATCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAAt                                                       >  2:380217/1‑251 (MQ=255)                                                                                                                                                                                         atatTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACCCTTCAATGTCATATTTAAATGAg                                                    <  2:12448/251‑1 (MQ=255)                                                                                                                                                                                                                aGTCCCCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTGTCTCAgg                             >  2:64372/1‑251 (MQ=255)                                                                                                                                                                                                                   ccTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGaa                                                                                                                                                                      >  1:216956/1‑111 (MQ=255)                                                                                                                                                                                                                   ccTCGGGTCGTTGTCC‑ACTCTTTTGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGaa                                                                                                                                                                      <  2:216956/111‑1 (MQ=255)                                                                                                                                                                                                                                         tttGGCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTGTCTCAGGAAAGTTGATTAATCTACACTcc       <  2:9740/249‑1 (MQ=255)                                                                                                                                                                                                                                            ggCATATCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTGTCTCAGGAAAGTTGATTAATCTACACTCCGTCtt  <  2:251907/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       AGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCC‑ACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGG‑TTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTGTCTCAGGAAAGTTGATTAATCTACACTCCGTCTT  >  NC_000913/316280‑316764

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