Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,204,617	G→A	intergenic (‑28/‑77)	folB ← / → plsY	dihydroneopterin aldolase/putative glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,204,617	0	G	A	100.0%	113.5 / NA	34	intergenic (‑28/‑77)	folB/plsY	dihydroneopterin aldolase/putative glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/16);  total (18/16)

GCGCCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCGGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTATCGGTGGCAAGGGAGCAGCCGTAGC  >  NC_000913/3204388‑3204863                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       gcgcCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAAtg                                                                                                                                                                                                                                   <  2:159664/251‑1 (MQ=255)  cgcCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAAtgt                                                                                                                                                                                                                                  <  1:8184/251‑1 (MQ=255)    cccccTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGc                                                                                                                                                                                                                                <  1:259515/251‑1 (MQ=255)               tGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCcac                                                                                                                                                                                                                     >  1:294350/1‑251 (MQ=255)                        accGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATc                                                                                                                                                                                                            >  2:68848/1‑251 (MQ=255)                               cTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAac                                                                                                                                                                                                     >  1:214836/1‑251 (MQ=255)                                          gCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCa                                                                                                                                                                                          <  2:260206/251‑1 (MQ=255)                                           cGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCa                                                                                                                                                                                          >  2:276044/1‑250 (MQ=255)                                              accTGAGGAAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAAc                                                                                                                                                                                      <  2:34247/249‑1 (MQ=255)                                                      aaaTCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAAGTTTCAAAACGGAACAGc                                                                                                                                                                              <  1:273067/250‑1 (MQ=255)                                                                                         gTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGTAATGATCCTCATc                                                                                                                                           <  1:51230/251‑1 (MQ=255)                                                                                                      aTTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGATTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCg                                                                                                                              <  1:195135/251‑1 (MQ=255)                                                                                                            atatCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTTGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCa                                                                                                                        <  2:118258/251‑1 (MQ=255)                                                                                                              atCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCAt                                                                                                                       <  1:126314/250‑1 (MQ=255)                                                                                                                cGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTc                                                                                                                    >  1:293380/1‑251 (MQ=255)                                                                                                                   acacTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTcc                                                                                                                   >  1:288920/1‑249 (MQ=255)                                                                                                                     acTAACGTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTg                                                                                                               <  2:288920/251‑1 (MQ=255)                                                                                                                                 tCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGttt                                                                                                   >  2:104679/1‑251 (MQ=255)                                                                                                                                    aTGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGcc                                                                                                <  1:221163/251‑1 (MQ=255)                                                                                                                                           gTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTtgtg                                                                                         >  2:220594/1‑251 (MQ=255)                                                                                                                                              cccAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTgg                                                                                      <  1:272173/251‑1 (MQ=255)                                                                                                                                                       aaaCACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCgg                                                                             >  2:28328/1‑250 (MQ=255)                                                                                                                                                                   tggtgATTACCGAAAGTTGCTCTATAAATAAAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGCTTCATCGTTGACAACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAAcc                                                                   >  1:276026/1‑249 (MQ=255)                                                                                                                                                                     gtgATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCg                                                               >  2:27681/1‑251 (MQ=255)                                                                                                                                                                        aTTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCt                                                            >  1:70289/1‑251 (MQ=255)                                                                                                                                                                                                    tatCCATCACGTCCTGCATGCTTTTTGGCTAACCAGATACCACTTCAGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACAAATGTGtt                                >  1:281597/1‑251 (MQ=255)                                                                                                                                                                                                       ccATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTTCCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACg                             <  1:198514/251‑1 (MQ=255)                                                                                                                                                                                                       ccATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACg                             >  2:45541/1‑251 (MQ=255)                                                                                                                                                                                                           cACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTATc                         <  1:276044/251‑1 (MQ=255)                                                                                                                                                                                                              gTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTATCGGt                      >  2:305908/1‑251 (MQ=255)                                                                                                                                                                                                                         tttttGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCCATGTGTTACGTATCGGTGGCAAGGGagc           >  2:77006/1‑251 (MQ=255)                                                                                                                                                                                                                            ttgacTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTATCGGTGGCAAGGGagcagc        >  1:213882/5‑251 (MQ=255)                                                                                                                                                                                                                                 tAACCAGATACCACTTCCGACGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATGCTGGTTTGCCGCGTGTGAGGGCTGCCCGATCCGCGAACCAGCAGCTCCGGCAATCCAGGGGCAACCAGTGTGTTACGTATCGGTGGCAAGGCAGCAGCCGt     >  2:61577/1‑249 (MQ=255)                                                                                                                                                                                                                                  aaCCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTATCGGTGGCAAGGGAGCAGCCGTAGc  <  2:289860/251‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       GCGCCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCGGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTATCGGTGGCAAGGGAGCAGCCGTAGC  >  NC_000913/3204388‑3204863

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