Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,624,673	T→A	K192I (AAA→ATA)	malA ←	alpha‑glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,624,673	0	T	A	100.0%	145.4 / NA	46	K192I (AAA→ATA)	malA	alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/16);  total (30/16)

GTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTTTATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCCA  >  CP000730/1624538‑1624810                                                                                                                                        |                                                                                                                                          gTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAGTTGCATCTACTCTAAATCCATCAATGCCTATATCa                                                                                                                                      >  1:671767/1‑141 (MQ=255)      tataCCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAacca                                                                                                                                 <  1:737987/141‑1 (MQ=255)               cTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTCATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTtcat                                                                                                                        >  1:362723/1‑141 (MQ=255)                     cATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATGTCAAACCACCAGTTCATCATTTc                                                                                                                  >  1:554025/1‑141 (MQ=255)                               tCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCa                                                                                                        >  1:261707/1‑141 (MQ=255)                               acAAATGCTGAGGCATATATTTAAACCTCAGGTACAGGTAAGGCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTCAATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCa                                                                                                        <  1:198607/140‑1 (MQ=255)                                      cTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTaa                                                                                                 <  1:87940/141‑1 (MQ=255)                                               atGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCACCATTTCAAATACAGCATCTCTAACTTCCGGAt                                                                                        >  1:1013361/1‑141 (MQ=255)                                                   tttACCCTCCGGTACAGGTAAGTCCCCCGCTTCAACCGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTAcc                                                                                    <  1:1141278/141‑1 (MQ=255)                                                   tttACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTAcc                                                                                    <  1:450831/141‑1 (MQ=255)                                                   tttACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTAcc                                                                                    <  1:849537/141‑1 (MQ=255)                                                       cccTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCaa                                                                                >  1:130088/1‑141 (MQ=255)                                                       cccTCAGGTACAGGTAAGTCACCCGCTGCAAACGTCTTCTTAATGTGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCaa                                                                                >  1:641158/1‑141 (MQ=255)                                                        ccTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAGATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAAt                                                                               >  1:164194/1‑141 (MQ=255)                                                        ccTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAAt                                                                               >  1:2934/1‑141 (MQ=255)                                                        ccTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCAGCATTTCCAATACAGCATCTATAACTTCCGGATTACCCCAAt                                                                               >  1:29791/1‑141 (MQ=255)                                                           caggtacaggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCCTATAACCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCa                                                                            >  1:35148/1‑141 (MQ=255)                                                           caggtacaggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCCCCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCa                                                                            >  1:35444/1‑141 (MQ=255)                                                                 caggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTGAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAg                                                                      >  1:476291/1‑141 (MQ=255)                                                                 caggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAg                                                                      >  1:50743/1‑141 (MQ=255)                                                                 caggtaAGTCACCCGCGTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCGCTGACTTCCGGATTCACCCAATTCAACTCAg                                                                      >  1:189814/1‑141 (MQ=255)                                                                   ggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGgt                                                                    >  1:493463/1‑141 (MQ=255)                                                                   ggtaAGTCAACCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGgt                                                                    <  1:448633/141‑1 (MQ=255)                                                                        gTCACCCGCTTCAAACTTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCGATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGttt                                                               <  1:1179846/141‑1 (MQ=255)                                                                           aCCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTCCTCTAACTCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGtttttt                                                            >  1:230291/1‑141 (MQ=255)                                                                             ccGCTTCAAACGTCTTCTTAATATACGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTAc                                                          >  1:870017/1‑141 (MQ=255)                                                                                    aaaCGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATaa                                                   >  1:1173457/1‑141 (MQ=255)                                                                                      aCGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAAt                                                 >  1:437271/1‑141 (MQ=255)                                                                                             cttAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAaata                                          >  1:419746/1‑141 (MQ=255)                                                                                                aaTATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAataata                                       <  1:679046/141‑1 (MQ=255)                                                                                                 atatGCGTAACTGCATCTACACTAAATCCACTATTGCCTATATCAAACCACCAGGTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAatat                                      <  1:103754/141‑1 (MQ=255)                                                                                                       gTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCa                                <  1:539203/141‑1 (MQ=255)                                                                                                          aTTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCAATTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCCGTa                             <  1:858248/141‑1 (MQ=255)                                                                                                             gCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTa                          >  1:790009/1‑141 (MQ=255)                                                                                                             gCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTa                          >  1:894976/1‑141 (MQ=255)                                                                                                              cATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAg                         >  1:1051280/1‑141 (MQ=255)                                                                                                              cATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAg                         >  1:379999/1‑141 (MQ=255)                                                                                                                tCTACTCTAAATCCATCAATGCCTATATCAAAGCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATAAGGTTGTTTTTTACTGAATAAATGGAAGTAATATTACGAAGTATTTGca                       >  1:528004/1‑141 (MQ=255)                                                                                                                tCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGca                       <  1:496980/141‑1 (MQ=255)                                                                                                                tCTACCCTAACTCCATCAATTCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGca                       <  1:43457/141‑1 (MQ=255)                                                                                                                      ctAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCatat                 <  1:984662/141‑1 (MQ=255)                                                                                                                      ctAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCatat                 >  1:1091707/1‑141 (MQ=255)                                                                                                                      ctAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTACTCAGTATTAGCATCatat                 >  1:900937/1‑141 (MQ=255)                                                                                                                       tAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACCGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCGGTATTAGCATCATAtt                >  1:702714/1‑141 (MQ=255)                                                                                                                           tCCATCAATGCCTATATCACACCACCAGTTCATCATTTCAAATACAGCATCTCTCACTTCCGGCTTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATACATGGAAATAATATTGCTCGGTCTTAGCATCATTTTTCCa            >  1:356692/1‑141 (MQ=255)                                                                                                                                     ccTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCCa  <  1:527884/141‑1 (MQ=255)                                                                                                                                        |                                                                                                                                          GTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTTTATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCCA  >  CP000730/1624538‑1624810

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