Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,781,017	G→A	*397* (TAG→TAA)	lldD →	L‑lactate dehydrogenase, FMN‑linked

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,017	0	G	A	100.0%	91.8 / NA	30	*397* (TAG→TAA)	lldD	L‑lactate dehydrogenase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/21);  total (9/21)

GATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCG  >  NC_000913/3780888‑3781126                                                                                                                                  |                                                                                                              gaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTgcg                                                                                                      <  1:545110/139‑1 (MQ=255) gaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTgcg                                                                                                      <  2:413249/139‑1 (MQ=255)       aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                <  2:95137/139‑1 (MQ=255)       aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                <  2:696834/139‑1 (MQ=255)       aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                >  1:668902/1‑139 (MQ=255)            cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCAccc                                                                                           <  1:428702/139‑1 (MQ=255)               tGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTaa                                                                                        >  2:641095/1‑139 (MQ=255)                   gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCct                                                                                    <  1:661462/139‑1 (MQ=255)                   gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCct                                                                                    <  2:146686/139‑1 (MQ=255)                     tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCctct                                                                                  <  2:366925/139‑1 (MQ=255)                                                 tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                      <  2:324654/139‑1 (MQ=255)                                                  aCGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGTCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGTGGACCGTTCGGCGCt                                                     <  1:535833/139‑1 (MQ=255)                                                  aCGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                     <  2:19657/139‑1 (MQ=255)                                                       aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatg                                                <  2:60187/139‑1 (MQ=255)                                                        gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCct                                                                                    >  2:100550/1‑102 (MQ=255)                                                        gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCct                                                                                    <  1:100550/102‑1 (MQ=255)                                                                  tGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCctct                                                                                  <  2:374004/94‑1 (MQ=255)                                                                  tGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCctct                                                                                  >  1:374004/1‑94 (MQ=255)                                                                    cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActc                                   <  2:483280/139‑1 (MQ=255)                                                                       gggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                >  1:515861/1‑139 (MQ=255)                                                                                       tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                         >  1:667121/1‑130 (MQ=255)                                                                                       tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                         <  2:667121/130‑1 (MQ=255)                                                                                           tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGtt                                                            <  2:778981/91‑1 (MQ=255)                                                                                           tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGtt                                                            >  1:778981/1‑91 (MQ=255)                                                                                                     gCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg  <  2:515861/139‑1 (MQ=255)                                                                                                     gCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg  <  2:668902/139‑1 (MQ=255)                                                                                                                        aTGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                      >  1:800358/1‑52 (MQ=255)                                                                                                                        aTGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                      <  2:800358/52‑1 (MQ=255)                                                                                                                              cATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                      >  2:143660/1‑62 (MQ=255)                                                                                                                              cATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                      <  1:143660/62‑1 (MQ=255)                                                                                                                                  |                                                                                                              GATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCG  >  NC_000913/3780888‑3781126

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