Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,971	G→A	A22V (GCT→GTT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,971	0	G	A	100.0%	130.9 / NA	33	A22V (GCT→GTT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/21);  total (12/21)

AACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCG  >  REL606/3497899‑3498042                                                                         |                                                                        aaCGTTGCGATCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACg                                                                       >  2:547417/1‑75 (MQ=255)      tGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATaga                                                                  >  2:217385/1‑75 (MQ=255)                     aTGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTc                                                   <  2:1419304/75‑1 (MQ=255)                      tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCt                                                  <  1:757604/75‑1 (MQ=255)                      tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCt                                                  <  2:1229491/75‑1 (MQ=255)                      tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCt                                                  >  2:804397/1‑75 (MQ=255)                       gTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTa                                                 >  2:1098344/1‑75 (MQ=255)                        tttATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTAc                                                <  1:355253/75‑1 (MQ=255)                        tttATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTAc                                                <  2:1548471/75‑1 (MQ=255)                        tttATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTAc                                                <  1:141320/75‑1 (MQ=255)                                  acgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTg                                      <  1:217385/75‑1 (MQ=255)                                    gacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTggg                                    <  1:30391/75‑1 (MQ=255)                                    gacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTggg                                    <  1:913962/75‑1 (MQ=255)                                     acCGGGTCCTTTCCAATCGTAAACATCACCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGc                                   >  1:563488/1‑75 (MQ=255)                                        gggTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGAt                                >  1:18231/1‑75 (MQ=255)                                        gggTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGAt                                <  2:6830/75‑1 (MQ=255)                                         ggTCCTTTCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATg                               <  1:810939/75‑1 (MQ=255)                                              tttCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAatat                          >  1:1269487/1‑75 (MQ=255)                                                tCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGt                        <  2:836207/75‑1 (MQ=255)                                                tCCAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGt                        <  2:795107/75‑1 (MQ=255)                                                tCCAAGCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGt                        <  1:1271292/75‑1 (MQ=255)                                                  cAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTaa                      <  2:1563023/75‑1 (MQ=255)                                                  cAATCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTaa                      <  1:479480/75‑1 (MQ=255)                                                    aTCGTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAAc                    <  2:827417/75‑1 (MQ=255)                                                       gTAAACATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGa                 >  2:540228/1‑75 (MQ=255)                                                         aaaCATCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGc               >  2:785375/1‑75 (MQ=255)                                                              tCAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCa          <  2:171825/75‑1 (MQ=255)                                                               cAGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCat         >  1:476831/1‑75 (MQ=255)                                                                aGCTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCata        <  2:202745/75‑1 (MQ=255)                                                                  cTTGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCAtata      >  2:1593093/1‑75 (MQ=255)                                                                   ttGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAg     <  2:378081/75‑1 (MQ=255)                                                                    tGTAAACGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAgg    <  2:286904/75‑1 (MQ=255)                                                                      tAAACGATAGAGTGCTTAAGGGATTCTACGCTAAGCGAGGGCGATGAATATGTAAACGGAGCATTCATATAGGCg  >  2:957517/1‑75 (MQ=255)                                                                         |                                                                        AACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCG  >  REL606/3497899‑3498042

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