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

AAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC  >  NC_000913/3780881‑3781145                                                                                                                                         |                                                                                                                                 aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATc                                                                                                                                >  2:494351/1‑139 (MQ=255)            aaaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGcccc                                                                                                                     >  1:718731/1‑139 (MQ=255)              aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAATAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGTGCCCCCt                                                                                                                   <  1:962227/139‑1 (MQ=255)              aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  2:691503/139‑1 (MQ=255)              aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  1:370635/139‑1 (MQ=255)              aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  1:360788/139‑1 (MQ=255)              aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  1:817986/139‑1 (MQ=255)              aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  1:980109/139‑1 (MQ=255)              aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  1:991713/139‑1 (MQ=255)              aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  2:1313012/139‑1 (MQ=255)                   cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCAccc                                                                                                              <  2:1189674/139‑1 (MQ=255)                         cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGccccc                                                                                                                    >  1:31000/1‑127 (MQ=255)                         cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGccccc                                                                                                                    <  2:31000/127‑1 (MQ=255)                          gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCct                                                                                                       <  1:499721/139‑1 (MQ=255)                               cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTccc                                                                                                  <  2:1258279/139‑1 (MQ=255)                                        aaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                         >  1:551136/1‑139 (MQ=255)                                                    aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCg                                                                             <  2:995218/139‑1 (MQ=255)                                                       ttACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                          <  1:812913/139‑1 (MQ=255)                                                       ttACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                          <  2:718731/139‑1 (MQ=255)                                                       ttACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                          <  2:760779/139‑1 (MQ=255)                                                       gtACGCGAGTTTCTCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGTCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                          <  2:395741/138‑1 (MQ=39)                                                        tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                         <  2:551136/139‑1 (MQ=255)                                                            cAAGATTCGCTGGTGCAGGGGCT‑GGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                    >  1:681199/1‑139 (MQ=255)                                                               gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCctct                                                                                                     <  2:1349087/104‑1 (MQ=255)                                                               gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCctct                                                                                                     >  1:1349087/1‑104 (MQ=255)                                                                 ttCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAc                                                                >  2:555201/1‑139 (MQ=255)                                                                         tGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAc                                                        <  1:176003/139‑1 (MQ=255)                                                                           cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGGATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActc                                                      <  1:790642/139‑1 (MQ=255)                                                                              gggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                          <  1:1024397/116‑1 (MQ=255)                                                                              gggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                          >  2:1024397/1‑116 (MQ=255)                                                                                gCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                    <  2:94814/104‑1 (MQ=255)                                                                                gCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                    >  1:94814/1‑104 (MQ=255)                                                                                   gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGtt                                                                                                                             >  1:1065029/1‑60 (MQ=255)                                                                                   gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGtt                                                                                                                             <  2:1065029/60‑1 (MQ=255)                                                                                               gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                     >  2:1057084/1‑104 (MQ=255)                                                                                               gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                     <  1:1057084/104‑1 (MQ=255)                                                                                                  tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   >  2:1020809/1‑55 (MQ=255)                                                                                                  tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCt                                                                                                                   <  1:1020809/55‑1 (MQ=255)                                                                                                  tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAAccc                                                                                                        >  1:200816/1‑66 (MQ=255)                                                                                                  tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAAccc                                                                                                        <  2:200816/66‑1 (MQ=255)                                                                                                  tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                               >  2:1034168/1‑139 (MQ=255)                                                                                                    cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGccccc                                                                                                                    >  2:1082793/1‑52 (MQ=255)                                                                                                    cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAATCGTTTGCGccccc                                                                                                                    <  1:1082793/52‑1 (MQ=255)                                                                                                                 cATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                            >  2:1080683/1‑79 (MQ=255)                                                                                                                 cATGGCGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                            <  1:1080683/79‑1 (MQ=255)                                                                                                                      cGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                         <  2:489909/77‑1 (MQ=255)                                                                                                                      cGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                         >  1:489909/1‑77 (MQ=255)                                                                                                                      cGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCa                                       <  2:149090/111‑1 (MQ=255)                                                                                                                      cGAAAGGGAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCa                                       >  1:149090/1‑111 (MQ=255)                                                                                                                             gAATGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg    <  2:681199/139‑1 (MQ=255)                                                                                                                               aTGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTccc                                                                                                  >  1:797452/1‑43 (MQ=255)                                                                                                                               aTGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTccc                                                                                                  <  2:797452/43‑1 (MQ=255)                                                                                                                               aTGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc  <  1:555201/139‑1 (MQ=255)                                                                                                                                tGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                          >  1:1140090/1‑66 (MQ=255)                                                                                                                                tGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                          <  2:1140090/66‑1 (MQ=255)                                                                                                                                tGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                        >  2:731020/1‑68 (MQ=255)                                                                                                                                tGCGGCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                        <  1:731020/68‑1 (MQ=255)                                                                                                                                   ggCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                            <  2:1206364/109‑1 (MQ=255)                                                                                                                                   ggCATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                            >  1:1206364/1‑109 (MQ=255)                                                                                                                                     cATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGtt                                                                               >  2:511541/1‑56 (MQ=255)                                                                                                                                     cATAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGtt                                                                               <  1:511541/56‑1 (MQ=255)                                                                                                                                         |                                                                                                                                 AAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC  >  NC_000913/3780881‑3781145

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