Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,707,947	C→A	intergenic (‑242/+669)	dppA ← / ← proK	dipeptide ABC transporter periplasmic binding protein/tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,707,947	0	C	A	95.7%	60.5 / ‑4.0	24	intergenic (‑242/+669)	dppA/proK	dipeptide ABC transporter periplasmic binding protein/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (10/12);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATAAGTAACGCAAAAAAAG  >  NC_000913/3707702‑3708173                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                   gCATTATTCTGGTCCAAGTGGGATGTTTGTTGTTTTAACCCTGTGCAGTGGGTTGTCGCTGCCTGGCCTTTTTGGTTTGTTGCCCGGTCGGGTTAACGTTATGAGGTGGGGATGCCTTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCGACAATCTGTCATCAGAATGTGAAAACGTCAATACAGCCCACCGGGATTTCCACCAACGTTGAGAATCCAcaaacaaa                                                                                                                                                                                                                               <  2:9424/251‑1 (MQ=255)           gCTCCAATTAAGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACACTCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAAAAAAGATTaaaca                                                                                                                                                                                                                      >  2:292898/1‑248 (MQ=255)             tCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAACATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAAc                                                                                                                                                                                                                   >  1:55028/1‑251 (MQ=255)                aaTTGTGATGTTTGTTGTTTTAACCCTTGTCAGTGGGTTGTCGCTGCCTGACCTTTTT‑GTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCTGATCCTCCATAAAATGCCCCTCGTGACCGACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCa                                                                                                                                                                                                               <  1:274298/251‑1 (MQ=255)                    gtgATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAAc                                                                                                                                                                                                            >  1:152736/1‑251 (MQ=255)                        tgtttgttGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACTGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCACCGGTGAGAATCCACAAACAAAGATTACAAAAACTTCAAACAGCt                                                                                                                                                                                                        >  2:91830/1‑251 (MQ=255)                                   ttAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACCAAGATTAAAAAAACTTCAAACAGCTATTTgcagca                                                                                                                                                                                              >  1:173413/1‑250 (MQ=255)                                      atccTTTGCAGTGGGTTGTCGCTGCCTGATCTTTTTTGGTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAAGGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAgaga                                                                                                                                                                                          <  1:70433/249‑1 (MQ=255)                                        ccTTTGCAGTGGGTTGTCGCTGCCTTACCGTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTTTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGaaa                                                                                                                                                                                        <  2:141322/251‑1 (MQ=255)                                             gCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtgt                                                                                                                                                                                   <  2:75819/251‑1 (MQ=255)                                                                              ttttgCCCGGTCGGGTCAACGTTATGAGGCGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTa                                                                                                                                                   >  1:232879/1‑250 (MQ=255)                                                                                                       gagggggtgATGCAGTACTAATTAACATCAGTGTAACAGTCCGGAGACCACATAAAATGCTCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAAAGTCAATACAGTCATCCTGGATTTACACCAATGGTGAGAATCCACAAACAAAGAGTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAAGATGACCAGAACAATCAGCGTAATATTTAGCAACATTTGGTGATTAAg                                                                                                                         <  2:70372/243‑1 (MQ=255)                                                                                                                       aCTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTa                                                                                                                           >  1:232782/1‑233 (MQ=255)                                                                                                                       aCTAATTAACATCAGAGTAACAGTGCGGATCCTCCATAAAATGCCCCTAGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAACAAACTTACAACAGCTATTTGCAGCAGAGAAATTTGTGCTACACCCACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTAAATg                                                                                                         >  2:219445/1‑251 (MQ=255)                                                                                                                       aCTAAATAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTa                                                                                                                           <  2:232782/233‑1 (MQ=255)                                                                                                                           aTTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGAGTAAAAAAGCTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGaaaa                                                                                                     <  1:268001/251‑1 (MQ=255)                                                                                                                             aaaCATCCGTTTAACAGTGCGAATCATCCATAAAATGCCCCTCATCACCTACAAGCTGTCAACAGAATGTGAAAACGTCAAGACAGACAAACGAGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACAAGATCAATCAGCTTAATATTTAGCAACACTTGGTGACTAATGTTTATGTATTTAATGaaaaaa                                                                                                   <  2:125178/250‑1 (MQ=255)                                                                                                                                   cAGTGTAACAGTGCGGATCCGCCATAAAATGCCCCTCGTGACCTACAATCTTTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGATAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCATAACAATCAGCTTAATATTTAGCAACATTTGTTTATTAATTTTTATGCATTTAATGAAATAAATCTg                                                                                              <  2:37941/250‑1 (MQ=255)                                                                                                                                                                               aCAATCTGGCACAAGAATGTGAAAACGTCAATACAGCCAACCGGGATGTACCCCAACGGTGAGAATCCACAATCAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCATAGAAATTTTTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTGATGCATTTAATGAAAAAAAGCTGAGGAAAAGGTGGATATCTGAGCTATTAATGCCATGATCGCTaaa                                                  <  2:273812/250‑1 (MQ=255)                                                                                                                                                                                             gAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGTaaaa                                   <  2:153474/251‑1 (MQ=255)                                                                                                                                                                                                  tgAAAAAGTCAATACAGACAACCGGGATTTACACAAACGGTGAGAATCCACAAACAAAGATTAAAAAACCTTAAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTCAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAAAACCATGATCGCTAAAACGAACAGCGGAACATGca                               >  2:205617/1‑249 (MQ=255)                                                                                                                                                                                                                 aGCCAACCGGGATTTACACCAACGGTGAGAAACCACAAACAAAGATTAAAAAAACTACAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCAGTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCGAAAACGAACAGCGGAAAATGCCAGGGCCAGCCATAAGt               <  2:94011/251‑1 (MQ=255)                                                                                                                                                                                                                  gCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTACTGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATAAGt               >  2:235580/1‑250 (MQ=255)                                                                                                                                                                                                                      aCCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATAAGTAACGc          <  1:91830/251‑1 (MQ=255)                                                                                                                                                                                                                               tACACCAACGGTGAGAATCAACAGACAAAGATTAAAGAAACTTAAAACAGCTATTTGCAGAAGAGAAAATTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGCTTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATAAACGCTAGAACGCACAGCGGAAATTGCATGAGCCAACACTAAGTAACGCaaaaaacg  >  2:228623/1‑248 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                   GCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATAAGTAACGCAAAAAAAG  >  NC_000913/3707702‑3708173

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