Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,569,344	(T)6→5	intergenic (‑16/+2)	glgC ← / ← glgX	glucose‑1‑phosphate adenylyltransferase/limit dextrin alpha‑1,6‑glucohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,569,339	0	T	.	100.0%	107.8 / NA	26	intergenic (‑11/+7)	glgC/glgX	glucose‑1‑phosphate adenylyltransferase/limit dextrin alpha‑1,6‑glucohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/14);  total (12/14)

TCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTAT  >  NC_000913/3569259‑3569410                                                                                 |                                                                        tCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCt                                                                         <  2:621636/81‑1 (MQ=255)  cAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCtt                                                                        >  1:621686/1‑81 (MQ=255)   aGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCttt                                                                       <  2:621638/81‑1 (MQ=255)   aGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCttt                                                                       <  1:621640/81‑1 (MQ=255)       cAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCtttttat                                                                   >  1:621691/1‑79 (MQ=255)         aCAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCtttttatca                                                                 <  1:621641/81‑5 (MQ=255)                 cAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTgg                                                         >  2:621681/1‑80 (MQ=255)                         gCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAacacata                                                 >  2:621697/1‑78 (MQ=255)                             gcgcgcCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAAt                                             <  1:621642/80‑1 (MQ=255)                               gcgcCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCg                                          <  1:621643/81‑1 (MQ=255)                               gcgcCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCg                                          <  1:621644/81‑1 (MQ=255)                                 gcCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCgtg                                        >  2:621698/1‑81 (MQ=255)                                      cATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAg                                   >  2:621701/1‑81 (MQ=255)                                      cATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAg                                   <  2:621644/81‑1 (MQ=255)                                        ttAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGt                                 <  1:621646/81‑1 (MQ=255)                                          aaGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTcc                               >  1:621696/1‑81 (MQ=255)                                              gATCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGc                           >  1:621688/1‑81 (MQ=255)                                               aTCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGc                           <  2:621646/80‑1 (MQ=255)                                                tCGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGc                           >  2:621696/1‑79 (MQ=255)                                                 cGTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGc                           >  1:621706/1‑78 (MQ=255)                                                  gTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAg                        <  2:621648/80‑1 (MQ=255)                                                  gTTCTTCTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGa                       <  1:621648/81‑1 (MQ=255)                                                   ttcttcTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGAc                      <  2:621649/81‑1 (MQ=255)                                                   ttcttcTCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGAc                      <  1:621651/81‑1 (MQ=255)                                                      ttctCTAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGc                   >  1:621692/1‑81 (MQ=255)                                                          ctAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGt                >  1:621695/1‑80 (MQ=255)                                                           tAAACTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCg                 <  2:621651/78‑1 (MQ=255)                                                              aCTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAGTCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCa           <  1:621653/81‑1 (MQ=255)                                                               cTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCAc          >  1:621690/1‑81 (MQ=255)                                                               cTAACCATGACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCAc          <  1:621654/81‑1 (MQ=255)                                                                       gACTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTAt  >  2:621703/1‑81 (MQ=255)                                                                        aCTAACTCC‑TTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTAt  <  2:621654/80‑1 (MQ=255)                                                                                 |                                                                        TCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTAT  >  NC_000913/3569259‑3569410

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