Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,843,339	C→T	T323T (ACC→ACT)	gdhA →	glutamate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,843,339	0	C	T	96.7%	96.0 / NA	30	T323T (ACC→ACT)	gdhA	glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/12);  minor base G (0/1);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGTCGCCCAGTACGCTATCGAAAAAGCGATGGAATTTGGTGCTCGTGTGATCACTGCGTCAGACTCCAGCGGCACTGTAGTTGATGAAAGCGGATTCACGAAAGAGAAACTGGCACGTCTTATCGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACCCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGTTGACGCACGTTTGCATC  >  NC_000913/1843096‑1843580                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  cGTCGCCCAGTACGCTATCGAAAAAGCGATGGAATTTGGTGCTCGTGGGATCACTGCGTCAGACTCCAGCGGCACTGTAGTTGATGAAAGCGGATTCACGAAAGAGAAACTGGCACGTCTTATCGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATg                                                                                                                                                                                                                                            >  1:171320/1‑251 (MQ=255)                   cGAAAAAGCGATGGAATTTGGTGCTCGTGTGATCACTGCGTCAGACTCCAGCGGCACTGTAGTTGATGAAAGCGGATTCACGAAAGAGAAACTGGCACGTCTTATCGAAATCAAAGCCAGCCGCGATGGGCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGAcgcc                                                                                                                                                                                                                           <  2:183658/250‑1 (MQ=255)                         aGCGATGGAATTTGGTGCTCGTGTGATCACTGCGTCAGACTCCAGCGGCACTGTAGTTGATGAAAGCGGATTCACGAAAGAGAAACTGGCACGTCTTATCGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATc                                                                                                                                                                                                                    >  2:193185/1‑251 (MQ=255)                                gAATTTGGTGCTCGTGTGATCACTGCGTCAGACTCCAGCGGCACTGTAGTTGATGAAAGCGGATTCACGAAAGAGAAACAGGCACGTCTTATCGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTAt                                                                                                                                                                                                             <  2:293707/251‑1 (MQ=255)                                                      ctgcgcccgACTCCAGCGGCACTGTAACTGATGAACGCGAATTCAAGAAAGAGAAACTGCCACGTCTTATCGAAATCAAAGACAGCCGCGCTGGTAGAGTGGCAGCTTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGACACTCAGAATGAACTGGATGTTGACGCCTCGAATCAGCTTATCGCTAATGGCGTTAAAGCNga                                                                                                                                                                                        >  2:372251/9‑249 (MQ=255)                                                                   cAGCGGCACTGTAGTTGATGAAAGCGGATTCACGAAAGAGAAACTGGCACGTCTTATCGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCa                                                                                                                                                                           >  1:367427/1‑250 (MQ=255)                                                                             gTAGTTGATGAAAGCGGATTCACGAAAGAGAAACTGGCACGTCTTATCGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCg                                                                                                                                                                  <  1:308988/249‑1 (MQ=255)                                                                                    attaaaGCGGATTCACGAAAGAGAACCTGGCACGTCTTATCGAAATCAAAGCCAGCCGCTATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTGGCGCCACGCAGAATGAACTGGATGTTGACGCCGCGCAGCAGTTTAGCGCTAATGGCGTTAAAGCCGTCGCCGAAGTGCCAAATATGCCGACCACCATc                                                                                                                                                         <  2:184166/248‑1 (MQ=255)                                                                                                  aCGAAAGAGAAACTGGCACGTCTTATCGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACt                                                                                                                                           <  2:132206/251‑1 (MQ=255)                                                                                                                     gTCTTATCGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTACGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGt                                                                                                                         >  1:223946/1‑250 (MQ=255)                                                                                                                            cGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTAttt                                                                                                                  <  1:86595/250‑1 (MQ=255)                                                                                                                             gAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGc                                                                                                                >  1:132926/1‑251 (MQ=255)                                                                                                                               aaTCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCAc                                                                                                              >  2:277247/1‑251 (MQ=255)                                                                                                                                               gATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGACTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTa                                                                                               >  2:156842/1‑250 (MQ=255)                                                                                                                                                              gATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGACGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCGAATGCTGGTGGCGTCGc                                                                               >  2:137838/1‑251 (MQ=255)                                                                                                                                                               aTTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGCCGCt                                                                              >  2:214103/1‑251 (MQ=255)                                                                                                                                                                    gCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATc                                                                         <  2:182084/251‑1 (MQ=255)                                                                                                                                                                      cAAATAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCgg                                                                       <  2:254332/251‑1 (MQ=255)                                                                                                                                                                        aaGAATTTGGTCTGGACTATCTCGAAGGCCAACAGCCGTGGTCTATACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGACGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAGGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGCCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCgtgc                                                                     >  2:410715/1‑248 (MQ=255)                                                                                                                                                                                ggtctggtctATCTCGAAGCCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGTGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGTTAAAGCGGCTAATGCTGGTGGCGTGGCTACCTCGGGCCTGGaaa                                                              <  2:327049/250‑1 (MQ=255)                                                                                                                                                                                              cGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGAGGTGGACGCCGCGCATCAGCTTATCGGTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCAACGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGcg                                                <  2:77450/250‑2 (MQ=255)                                                                                                                                                                                                                            tGATATCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGtt                  <  2:367427/250‑1 (MQ=255)                                                                                                                                                                                                                              atatCGCCCTGCCTTGCGCCACTCAGAATGAACTGTATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTAGTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGTTGAc               <  1:410047/251‑1 (MQ=255)                                                                                                                                                                                                                              atatCGCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGTTGAc               >  1:182563/1‑251 (MQ=255)                                                                                                                                                                                                                                   gcacTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATAGGGCCTGGAAATGGCACAAAACTCTGCGAGCCTGGGCTGGAAAGCCGAGAAAGTTGACGCACg          >  1:203840/4‑251 (MQ=255)                                                                                                                                                                                                                                   gCCCTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGTTGACGCACg          >  1:203831/1‑251 (MQ=255)                                                                                                                                                                                                                                     ccTGCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGTTGACGCACGtt        <  2:85491/251‑1 (MQ=255)                                                                                                                                                                                                                                       tGCCTTGCGCCACTAAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCCCAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGTTGACGCACGTTTg      >  1:357953/1‑251 (MQ=255)                                                                                                                                                                                                                                        gCCTTGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGTTGACGCACGTTTGc     >  1:117674/1‑251 (MQ=255)                                                                                                                                                                                                                                            tGCGCCACTCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGAGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGTTGACGCACGTTTGcatc  >  1:308703/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  CGTCGCCCAGTACGCTATCGAAAAAGCGATGGAATTTGGTGCTCGTGTGATCACTGCGTCAGACTCCAGCGGCACTGTAGTTGATGAAAGCGGATTCACGAAAGAGAAACTGGCACGTCTTATCGAAATCAAAGCCAGCCGCGATGGTCGAGTGGCAGATTACGCCAAAGAATTTGGTCTGGTCTATCTCGAAGGCCAACAGCCGTGGTCTCTACCGGTTGATATCGCCCTGCCTTGCGCCACCCAGAATGAACTGGATGTTGACGCCGCGCATCAGCTTATCGCTAATGGCGTTAAAGCCGTCGCCGAAGGGGCAAATATGCCGACCACCATCGAAGCGACTGAACTGTTCCAGCAGGCAGGCGTACTATTTGCACCGGGTAAAGCGGCTAATGCTGGTGGCGTCGCTACATCGGGCCTGGAAATGGCACAAAACGCTGCGCGCCTGGGCTGGAAAGCCGAGAAAGTTGACGCACGTTTGCATC  >  NC_000913/1843096‑1843580

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