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%	92.0 / NA	29	K192I (AAA→ATA)	malA	alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (17/12);  total (17/12)

GTAACCAAGTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTTTATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCCATTA  >  CP000730/1624530‑1624813                                                                                                                                                |                                                                                                                                             gTAACCAAGTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATc                                                                                                                                          >  1:276745/1‑148 (MQ=255) gTAACCAAGTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATc                                                                                                                                          >  1:27702/1‑148 (MQ=255) gTAACCAAGTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATc                                                                                                                                          >  2:254832/1‑148 (MQ=255)      cAAGTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAacc                                                                                                                                     <  1:239985/148‑1 (MQ=255)               ataCCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTtca                                                                                                                            <  1:201614/148‑1 (MQ=255)                taCCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTtcat                                                                                                                           <  2:345211/148‑1 (MQ=255)                      acTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTc                                                                                                                     >  2:95950/2‑148 (MQ=255)                       cTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCa                                                                                                                    <  1:158001/148‑1 (MQ=255)                          aTTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATAGAAGCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAAt                                                                                                                 >  1:91945/1‑148 (MQ=255)                                atCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCa                                                                                                           >  1:354670/1‑148 (MQ=255)                                         aaaTGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACt                                                                                                  >  1:365062/1‑148 (MQ=255)                                         aaaTGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACt                                                                                                  >  1:365058/1‑148 (MQ=255)                                            tGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTcc                                                                                               <  1:171045/148‑1 (MQ=255)                                            tGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTcc                                                                                               <  1:152035/148‑1 (MQ=255)                                             gCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCg                                                                                              >  1:313570/1‑148 (MQ=255)                                             gCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCg                                                                                              >  1:359789/1‑148 (MQ=255)                                             gCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCg                                                                                              >  1:200505/1‑148 (MQ=255)                                              cTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCgg                                                                                             <  1:287938/148‑1 (MQ=255)                                                         gTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAAt                                                                                  <  1:118468/148‑1 (MQ=255)                                                                     ggtacaggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTGATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGgtt                                                                      >  1:300338/1‑148 (MQ=255)                                                                        acaggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGgttgtt                                                                   >  2:30002/1‑148 (MQ=255)                                                                                                  cttcttAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAatat                                         >  2:55659/1‑148 (MQ=255)                                                                                                  cttcttAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAatat                                         >  2:123584/1‑148 (MQ=255)                                                                                                            tGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTAt                               <  2:25000/148‑1 (MQ=255)                                                                                                                              ctAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATg             >  1:247610/1‑148 (MQ=255)                                                                                                                                   tCCATCAATGCCTATATCAAACCAACAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGAt        >  2:99706/1‑148 (MQ=255)                                                                                                                                      aTCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCCa     <  1:232917/148‑1 (MQ=255)                                                                                                                                         aaTGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCCATTa  <  2:110680/148‑1 (MQ=255)                                                                                                                                         aaTGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCCATTa  <  2:344999/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                             GTAACCAAGTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTTTATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCCATTA  >  CP000730/1624530‑1624813

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