Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,340,953	C→T	G226R (GGG→AGG)	melR ←	DNA‑binding transcriptional dual regulator MelR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,340,953	0	C	T	85.3%	73.8 / 2.2	34	G226R (GGG→AGG)	melR	DNA‑binding transcriptional dual regulator MelR
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base T (14/15);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GATATCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTGCTCGTTCGGACT  >  NC_000913/4340710‑4341199                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       gataTCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGc                                                                                                                                                                                                                                                  <  1:200094/250‑1 (MQ=255)           ttAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGc                                                                                                                                                                                                                                         <  1:49866/249‑1 (MQ=255)                               ggcgCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGcc                                                                                                                                                                                                                   <  1:334382/251‑1 (MQ=255)                                         gCTAAGTTTGCGGTATTGTTGCTGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGa                                                                                                                                                                                                         <  1:315136/251‑1 (MQ=255)                                           tAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATg                                                                                                                                                                                                       >  1:118724/1‑251 (MQ=255)                                                    ggTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCTCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGcc                                                                                                                                                                                              >  1:94140/1‑251 (MQ=255)                                                          gtggcgtctaCATGCCGAAAGTTTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATTCTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCAGCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGCTCa                                                                                                                                                                                        <  2:119381/242‑1 (MQ=255)                                                               ggcgACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGtt                                                                                                                                                                                   >  1:333525/1‑251 (MQ=255)                                                                 cgACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGGAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGCCGTGCTCAGCCACATCGTTGATGGTCAGCGCCTACTCATAGtttt                                                                                                                                                                                 >  2:4936/1‑251 (MQ=255)                                                                           aCATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATaaa                                                                                                                                                                       <  1:116/251‑1 (MQ=255)                                                                                       aaCGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTg                                                                                                                                                           >  1:209041/1‑251 (MQ=255)                                                                                        acggtgtGTAAAAACGAATACTTTAACGAAAGCCGTCCGTCAGGGAAATATCGAGAATTCTTTTATCGGTAGCGCTCAGTAACGCGCGAACGTGGTTGATGCTCATCGCGGTAATGTACTGTTTCATCGTCAATGGCATGACCCGCTGGAATAGCCTCATTGCATAGTTGGCGGTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTgg                                                                                                                                                          <  2:235241/245‑1 (MQ=255)                                                                                          gTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCt                                                                                                                                                        >  1:137208/1‑251 (MQ=255)                                                                                                   aatcGACTACGTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGTCTAACATaaaa                                                                                                                                               <  1:214539/248‑1 (MQ=255)                                                                                                   aaaCGTCTATTTGAACGAAAGCCTGCCCGCAGGGCAATATCGGGAATTCTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGGTCATCGCTGTAATGTACTGTTTCATCGTTAATTGCATGATCATCTGTAATAGCCTCATTGCATAGTTGGCGTTAAGTGTGACTTGCTCAGCAACATCGTTGATGGTCAGCGCCTGATCATAGTTTACGGCAATAAAGCCCAGCATCTCGCTAACATaaaa                                                                                                                                               <  2:263959/251‑1 (MQ=255)                                                                                                    aaCGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATAGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCTCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATTGTTTTAGGCACTAAAGCCCAGAATCTGGCTAACATaaaa                                                                                                                                               >  2:336191/1‑250 (MQ=255)                                                                                                                     aaGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCATTGATGGTCAGCGCCTGATCATAGTTTTCGGCAACAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGAc                                                                                                                             >  1:89553/1‑251 (MQ=255)                                                                                                                        ccTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGAATGGCGCGAGACGCt                                                                                                                          >  1:150981/1‑251 (MQ=255)                                                                                                                              gtctgggCAATATCGAGTATCCTTTTAGCGGTATCTCCCTTTAACGCGCGAACGTGGTTGAGGCGCTTCGCGGTAATGTACTGTTTCAGCGTCAATTGCACGACCCGCTGGAATATCCCCATTGCATAGTTGTCGTTAAGTTTGACGGGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCTGCATCTGGCTAACATAAAATTGCGCATGGCGCGATACGCTGttttt                                                                                                                    <  2:94140/247‑1 (MQ=255)                                                                                                                                cAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTtgt                                                                                                                  <  1:236684/251‑1 (MQ=255)                                                                                                                                              gaATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGtt                                                                                                    >  1:82929/1‑251 (MQ=255)                                                                                                                                                                               gAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAgg                                                                   >  1:156394/1‑251 (MQ=255)                                                                                                                                                                                       ttGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCg                                                           >  1:128380/1‑251 (MQ=255)                                                                                                                                                                                                            tgctgtttcatcgtcatgtgcgattcccGCTGGAATATCCTCATTGCATATTTGGCGTTAAGTTTGTCTTGCACAGCCACATCGTTGATGTTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGTTAACATAAAATTGCGTATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAAttt                                       <  2:184924/225‑1 (MQ=255)                                                                                                                                                                                                                  ttCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTACGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGTCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCaa                                 <  1:96069/250‑1 (MQ=255)                                                                                                                                                                                                                   tCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATg                               >  1:76305/1‑251 (MQ=255)                                                                                                                                                                                                                     aTCGTCTATTGCATGACCCGCTGGAATATCCTCATGGCATAGTTGGCGGTAAGTTTGACTTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCTAATTTCATCAATGGc                             <  1:239229/251‑1 (MQ=255)                                                                                                                                                                                                                     aTCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGc                             >  1:26947/1‑251 (MQ=255)                                                                                                                                                                                                                               gCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCg                    <  1:336191/250‑1 (MQ=255)                                                                                                                                                                                                                               gCATGACCCGCGGGAATATCCTAATGTCATAGTGGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCg                    <  2:209041/250‑1 (MQ=255)                                                                                                                                                                                                                                 aTGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAAt                 >  1:128247/1‑251 (MQ=255)                                                                                                                                                                                                                                   gACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAAttt               <  1:317106/251‑1 (MQ=255)                                                                                                                                                                                                                                       cGCTTGAATATCCTCAGTGCATATTTGGCGTTAAGTTTGACTTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTGCTc           <  2:113514/251‑1 (MQ=255)                                                                                                                                                                                                                                                atCCCCATTGCAAAATTGGCGTTAAGTTTGACATGCTCAGCCACATCGTTTATGGTTAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCGGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTATTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGTTAAATCGCTTGAGCATCAGGCCAATTGCATCAATGGCGAGCTGGCGAATTTGCTCGTTCGGACt  <  2:116758/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       GATATCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTGCTCGTTCGGACT  >  NC_000913/4340710‑4341199

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