Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,828,529	T→C	G402G (GGA→GGG)	astD ←	aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,828,529	0	T	C	95.0%	58.9 / ‑3.8	20	G402G (GGA→GGG)	astD	aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (8/11);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

CCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGT  >  NC_000913/1828407‑1828660                                                                                                                           |                                                                                                                                    ccgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTAtt                                                                                                                    >  2:54217/1‑140 (MQ=255)                 tGGCAGCACCGGTAAGCCGTTTATTCCAGTTCACACTCCACGCCCGCGACTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCt                                                                                                    <  2:196921/140‑1 (MQ=255)                        aCCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGa                                                                                             <  1:621790/140‑1 (MQ=255)                                     ttGTTCCAGTTAACAATCCCCGGCCGCGCTTCCAGCACCAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACg                                                                                <  1:298165/140‑1 (MQ=255)                                            aGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGAcg                                                                         <  2:464477/140‑1 (MQ=255)                                              ttAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGAcgcg                                                                       >  1:571629/1‑140 (MQ=255)                                                          gcccgcGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTcc                                                           >  1:386531/1‑140 (MQ=255)                                                             cgcgCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGaa                                                        <  2:571629/140‑1 (MQ=255)                                                              gcgcCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAac                                                       <  2:392471/140‑1 (MQ=255)                                                               cgcCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAGGTATCATGACGACAGACGCGCCATAACGGTCCGAaca                                                      >  1:274263/1‑140 (MQ=255)                                                               cgcCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAaca                                                      >  2:29711/1‑140 (MQ=255)                                                                      aGCAACAGTTGATCGAACTTTTTCCTCTCCGGGGAAACCAGACCGCAAGAGAGCCCGAAGCG‑GTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCa                                              <  1:225438/140‑1 (MQ=255)                                                                            aGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTgg                                         <  1:511882/140‑1 (MQ=255)                                                                              ttGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTa                                       >  1:611529/1‑140 (MQ=255)                                                                                     aaCTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGc                                >  1:632749/1‑140 (MQ=255)                                                                                            cccGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCt                         <  1:54217/140‑1 (MQ=255)                                                                                              cGCTCGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGt                       <  1:474206/140‑1 (MQ=255)                                                                                                  cGGGGGAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTATCAtt                   <  1:548922/140‑1 (MQ=255)                                                                                                    ggggAAACCAGACCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTc                 >  2:331771/1‑140 (MQ=255)                                                                                                               aCCGCAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCAGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCg      >  1:263796/1‑140 (MQ=255)                                                                                                                   cAAGAGAGCCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTTCGCCAGCAACGACTGTCATTTGAATGATCCCCGGCGt  >  1:278513/1‑140 (MQ=255)                                                                                                                           |                                                                                                                                    CCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGT  >  NC_000913/1828407‑1828660

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