Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,337,446	C→A	D162Y (GAC→TAC)	adiY ←	DNA‑binding transcriptional activator AdiY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,337,446	0	C	A	93.8%	47.7 / ‑5.7	16	D162Y (GAC→TAC)	adiY	DNA‑binding transcriptional activator AdiY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (6/9);  minor base T (1/0);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATGCTGTGTTGCGCGGCAACGCGAAAGACGTCTGGCGTCAGGCATGGCAGGGTGAGTATGGGCGTAGCGCTACGTTG  >  NC_000913/4337209‑4337689                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    gcggcgTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCg                                                                                                                                                                                                                                        <  2:166868/251‑1 (MQ=255)   ggcgTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCGCt                                                                                                                                                                                                                                      <  1:262859/251‑1 (MQ=255)                 gacgTCTTTAAAGGCAGAAATAATGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCCGGTGGCGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCCTCCTAAGATTCTAGTTTTTGTGAATATCGCTTTCAATAATTTGAt                                                                                                                                                                                                                        >  2:179789/4‑251 (MQ=255)                           aaGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTc                                                                                                                                                                                                              <  2:36166/251‑1 (MQ=255)                                                                                 gagaTATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGAGTCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACATAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTc                                                                                                                                                        <  1:36123/251‑1 (MQ=255)                                                                                  agaTATTTTTACCGTCCATCATTAATTCATGTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACGAAGAGTCCATTATTTGTGAATATCTGTTTCAATAAGTTGAGAAGCGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCg                                                                                                                                                       <  1:153193/251‑1 (MQ=255)                                                                                                         aagtCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGGTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGGTTCCAGTAATTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTATCATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGaa                                                                                                                                <  1:180547/248‑1 (MQ=255)                                                                                                                  aCGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTTAATAATAATgc                                                                                                                       <  2:174712/251‑1 (MQ=255)                                                                                                                    ggCATAACGCATGCGGCAGGTGGTGATTATTTGGCTTTAACTGGTGTTTTCGTTTTTCACTTTTTTCTTTAACAGACTTGGGCTAAGCAATAAACACCTGGCTACCATCCTAGGATTCCAGTATTTGTGATTATCCCTTTCAAAAATTTGATAAACGCTGTCACTTACCCAATTACGAAACTTATACACCATTAATGAAACGAATTTTTTACTGTCGAAAAAACGGGAAAGCACGGTGattataaatccgg                                                                                                                     >  1:285287/1‑239 (MQ=255)                                                                                                                                     aGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCt                                                                                                    >  1:326141/1‑251 (MQ=255)                                                                                                                                          gtgATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTc                                                                                               >  1:5791/1‑251 (MQ=255)                                                                                                                                           tgATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGAATTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCACTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTc                                                                                               >  1:159357/1‑250 (MQ=255)                                                                                                                                                attTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTc                                                                                         >  2:305036/1‑251 (MQ=255)                                                                                                                                                  ttGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTGTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGGGTTCGTTTCTTTTCTTGCTCATTTTCt                                                                                        <  2:40596/250‑1 (MQ=255)                                                                                                                                                                                                                                    aGATTCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATGCTGTGTTGCGCGGCAACGCGAAAGACGTCTGGCGTCAGGCATGGCAGGGTGAGTATGGGCGTAGCGCTACg     >  1:83153/1‑251 (MQ=255)                                                                                                                                                                                                                                        tCCAGTATTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATGCTGTGTTGCGCGGCAACGCGAAAGACGTCTGGCGTCAGGCATGGCAGGGTGAGTATGGGCGTAGCGCTACGTTg  <  1:179789/250‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    GCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGTACTGTTGTAGCCGCAGGACTGTGATACCTGAGAGATATTTTTACCGTCCATCATTAATTCATTTACGGCATAACGCATGCGGCAGGTGGTGATTATTTGGCTATAACTGGTGTTTTCGCTTTTCAACTTTTTCTTTAACAGACTTGGGCTAAGACATAAACAACTGGCTACCATACTAAGATTCCAGTCTTTGTGAATATCGCTTTCAATAATTTGATAAACGCTGTCACTTACACAATTACGTAACATATACATCATTAATGAAACGAATTTTTTACTGTCGAGAAAACGGGATAGCACCGTGAATAATAATGCGCGTGTTCGTTCCTTTTCTGACTCAGTTTCTGCGGGCATCATGCTGTGTTGCGCGGCAACGCGAAAGACGTCTGGCGTCAGGCATGGCAGGGTGAGTATGGGCGTAGCGCTACGTTG  >  NC_000913/4337209‑4337689

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