Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,486,205	C→A	T29K (ACG→AAG)	crp →	DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,205	0	C	A	96.2%	76.6 / NA	26	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (13/12);  total (13/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 = 1.00e+00

GTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCT‑ATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTG  >  NC_000913/3485962‑3486437                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           gTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTAtt                                                                                                                                                                                                                                     <  1:79450/251‑1 (MQ=255)      aCATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGCCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCNcc                                                                                                                                                                                                                                 >  2:260759/1‑250 (MQ=255)        aTTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAgg                                                                                                                                                                                                                              <  1:132164/251‑1 (MQ=255)                   gCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCCGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCCAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCACGAGCAAGCTTATTCCCCAGGGTGAAAAAGCg                                                                                                                                                                                                                   >  1:128724/1‑251 (MQ=255)                    cATGTATGCAAAGGACGTCACATTATCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGAGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGATAAAGCgg                                                                                                                                                                                                                  >  2:131387/1‑251 (MQ=255)                                         aTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAAAAGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAACAGGGTGAAAAAGCGGAAACGCTGTACTACATCGtt                                                                                                                                                                                              >  1:285754/1‑250 (MQ=255)                                             ccttGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGAACTCCAGAGACAGCGGCGTTATATGGCTCTGGAGAAAGTTTATAACAGAGGATAACAGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTATTGTCTCATTGCCACATTCATAAGAACCCATCCAAGAGCAAGCTTATCCACCAGGGTGAAAAAGCGGAAACGCTGAACTAAATCgatcaagg                                                                                                                                                                                         >  2:149159/4‑244 (MQ=255)                                              cGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAACAGCGGAAACGCTGTACTACATCGTTAAAGGc                                                                                                                                                                                        >  2:146304/1‑251 (MQ=255)                                               gTGCATTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGct                                                                                                                                                                                       <  1:227892/251‑1 (MQ=255)                                                         gTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGTCAATCCAGAGACAGCGGCGTTATCTGTCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGt                                                                                                                                                                             <  1:242208/251‑1 (MQ=255)                                                                       cccAGGAAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTAAAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACg                                                                                                                                                               <  2:112812/251‑1 (MQ=255)                                                                                    aaGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGCCGAAGAGGGTaa                                                                                                                                                     >  1:283544/1‑248 (MQ=255)                                                                                     atcATATGTCGCCAATCCAGAGACAGCGGCGTTAACTGGCTCTGAAGAAAGCTTAAAACAGAGGATAACCTCGCAGGTTGCTTGGCAAACCACGAACAGACCCGACTCTCGAATGGTTCGTGTCTCATTGCCACAGGCAGAAGTTGCTATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTTGGACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGaa                                                                                                                                                 <  1:113935/249‑1 (MQ=255)                                                                                            ttCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATcc                                                                                                                                          >  1:18614/1‑251 (MQ=255)                                                                                                     ccAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAACCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCT‑At                                                                                                                                 <  2:213095/250‑1 (MQ=255)                                                                                                     ccAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGGGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCT‑At                                                                                                                                 <  2:101222/250‑1 (MQ=255)                                                                                                                              gcggAGAAAGCTGATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACGCGACTCTCGAATGGCTCTTGTCTCATTGCTACATTCATAAGTACACATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGAGCTGATCAAAGACGAAGAGGGTAAAGACATGATCCTATCTT‑ATCTGAGTCAGGGTGATTTAATTGGCg                                                                                                        >  2:144324/3‑250 (MQ=255)                                                                                                                                          tataACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCT‑ATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGt                                                                                            <  2:1925/250‑1 (MQ=255)                                                                                                                                               cAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGc                                                                                                                                                                                                                    >  2:111630/1‑126 (MQ=255)                                                                                                                                               cAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGc                                                                                                                                                                                                                    <  1:111630/126‑1 (MQ=255)                                                                                                                                                           gcgcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAAGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCT‑ATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAAc                                                                          >  2:76454/1‑251 (MQ=255)                                                                                                                                                                                gCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCT‑ATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGT‑ACGTGCGa                                                    <  1:172619/251‑1 (MQ=255)                                                                                                                                                                                                   gAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCATAGGCGCAGAGGGTAAAGAAATGATCCTATCCT‑ATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGGTAAGTGGc                                 >  1:139811/1‑251 (MQ=255)                                                                                                                                                                                                               tctctttcccACATTCATAAGTACCCATCCAAGAACACTCTTTTTCACAGGGGTGAAAAAGGGGAAACCCTGTACTACATCTTTAAGGGCTTTGTGGCGGTGTTGAACAAAGACGAGGAGGTTAAAGAAATG‑TCCCCCCCTTAGCTAAATCAGGGTATTTTTATTGGCAAACTGGGCTTTTTTGAAGAGGTCCAGGAAGTTGGCTCAGGGTT‑CCGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGt                      <  1:139306/242‑1 (MQ=255)                                                                                                                                                                                                                                  agtccccaacaaaGAGCAATCTTATTCACCGGGGGGAAAAAGCGGAACCGCGGTACTCCAACTTAAAAGGCTCGGTGGCAGTGCGAAGAAAAGAGAAAGAGGGTAAAGAAAGGTTCCTCTCCT‑ATCTGAATCGGGGTGATTTTTTGGGCGAACTGGGCTTGTTTGAGGAGGGCCAGGAACGTAGCGC‑TGGGTTACGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTg  <  2:72399/240‑1 (MQ=255)                                                                                                                                                                                                                                  aGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCT‑ATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTg  >  1:282197/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           GTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCT‑ATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTG  >  NC_000913/3485962‑3486437

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