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

GATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC  >  NC_000913/3780888‑3781145                                                                                                                                   |                                                                                                                                  gaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTgcg                                                                                                                          <  2:909404/139‑1 (MQ=255)  aTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTgcgc                                                                                                                         <  1:364001/139‑1 (MQ=255)       aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCt                                                                                                                    <  1:829441/139‑1 (MQ=255)            cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCAccc                                                                                                               >  1:109473/1‑139 (MQ=255)                     tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCctct                                                                                                      >  2:383514/1‑139 (MQ=255)                     tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGACTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCctct                                                                                                      <  2:900639/139‑1 (MQ=255)                        cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTccc                                                                                                   <  1:864335/139‑1 (MQ=255)                          ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTc                                                                                                 >  2:918432/1‑139 (MQ=255)                            cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                               <  1:1074105/139‑1 (MQ=255)                                        cAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAcc                                                                                   <  2:109473/139‑1 (MQ=255)                                             aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCg                                                                              <  1:383514/139‑1 (MQ=255)                                             aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCg                                                                              <  1:466075/139‑1 (MQ=255)                                               aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                            <  1:805930/139‑1 (MQ=255)                                                 tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                          <  1:918432/139‑1 (MQ=255)                                                 tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                          <  1:859308/139‑1 (MQ=255)                                                     cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCt                                                                                                                    <  2:554220/93‑1 (MQ=255)                                                     cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCt                                                                                                                    >  1:554220/1‑93 (MQ=255)                                                                            gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCAcc                                                                                                                <  1:62866/74‑1 (MQ=255)                                                                            gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCAcc                                                                                                                >  2:62866/1‑74 (MQ=255)                                                                            gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                           <  1:461285/111‑1 (MQ=255)                                                                            gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                           >  2:461285/1‑111 (MQ=255)                                                                            gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGa                                              >  1:1055317/1‑139 (MQ=255)                                                                                 aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACACCg                                         <  1:945433/139‑1 (MQ=255)                                                                                       tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                           >  1:489434/1‑100 (MQ=255)                                                                                       tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                           <  2:489434/100‑1 (MQ=255)                                                                                       tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACACCGCActct                                   >  2:973215/1‑139 (MQ=255)                                                                                         ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑c                                                        >  2:1063969/1‑116 (MQ=255)                                                                                         ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑c                                                        <  1:1063969/116‑1 (MQ=255)                                                                                                  cTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGAGAACGGTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACAACGCACTCTTTTTTTCTccc                        <  2:539275/139‑1 (MQ=255)                                                                                                   tGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCt                                                                                                                    >  1:145586/1‑47 (MQ=255)                                                                                                   tGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCt                                                                                                                    <  2:145586/47‑1 (MQ=255)                                                                                                     gCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCt                                                                                                                    <  2:996516/45‑1 (MQ=255)                                                                                                     gCTCCCATGGCGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCt                                                                                                                    >  1:996516/1‑45 (MQ=255)                                                                                                         ccATGGCGAAAGGGAATGCGGCATTATTCG‑TTGCCCCCCCTTCACCTAACCCCTCCCCTCAGGGGAAGGGGAACGTTCGGCGCTGTAAGTTCCTCCTCAGCTC‑GAAAAGGCACCGGACCCTTTTTTTCTCCCTCGcccc                 <  2:1060517/139‑1 (MQ=255)                                                                                                               cGAAAGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg           >  1:968345/1‑139 (MQ=255)                                                                                                                 aaaGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                               >  2:644957/1‑54 (MQ=255)                                                                                                                 aaaGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                               <  1:644957/54‑1 (MQ=255)                                                                                                                 aaaGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACACCGCa                                       >  2:844914/1‑109 (MQ=255)                                                                                                                 aaaGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACACCGCa                                       <  1:844914/109‑1 (MQ=255)                                                                                                                   aGGGAATGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga       <  1:973215/139‑1 (MQ=255)                                                                                                                        aTGCGGCA‑TAATCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc  >  2:779299/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                  GATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCA‑TAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCA‑CTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCC  >  NC_000913/3780888‑3781145

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