Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,222,536	A→C	intergenic (‑49/+268)	arpA ← / ← iclR	ankyrin repeat protein/transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,222,536	0	A	C	100.0%	147.7 / NA	40	intergenic (‑49/+268)	arpA/iclR	ankyrin repeat protein/transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/19);  total (21/19)

TTCTCTGAGAAGCTTTTTAATAGAGGCGTCGCCAGCTCCTTGCCAGAAAATTTATCCTCGAGTTCTTTATAAAACAATTCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACAGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACCCTTATTC  >  NC_000913/4222327‑4222750                                                                                                                                                                                                                  |                                                                                                                                                                                                                       ttCTCTGAGAAGCTTTTTAATAGAGGCGTCGCCAGCTCCTTGCCAGAAAATTTATCCTCGAGTTCTTTATAAAACAATTCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAAc                                                                                                                                                                                              >  1:390882/1‑236 (MQ=255)        agaAGCTTTTTAATAGAGGCGTCGCCAGCTCCTTGCCAGAAAATTTATCCTCGAGTTCTTTATAAAACAATTCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGaa                                                                                                                                                                                       <  2:142642/236‑1 (MQ=255)           aGCTTTTTAATAGAGGCGTCGCCAGCTCCTTGCCAGAAAATTTATCCTCGAGTTCTTTATAAAACAATTCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTTATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAATTAAAGAAata                                                                                                                                                                                    >  2:59673/1‑236 (MQ=255)                          gCGTCGCCAGCTCCTTGCCAGAAAATTTATCCTCGAGTTCTTTATAAAACAATTCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGaaa                                                                                                                                                                                                                  >  2:418229/1‑191 (MQ=255)                          gCGTCGCCAGCTCCTTGCCAGAAAATTTATCCTCGAGTTCTTTATAAAACAATTCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGaaa                                                                                                                                                                                                                  <  1:418229/191‑1 (MQ=255)                             tCGCCAGCTCCTTGCCAGAAAATTTATCCTCGAGTTCTTTATAAAACAATTCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAAtct                                                                                                                                                                  >  1:594220/1‑236 (MQ=255)                                                                          tcAATTCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGAtt                                                                                                                     <  1:756470/235‑1 (MQ=255)                                                                               tCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGtt                                                                                                                <  2:709287/236‑1 (MQ=255)                                                                                     aGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTaataat                                                                                                          >  1:618031/1‑236 (MQ=255)                                                                                              tgtTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGt                                                                                                  <  1:354343/235‑1 (MQ=255)                                                                                                 tCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGa                                                                                              >  2:153714/1‑236 (MQ=255)                                                                                                        tCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTaa                                                                                                       <  1:677347/220‑1 (MQ=255)                                                                                                        tCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTaa                                                                                                       >  2:677347/1‑220 (MQ=255)                                                                                                             gCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAaat                                                                                   >  1:325458/1‑235 (MQ=255)                                                                                                                   ttattaATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCt                                                                             <  2:594220/235‑1 (MQ=255)                                                                                                                      ttaATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGt                                                                         >  1:28657/1‑236 (MQ=255)                                                                                                                           atTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTtaata                                                                    >  1:156946/1‑236 (MQ=255)                                                                                                                           atTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTtaata                                                                    >  1:651061/1‑236 (MQ=255)                                                                                                                              tGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGa                                                                           >  2:87860/1‑226 (MQ=255)                                                                                                                              tGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGa                                                                           <  1:87860/226‑1 (MQ=255)                                                                                                                               gCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAg                                                                <  2:651061/236‑1 (MQ=255)                                                                                                                                             aCGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGtgtgtg                                                          <  1:153714/228‑1 (MQ=255)                                                                                                                                             aCGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGaa                                                        <  2:390882/230‑1 (MQ=255)                                                                                                                                                           tAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTAtt                                    >  2:464953/1‑236 (MQ=255)                                                                                                                                                                 tttAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCg                              >  2:168743/1‑236 (MQ=255)                                                                                                                                                                          ctcATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCa                           >  2:456051/1‑230 (MQ=255)                                                                                                                                                                          ctcATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCa                           <  1:456051/230‑1 (MQ=255)                                                                                                                                                                           tcATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGAt                    <  2:156946/236‑1 (MQ=255)                                                                                                                                                                                  ttATTGATGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGa             <  2:325458/236‑1 (MQ=255)                                                                                                                                                                                         tGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTaa                                                                                       >  2:231666/1‑155 (MQ=255)                                                                                                                                                                                         tGTCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTaa                                                                                       <  1:231666/155‑1 (MQ=255)                                                                                                                                                                                           tCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGaa                                                        >  2:40133/1‑184 (MQ=255)                                                                                                                                                                                           tCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGaa                                                        <  1:40133/184‑1 (MQ=255)                                                                                                                                                                                           tCTTGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGCAACCCTTAt    >  2:218081/1‑236 (MQ=255)                                                                                                                                                                                             ttGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCgg                                                                                                                                         >  2:405497/1‑101 (MQ=255)                                                                                                                                                                                             ttGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCgg                                                                                                                                         <  1:405497/101‑1 (MQ=255)                                                                                                                                                                                             ttGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGacac          <  2:180428/228‑1 (MQ=255)                                                                                                                                                                                             ttGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGacac          >  1:180428/1‑228 (MQ=255)                                                                                                                                                                                             ttGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACCCTTATTc  >  1:431915/1‑236 (MQ=255)                                                                                                                                                                                             ttGAATTTTAACTATCAGAACCGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACCCTTATTc  <  1:464953/236‑1 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                                       TTCTCTGAGAAGCTTTTTAATAGAGGCGTCGCCAGCTCCTTGCCAGAAAATTTATCCTCGAGTTCTTTATAAAACAATTCACTCAGGGTTTGGTGTTCATTTGTCTGGGCTGTATTATTAATATTTGCAGAGAAAGAACTACGAGGAATACGAGTAATCATTTAAATATCTCATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACAGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGATCTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACCCTTATTC  >  NC_000913/4222327‑4222750

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