Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,107,146	T→C	Q670Q (CAA→CAG)	speC ←	ornithine decarboxylase, biosynthetic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,107,146	0	T	C	100.0%	78.9 / NA	24	Q670Q (CAA→CAG)	speC	ornithine decarboxylase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/10);  total (14/10)

TCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGTTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAGGCAACGCCCCTTCTGCCG  >  NC_000913/3107069‑3107222                                                                              |                                                                             tCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTg                                                                            >  2:513863/1‑80 (MQ=255)     tataAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCg                                                                       >  1:513890/1‑81 (MQ=255)       taAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGcc                                                                     <  2:513825/81‑1 (MQ=255)            acCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGccccacc                                                                <  1:513827/81‑1 (MQ=255)                  cAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGAc                                                          >  1:513870/1‑81 (MQ=255)                   aGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACt                                                         >  1:513860/1‑81 (MQ=255)                     cTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTc                                                       >  2:513857/1‑81 (MQ=255)                          gCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCgg                                                  <  2:513827/81‑1 (MQ=255)                          gCGAAAATCCCGGCAACAAATACACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCAgg                                                  >  2:513864/1‑81 (MQ=255)                             aaaaTCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTa                                               >  1:513871/1‑81 (MQ=255)                                     ggCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCaa                                       <  1:513829/81‑1 (MQ=255)                                        aacaaATTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCaacg                                     >  2:513875/1‑78 (MQ=255)                                            aaTTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGc                                <  2:513829/81‑1 (MQ=255)                                             aTTCACCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGcc                               <  1:513831/81‑1 (MQ=255)                                                 aCCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAgg                            <  1:513832/80‑1 (MQ=255)                                                 aCCCCTTCTTCCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGt                           >  1:513867/1‑81 (MQ=255)                                                         ttcCAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAgg                   >  1:513891/1‑81 (MQ=255)                                                            cAGTGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAGGCaa                >  1:513886/1‑81 (MQ=255)                                                               tGCAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAGGCAACGc             <  2:513832/81‑1 (MQ=255)                                                                 cAAGGAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAGGCAACGccc           <  2:513833/81‑1 (MQ=255)                                                                    ggAAATAACGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAGGCAACGCCCCtt        <  2:513834/81‑1 (MQ=255)                                                                          aaCGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAGGCAACGCCCCTTCTGCCg  >  2:513871/1‑81 (MQ=255)                                                                          aaCGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAGGCAACGCCCCTTCTGCCg  >  2:513873/1‑81 (MQ=255)                                                                           aCGCTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAGGCAACGCCCCTTCTGCCg  >  1:513878/1‑80 (MQ=255)                                                                              |                                                                             TCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGTTGAACCGCCCCACCCCAGACTTCCCCGGGTACCACGCAAAGCACGCCAGGTGGATAAGGCAACGCCCCTTCTGCCG  >  NC_000913/3107069‑3107222

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑