Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	798,688	C→T	P35S (CCG→TCG)	pgl →	6‑phosphogluconolactonase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	798,688	0	C	T	100.0%	148.1 / NA	38	P35S (CCG→TCG)	pgl	6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/19);  total (19/19)

TAATTTTAGCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGCCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCTGCC  >  NC_000913/798541‑798827                                                                                                                                                    |                                                                                                                                            tAATTTTAGCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCgg                                                                                                                                          >  2:591178/1‑151 (MQ=255)      ttAGCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCa                                                                                                                                      <  1:706727/150‑1 (MQ=255)        aGCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGt                                                                                                                                   <  2:178776/151‑1 (MQ=255)         gCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTAGGGGCAGGTg                                                                                                                                  >  1:674614/1‑151 (MQ=255)               aCAGTCAGTTGCTAATTGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTgaatt                                                                                                                              >  1:361204/1‑145 (MQ=255)                cagtcagtTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCg                                                                                                                            >  2:567653/1‑150 (MQ=255)                cagtcagtTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTCGTCGATGTGTCGGGGCAGGTGCAGCCg                                                                                                                            >  1:230305/1‑150 (MQ=255)                    cagtTGCTAAATGCAAAGGAGAATTCATGAAGCAAACAGTTTATAGCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGAtg                                                                                                                         <  2:305603/149‑1 (MQ=255)                     agtTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGAtggtg                                                                                                                      >  2:536398/1‑151 (MQ=255)                         gCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCa                                                                                                                  <  2:89839/151‑1 (MQ=255)                            aaaTGCAAAGGAGCATTTATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGGTGGTGGTCAg                                                                                                                 >  1:90242/1‑149 (MQ=255)                            aaaTGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGcc                                                                                                               >  1:554735/1‑151 (MQ=255)                               tGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCgg                                                                                                            >  2:652606/1‑151 (MQ=255)                                   aaGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACaa                                                                                                        <  1:370251/151‑1 (MQ=255)                                    aGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACaaa                                                                                                       <  2:122833/151‑1 (MQ=255)                                      gAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACg                                                                                                     <  2:475766/151‑1 (MQ=255)                                        gCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGAACAAACGtt                                                                                                   <  1:633999/151‑1 (MQ=255)                                          aTTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTAt                                                                                                 <  1:405753/151‑1 (MQ=255)                                                 aaGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTAt                                                                                           >  2:267171/1‑150 (MQ=255)                                                     aaaCAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTgg                                                                                      >  2:129644/1‑151 (MQ=255)                                                        cAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGgtgt                                                                                   <  2:212605/151‑1 (MQ=255)                                                                atCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTg                                                                           <  2:192103/151‑1 (MQ=255)                                                                    ccAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGtt                                                                       >  2:36550/1‑151 (MQ=255)                                                                      aGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTc                                                                     >  1:140412/1‑151 (MQ=255)                                                                       gCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTcg                                                                    >  2:713822/1‑151 (MQ=255)                                                                       gCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTcg                                                                    <  1:743217/151‑1 (MQ=255)                                                                             agagCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCt                                                              >  2:326981/1‑151 (MQ=255)                                                                               agCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTgg                                                            <  2:128395/151‑1 (MQ=255)                                                                                gCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTgg                                                            >  1:77860/1‑150 (MQ=255)                                                                                        aTTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGcgtatcgt                                                    >  2:730591/1‑150 (MQ=255)                                                                                                                gAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTg                            <  1:590180/150‑1 (MQ=255)                                                                                                                gAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTg                            <  2:712159/150‑1 (MQ=255)                                                                                                                  aGGCGCACTGACGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGAcc                         <  1:669198/151‑1 (MQ=255)                                                                                                                             cGCTGACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCaga              <  2:186189/151‑1 (MQ=255)                                                                                                                                 gACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTc           <  1:373339/150‑1 (MQ=255)                                                                                                                                 gACACAGGTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCt          >  1:140998/1‑151 (MQ=255)                                                                                                                                        gTTGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCTGc   >  2:487756/1‑151 (MQ=255)                                                                                                                                          tGTCGATGTGTCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCTGcc  <  2:683919/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                            TAATTTTAGCTGTTACAGTCAGTTGCTAAATGCAAAGGAGCATTCATGAAGCAAACAGTTTATATCGCCAGCCCTGAGAGCCAGCAAATTCACGTCTGGAATCTGAATCATGAAGGCGCACTGACGCTGACACAGGTTGTCGATGTGCCGGGGCAGGTGCAGCCGATGGTGGTCAGCCCGGACAAACGTTATCTCTATGTTGGTGTTCGCCCTGAGTTTCGCGTCCTGGCGTATCGTATCGCCCCGGACGATGGCGCACTGACCTTTGCCGCAGAGTCTGCGCTGCC  >  NC_000913/798541‑798827

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