Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,565,112	G→A	G490G (GGC→GGT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,565,112	0	G	A	100.0%	143.5 / NA	40	G490G (GGC→GGT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/20);  total (20/20)

GCGTGATTAACCATTGGGAAATCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGGCCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCACACCGTTGGT  >  NC_000913/3565038‑3565189                                                                           |                                                                              gCGTGATTAACCATTGGGAAATCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAgg                                                                         >  2:620536/1‑81 (MQ=255) gCGTGATTAACCATTGGGAAATCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAgg                                                                         >  1:620536/1‑81 (MQ=255)          aCCATTGGGAAATCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTcc                                                                <  2:620483/81‑1 (MQ=255)          aCCATTGGGAAATCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTc                                                                 >  2:620558/1‑80 (MQ=255)            cATTGGGAAATCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAg                                                              >  1:620549/1‑81 (MQ=255)                gggAAATCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCAc                                                           <  1:620485/80‑1 (MQ=255)                   aaaTCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCt                                                       >  1:620550/1‑81 (MQ=255)                   aaaTCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCt                                                       >  1:620546/1‑81 (MQ=255)                              ttgttgCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCTGGTGTTCGTCCAGCACGGCTGAAAGCGATgg                                            >  2:620542/1‑81 (MQ=255)                              ttgttgCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATgg                                            >  2:620539/1‑81 (MQ=255)                                 ttgCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGtt                                         <  1:620486/81‑1 (MQ=255)                                   gCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTc                                        <  2:620486/80‑1 (MQ=255)                                       cTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCTATGGGTTCGCTAc                                   <  1:620489/81‑1 (MQ=255)                                       cTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTAc                                   <  2:620487/81‑1 (MQ=255)                                          aTTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCgc                                >  1:620570/1‑81 (MQ=255)                                             aaGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCgcca                              <  2:620489/80‑1 (MQ=255)                                              aGCAGGCTAAGTTCGGTGCGCCAGTTACTACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCgccag                             <  1:620491/80‑1 (MQ=255)                                              aGCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCgccagc                            >  1:620571/1‑81 (MQ=255)                                               gCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCATCACGGCTGAAAGCGATGGGTTCGCTACCgccagc                            >  1:620562/1‑80 (MQ=255)                                               gCAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCgccagcc                           <  2:620491/81‑1 (MQ=255)                                                cAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCgccagcca                          <  1:620494/81‑1 (MQ=255)                                                cAGGCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCgccagcca                          <  2:620492/81‑1 (MQ=255)                                                  ggCTAAGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCcagccagc                        >  2:620556/1‑81 (MQ=255)                                                       aGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGAcgc                   >  2:620575/1‑81 (MQ=255)                                                       aGGTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGAc                     <  1:620495/79‑1 (MQ=255)                                                         gTCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACgcgg                 <  1:620496/81‑1 (MQ=255)                                                          tCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACgcggc                <  2:620497/81‑1 (MQ=255)                                                          tCGGTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACgcggc                <  1:620497/81‑1 (MQ=255)                                                            ggTGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGt              >  1:620553/1‑81 (MQ=255)                                                              tGCGCCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTca            >  2:620554/1‑81 (MQ=255)                                                                cgcCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTcac           <  2:620499/80‑1 (MQ=255)                                                                cgcCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTcac           <  1:620499/80‑1 (MQ=255)                                                                 gcCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCacac         >  2:620553/1‑81 (MQ=255)                                                                 gcCAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCacac         >  2:620571/1‑81 (MQ=255)                                                                  ccAGTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCACAcc        <  2:620500/81‑1 (MQ=255)                                                                     gTTACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCACACcgtt     <  1:620502/81‑1 (MQ=255)                                                                      ttACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCACACcgttg    >  1:620566/1‑81 (MQ=255)                                                                      ttACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCACACcgttg    <  1:620503/81‑1 (MQ=255)                                                                      ttACGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCACACcgttg    >  2:620580/1‑81 (MQ=255)                                                                        aCGACCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCACACcgttggt  >  2:620560/1‑81 (MQ=255)                                                                           |                                                                              GCGTGATTAACCATTGGGAAATCACAGTGTTGTTGCAGCTCATTAAGCAGGCTAAGGTCGGTGCGCCAGTTACGGCCCAGGTGTTCGTCCAGCACGGCTGAAAGCGATGGGTTCGCTACCGCCAGCCAGCGACGCGGCGTCACACCGTTGGT  >  NC_000913/3565038‑3565189

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