Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	70,581	C→T	V65V (GTC→GTT)	araC →	DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	70,581	0	C	T	92.0%	74.3 / NA	25	V65V (GTC→GTT)	araC	DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/8);  minor base G (0/2);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGG  >  NC_000913/70338‑70821                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 gCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGAc                                                                                                                                                                                                                                           <  1:311188/251‑1 (MQ=255)       tGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGAGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGGAAGAAAGGTTATATTTTCAATCTCACCATTCGCGGGTAGGGGGGGGGGAAAAATCAGGGACGAGAATTTGTTTTCCGACCggggg                                                                                                                                                                                                                                     >  1:194830/1‑249 (MQ=255)                    tcttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGt                                                                                                                                                                                                                         <  2:21733/250‑1 (MQ=255)                                             aaaGTATGGCTGAAGCGCAAAATGTTCCCCTGCGGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACGGCTGGGAATGAAAGGTTATATTCTCAATCTCCCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTGTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTAc                                                                                                                                                                                               <  2:241955/251‑1 (MQ=255)                                                                                  gggATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGt                                                                                                                                                          >  1:365250/1‑251 (MQ=255)                                                                                           gTTTAATGCACATCTGGGGGCCGGTTTAACGCCAATTGAGGCCAACGGTTACCTCGATTTTTTTATCCACACAAAGCTGGGCATGAAAGGTTATATTCTCAATCTAACCATTCGCGGCCAGGGGGTGGTGAAAACTCAGGGACGAGAATTTTTTTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCCCTACGGTCGTCATCCTGAGGCTAGCGAATGGTATCACCAGTGGGTTTACt                                                                                                                                                 >  2:244974/1‑251 (MQ=255)                                                                                                    ccATCTGGTGGCGGGGTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGCAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCcgc                                                                                                                                        >  1:236637/1‑251 (MQ=255)                                                                                                                    ttAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGGATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATg                                                                                                                         >  2:9572/1‑250 (MQ=255)                                                                                                                      aaCGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAAt                                                                                                                      >  1:156838/1‑251 (MQ=255)                                                                                                                       ccgccaatagAGGCCAACGGTTATCTCAATTTTTTTATCCACCGACCGCTGGTAATGAAAGGTGAGATCCTCAATCTCCCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGAGGAGATTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGTGAGAGTCATCACCACGGTCGTCATCCGGAGGCTGGCGAATGGTATCACCAGGGGGTTTACGTTCTCCCGCGCTGCTTCTGGCATGAATg                                                                                                                     <  2:183673/242‑1 (MQ=255)                                                                                                                                gAGGCAAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTg                                                                                                            >  2:216782/1‑251 (MQ=255)                                                                                                                                  ggCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCATCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGc                                                                                                          >  1:127755/1‑251 (MQ=255)                                                                                                                                                             cgaccgacCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttct                                                                               >  1:335085/1‑251 (MQ=255)                                                                                                                                                               accgaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGAGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCAAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttcttt                                                                             >  1:352952/1‑251 (MQ=255)                                                                                                                                                                  gcccgCTGGGAATGAAATGTTATATTGTCAATCTCACCATTCGCGGTAAGGGGGTGGTTAAAAATCAGGGACGAGAATTTGTTTGCAGCCCGGGTGATATTTTTCTTTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGc                                                                          <  2:278174/249‑1 (MQ=255)                                                                                                                                                                  gaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGc                                                                          >  1:99460/1‑251 (MQ=255)                                                                                                                                                                             aTGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCGGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAgc                                                               <  2:365250/251‑1 (MQ=255)                                                                                                                                                                                                   tcACCGTTCGCGGTCAGGGGGTGGAGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGCTATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACGGGCCGTCATTATTTGCCAGTACGGGGTTCTTTCGCCCGGATGACGCGCACCAGCCGCATTTCAGCGAc                                         >  2:4794/1‑251 (MQ=255)                                                                                                                                                                                                    cACCATTCGCGGTCAGGGGGTGGTGAAAAAGCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGc                                            <  1:231181/247‑1 (MQ=255)                                                                                                                                                                                                              ggTCAGGGGGTGGGGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCa                              <  1:160335/251‑1 (MQ=255)                                                                                                                                                                                                                 cAGGGGGTGGTGCCAAATCAGGGACGAGCATTTGTTTGCCGACCGGGTGATATTTTGCTGTTACCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAAt                           >  2:173942/1‑251 (MQ=255)                                                                                                                                                                                                                  aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATc                          >  1:161567/1‑251 (MQ=255)                                                                                                                                                                                                                  aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATc                          >  1:139432/1‑251 (MQ=255)                                                                                                                                                                                                                                          gagaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAggg  <  1:125807/251‑1 (MQ=255)                                                                                                                                                                                                                                          gagaATTGGTGTGCCGCCCGGGTGATATTTGGCTTTTCCCGCCAGGAGAGATGCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAggg  <  1:319498/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 GCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGG  >  NC_000913/70338‑70821

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