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%	124.8 / NA	39	*397* (TAG→TAA)	lldD	L‑lactate dehydrogenase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/24);  total (15/24)

AAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC  >  NC_000913/3780884‑3781145                                                                                                                                      |                                                                                                                                 aaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGtt                                                                                                                             <  1:1249042/139‑1 (MQ=255)     gaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTgcg                                                                                                                         <  1:135713/139‑1 (MQ=255)     gaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTgcg                                                                                                                         <  2:1152586/139‑1 (MQ=255)           aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  2:1290236/139‑1 (MQ=255)           aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  1:1309656/139‑1 (MQ=255)           aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  1:431589/139‑1 (MQ=255)             tGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCa                                                                                                                 >  1:643110/1‑139 (MQ=255)                 gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCt                                                                                                             <  2:643110/139‑1 (MQ=255)                    gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAAc                                                                                                          >  2:1061842/1‑139 (MQ=255)                         tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCctct                                                                                                     <  2:66053/139‑1 (MQ=255)                               gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCAGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCAGCCTAACCCTCTCCCTCa                                                                                               <  2:878248/139‑1 (MQ=255)                                cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   >  2:1128351/1‑118 (MQ=255)                                cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  1:1128351/118‑1 (MQ=255)                                          aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                    <  2:1149788/139‑1 (MQ=255)                                                     tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                         <  1:355888/139‑1 (MQ=255)                                                         cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                     <  1:609995/139‑1 (MQ=255)                                                            gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCa                                                                                                                 <  1:1288389/92‑1 (MQ=255)                                                            gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCa                                                                                                                 >  2:1288389/1‑92 (MQ=255)                                                                             gCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                            <  2:657276/128‑1 (MQ=255)                                                                             gCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                            >  1:657276/1‑128 (MQ=255)                                                                               tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  2:319118/71‑1 (MQ=255)                                                                               tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   >  1:319118/1‑71 (MQ=255)                                                                               tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                         <  1:959584/113‑1 (MQ=255)                                                                               tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                         >  2:959584/1‑113 (MQ=255)                                                                                gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                           >  2:308243/1‑110 (MQ=255)                                                                                gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                           <  1:308243/110‑1 (MQ=255)                                                                                 ggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGtt                                                                               >  2:29647/1‑105 (MQ=255)                                                                                 ggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGtt                                                                               <  1:29647/105‑1 (MQ=255)                                                                                           tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActct                                   >  1:854194/1‑139 (MQ=255)                                                                                             ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                 <  2:1077969/107‑1 (MQ=255)                                                                                             ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                 >  1:1077969/1‑107 (MQ=255)                                                                                                gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCct                                                                                                       >  1:807122/1‑66 (MQ=255)                                                                                                gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCct                                                                                                       <  2:807122/66‑1 (MQ=255)                                                                                                  ggCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                         >  2:57998/1‑94 (MQ=255)                                                                                                  ggCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                         <  1:57998/94‑1 (MQ=255)                                                                                                              cATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCt                <  1:1061842/139‑1 (MQ=255)                                                                                                                           aaTGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                       <  1:928148/71‑1 (MQ=255)                                                                                                                           aaTGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                       >  2:928148/1‑71 (MQ=255)                                                                                                                            aTGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc  >  2:1498356/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                 AAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC  >  NC_000913/3780884‑3781145

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