Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,926,292	C→T	intergenic (+96/‑16)	queF → / → ppnN	7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,926,292	0	C	T	100.0%	111.1 / NA	34	intergenic (+96/‑16)	queF/ppnN	7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/15);  total (19/15)

TGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTG  >  NC_000913/2926045‑2926529                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              tGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCAt                                                                                                                                                                                                                                            >  2:77084/1‑251 (MQ=255)          tCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGttt                                                                                                                                                                                                                                   >  2:328329/1‑251 (MQ=255)               aaTGACCTGTTACGCTTCTGCCATCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTTGTTCTGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCATGTCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCtt                                                                                                                                                                                                                              <  1:367281/251‑1 (MQ=255)                            gCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGAGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATAtt                                                                                                                                                                                                                 >  2:76039/1‑251 (MQ=255)                                cTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGAGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGGCTGGTTCGGCAATAAATTTTTTCTAAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACAAATATTAGcc                                                                                                                                                                                                             >  1:62303/1‑251 (MQ=255)                                            aaaaTTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCa                                                                                                                                                                                                 <  1:47056/251‑1 (MQ=255)                                              aaTTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATg                                                                                                                                                                                               <  1:267135/251‑1 (MQ=255)                                                     cGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtgt                                                                                                                                                                                        >  1:204912/1‑251 (MQ=255)                                                                      cccGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTg                                                                                                                                                                       <  1:147866/251‑1 (MQ=255)                                                                             tGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGc                                                                                                                                                                >  1:327419/1‑251 (MQ=255)                                                                                cGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCtt                                                                                                                                                              >  1:216119/1‑250 (MQ=255)                                                                                    cTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACg                                                                                                                                                         >  2:14554/1‑251 (MQ=255)                                                                                         cATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCg                                                                                                                                                    <  1:166539/251‑1 (MQ=255)                                                                                                ccGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGGGGATATGCTTAAACGCACCGCcagcag                                                                                                                                             <  2:304298/251‑1 (MQ=255)                                                                                                         aGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATACATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACctct                                                                                                                                     >  1:240329/1‑250 (MQ=255)                                                                                                                      ttGTCCCATCGACCACAGGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGAATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACAGATGAGCCCGCTTGGCTCCATGGATATGTTGCCGCAGCGGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCg                                                                                                                        <  2:346951/250‑1 (MQ=255)                                                                                                                                 aCCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACt                                                                                                            >  2:24989/1‑251 (MQ=255)                                                                                                                                        gACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGTTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGt                                                                                                      >  2:87441/1‑250 (MQ=255)                                                                                                                                         aCTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTAc                                                                                                    <  1:112940/251‑1 (MQ=255)                                                                                                                                                              ttCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGAcc                                                                                <  1:62915/250‑1 (MQ=255)                                                                                                                                                              ttCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGAcc                                                                                <  1:317690/250‑1 (MQ=255)                                                                                                                                                                       ttGCGTGCTGGACACCAGCAGCAGGTTGTGAACAGTCATCCGGCAGGGCTATTGTAATCCAAGGGAATGACGATATTCGACTCATAAGGCGTTTTTCTTGCTTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATCTGCTTAAACGCACCGCCAGCAGCGACCTCTCTCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCAGGTAGTTTGACCGATAACAGc                                                                       >  2:25275/1‑250 (MQ=255)                                                                                                                                                                         gCGTGCTGGATTCACGCAGAAGGTTGTAAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATCAGGAGTTTTTCTTGATTACACATATTAGCACGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCaa                                                                     >  1:288913/1‑250 (MQ=255)                                                                                                                                                                            tGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACt                                                                                                            <  1:19860/208‑1 (MQ=255)                                                                                                                                                                            tGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACt                                                                                                            >  2:19860/1‑208 (MQ=255)                                                                                                                                                                             gCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGATTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGa                                                                  <  2:299335/249‑1 (MQ=255)                                                                                                                                                                             gCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGa                                                                  >  1:299335/1‑249 (MQ=255)                                                                                                                                                                                ggATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGAAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGCTTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGACAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTCGTTTGACCGATACCCGCAAAGAATTGc                                                             >  2:122243/1‑251 (MQ=255)                                                                                                                                                                                  aTTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTCGTTTGACCGATCACAGCACAGAATTGCTg                                                           >  2:90552/1‑251 (MQ=255)                                                                                                                                                                                     cACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGtct                                                        <  1:319392/251‑1 (MQ=255)                                                                                                                                                                                               gttgttAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGGCTCATAAGGAGTTTTTCGTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGAGAACAGCAAAGAATTGTTGTCTCGTTTTGaaa                                              <  2:25761/251‑1 (MQ=255)                                                                                                                                                                                                                   agagCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTg                          <  1:76039/248‑1 (MQ=255)                                                                                                                                                                                                                                      ggAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTa        >  2:15382/1‑250 (MQ=255)                                                                                                                                                                                                                                           gACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCATCTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTg  >  2:254077/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              TGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTG  >  NC_000913/2926045‑2926529

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