Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,536,308	G→A	V159M (GTG→ATG)	crr →	Enzyme IIA(Glc)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,536,308	0	G	A	100.0%	168.2 / NA	48	V159M (GTG→ATG)	crr	Enzyme IIA(Glc)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/19);  total (29/19)

GTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAG  >  NC_000913/2536173‑2536443                                                                                                                                        |                                                                                                                                        gTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATggg                                                                                                                                     <  2:424367/140‑1 (MQ=255)  tGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGt                                                                                                                                    <  2:151183/140‑1 (MQ=255)    aaaGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGa                                                                                                                                  <  2:519503/140‑1 (MQ=255)    aaaGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGa                                                                                                                                  <  1:90915/140‑1 (MQ=255)       gTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAAc                                                                                                                               <  2:108075/140‑1 (MQ=255)       gTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAAc                                                                                                                               <  1:61195/140‑1 (MQ=255)       gTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAAc                                                                                                                               <  1:509757/140‑1 (MQ=255)        ttGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAAcc                                                                                                                              <  2:116258/140‑1 (MQ=255)         tGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAAccc                                                                                                                             <  2:855873/140‑1 (MQ=255)          ggCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAAcccc                                                                                                                            >  1:338978/1‑140 (MQ=255)                aCTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTAt                                                                                                                      >  1:1462883/1‑140 (MQ=255)                  tGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATcc                                                                                                                    <  1:1152889/140‑1 (MQ=255)                      aTTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCAt                                                                                                                <  1:1486107/140‑1 (MQ=255)                        tGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCa                                                                                                              <  2:1212727/140‑1 (MQ=255)                        tGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCa                                                                                                              <  1:1008938/140‑1 (MQ=255)                         gAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACT‑‑‑‑AAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCaagaag                                                                                                         >  2:259026/1‑140 (MQ=255)                                    gCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCt                                                                                                  >  2:345366/1‑140 (MQ=255)                                        ctgctgGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGcc                                                                                              >  1:280336/1‑140 (MQ=255)                                           ctgGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCa                                                                                           <  2:269681/140‑1 (MQ=255)                                             ggAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGt                                                                                         >  1:506452/1‑140 (MQ=255)                                                agagAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGaa                                                                                      >  1:985683/1‑140 (MQ=255)                                                   gaAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGaaaaa                                                                                   >  1:118336/1‑140 (MQ=255)                                                     aaGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATg                                                                                 >  2:439216/1‑140 (MQ=255)                                                     aaGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATg                                                                                 >  1:1265820/1‑140 (MQ=255)                                                         cAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGgcgc                                                                             >  2:119877/1‑140 (MQ=255)                                                                  ccTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGgcg                                                                    >  2:347133/1‑140 (MQ=255)                                                                     gACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGTGCCCATCGGCGCCa                                                                 >  1:168389/1‑140 (MQ=255)                                                                      aCTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAt                                                                >  1:1319204/1‑140 (MQ=255)                                                                      aCTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAt                                                                >  2:1323383/1‑140 (MQ=255)                                                                       cTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtt                                                               >  1:1453217/1‑140 (MQ=255)                                                                        tCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACT‑‑‑‑AAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAttttttt                                                          <  1:259026/140‑1 (MQ=255)                                                                           ggTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtttttt                                                           >  1:923084/1‑140 (MQ=255)                                                                            gttgttATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAttttttt                                                          >  1:650418/1‑140 (MQ=255)                                                                                  aTCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCtt                                                    >  1:547072/1‑140 (MQ=255)                                                                                          cATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAgc                                            >  1:123156/1‑140 (MQ=255)                                                                                          cATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAgc                                            >  1:604719/1‑140 (MQ=255)                                                                                             ggACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAgcggc                                         >  1:1180539/1‑140 (MQ=255)                                                                                               aCGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCACATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGcggcgg                                       >  1:77880/1‑140 (MQ=255)                                                                                                   aaTCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCaaa                                   >  1:960787/1‑140 (MQ=255)                                                                                                    aTCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCaaaa                                  >  1:269399/1‑140 (MQ=255)                                                                                                             cTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCAt                         <  2:923084/140‑1 (MQ=255)                                                                                                             cTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCAt                         <  2:118336/140‑1 (MQ=255)                                                                                                             cTGATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCAt                         >  1:323361/1‑140 (MQ=255)                                                                                                               gATCAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCg                       >  2:581058/1‑140 (MQ=255)                                                                                                                  cAAACTGTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGctc                    <  2:338978/140‑1 (MQ=255)                                                                                                                        gTCCGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtgc              >  2:689785/1‑140 (MQ=255)                                                                                                                          ccGGTAGCGTAACCATGGGTGAAACCCCGGTTATCCACATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtgctg            >  2:1287256/1‑140 (MQ=255)                                                                                                                                    aaCCATGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAg  <  1:119877/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        GTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAG  >  NC_000913/2536173‑2536443

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