Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,489,590	T→C	R453R (CGT→CGC)	aldA →	aldehyde dehydrogenase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,489,590	0	T	C	96.3%	187.7 / ‑0.0	54	coding (1359/1440 nt)	aldA	aldehyde dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (25/27);  minor base . (1/1);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

TTACGGCCTGACCTCATCAATCTATACCCAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGTAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTT  >  NC_000913/1489455‑1489721                                                                                                                                        |                                                                                                                                    ttACGGCCTGACCTCATCAATCTATACCCAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAAt                                                                                                                                 <  1:183009/140‑1 (MQ=255)         tGACCTCATCAATCTATACCCAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTG‑AACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTg                                                                                                                        >  1:342716/1‑140 (MQ=255)                 tcaATCTATACCCAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCG‑CAATCCGGTATTGGCGGCGCa                                                                                                                >  1:749165/1‑140 (MQ=255)                  caATCTATACCCAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCa                                                                                                                <  1:178385/140‑1 (MQ=255)                     tCTATACCCAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGAt                                                                                                             >  2:456381/1‑140 (MQ=255)                     tCTATACCCAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTG‑AACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATg                                                                                                            <  2:342716/140‑1 (MQ=255)                          aCCCAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTaa                                                                                                        >  1:155857/1‑140 (MQ=255)                           cccAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGgcgc                                                                                                                 <  2:794678/130‑1 (MQ=255)                           cccAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGgcgc                                                                                                                 >  1:794678/1‑130 (MQ=255)                            ccAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAAc                                                                                                      <  1:251081/140‑1 (MQ=255)                              aaaaTCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACAt                                                                                                    >  1:323795/1‑140 (MQ=255)                               aaaTCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATg                                                                                                   >  2:884725/1‑140 (MQ=255)                               aaaTCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATg                                                                                                   >  2:763921/1‑140 (MQ=255)                               aaaTCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCG‑CAATCCGGTATTGGCGGCGCAGATGGTAAACATgg                                                                                                  <  2:749165/140‑1 (MQ=255)                                aaTCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATgg                                                                                                  <  2:266062/140‑1 (MQ=255)                                 aTCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGc                                                                                                 <  1:677336/140‑1 (MQ=255)                                   cTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCtt                                                                                               >  1:494816/1‑140 (MQ=255)                                      aaCGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCa                                                                                            <  2:155857/140‑1 (MQ=255)                                      aaCGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCa                                                                                            >  2:796490/1‑140 (MQ=255)                                           cgcgATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAat                                                                                       >  2:589129/1‑140 (MQ=255)                                                  aaaGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCa                                                                                <  1:456381/140‑1 (MQ=255)                                                  aaaGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCa                                                                                <  2:136725/140‑1 (MQ=255)                                                    aGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGa                                                                              >  2:728318/1‑140 (MQ=255)                                                            aaGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAggtg                                                                      >  2:854338/1‑140 (MQ=255)                                                             aGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAggtgg                                                                     <  1:854338/140‑1 (MQ=255)                                                                 cTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGAc                                                                             >  2:91553/1‑128 (MQ=255)                                                                 cTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGAc                                                                             <  1:91553/128‑1 (MQ=255)                                                                 cTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAggt                                                                       <  2:112540/134‑1 (MQ=255)                                                                 cTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAggt                                                                       >  1:112540/1‑134 (MQ=255)                                                                 cTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGttta                                                                 >  1:659417/1‑140 (MQ=255)                                                                 cTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGttta                                                                 >  1:631844/1‑140 (MQ=255)                                                                 cTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGttta                                                                 >  1:314793/1‑140 (MQ=255)                                                                     aGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGtttattta                                                             <  2:450725/140‑1 (MQ=255)                                                                     aGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGtttattta                                                             <  2:659417/140‑1 (MQ=255)                                                                     aGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGtttattta                                                             <  1:796490/140‑1 (MQ=255)                                                                     aGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGtttattta                                                             <  1:13882/140‑1 (MQ=255)                                                                               aaCTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAAt                                                   <  2:314793/140‑1 (MQ=255)                                                                               aaCTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAAt                                                   <  1:589129/140‑1 (MQ=255)                                                                                           ccGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAgg                                       >  2:813542/1‑140 (MQ=255)                                                                                           ccGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAgg                                       >  1:520225/1‑140 (MQ=255)                                                                                              tGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCgg                                    >  1:117902/1‑140 (MQ=255)                                                                                              tGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCgg                                    <  2:494816/140‑1 (MQ=255)                                                                                              tGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCgg                                    >  1:900237/1‑140 (MQ=255)                                                                                                              aTGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGc                    <  2:117902/140‑1 (MQ=255)                                                                                                              aTGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGc                    <  2:900237/140‑1 (MQ=255)                                                                                                               tGCAAGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGGTTTTTCCTCCGcc                   >  2:793201/1‑140 (MQ=255)                                                                                                                   aGGCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCtgtg               <  2:396634/140‑1 (MQ=255)                                                                                                                     gCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTgcg             <  2:520225/140‑1 (MQ=255)                                                                                                                     gCTTCCACGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTgcg             <  1:884725/140‑1 (MQ=255)                                                                                                                            cGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGcgcg           >  2:716923/1‑135 (MQ=255)                                                                                                                            cGCCGGATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGcgcg           <  1:716923/135‑1 (MQ=255)                                                                                                                                ggATGGCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGttt  >  2:395879/1‑140 (MQ=255)                                                                                                                                    ggCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTc         >  1:900881/1‑129 (MQ=255)                                                                                                                                    ggCGCAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTc         <  2:900881/129‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    TTACGGCCTGACCTCATCAATCTATACCCAAAATCTGAACGTCGCGATGAAAGCCATTAAAGGGCTGAAGTTTGGTGAAACTTACATCAACCGTGAAAACTTCGAAGCTATGCAAGGCTTCCACGCCGGATGGCGTAAATCCGGTATTGGCGGCGCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTT  >  NC_000913/1489455‑1489721

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