Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,685,561	T→C	intergenic (‑54/‑274)	glyA ← / → hmp	serine hydroxymethyltransferase/fused nitric oxide dioxygenase/dihydropteridine reductase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,685,561	0	T	C	100.0%	242.7 / NA	67	intergenic (‑54/‑274)	glyA/hmp	serine hydroxymethyltransferase/fused nitric oxide dioxygenase/dihydropteridine reductase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/34);  total (33/34)

TTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCT  >  NC_000913/2685427‑2685696                                                                                                                                       |                                                                                                                                        ttCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTg                                                                                                                                     >  1:1639273/1‑139 (MQ=255)   ccTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                   >  2:1616104/1‑139 (MQ=255)        cGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACg                                                                                                                              <  1:406497/139‑1 (MQ=255)        cGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACg                                                                                                                              <  1:1616104/139‑1 (MQ=255)            cTTTTTCCTTCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCTGTCCCGTTGGACAACGGCGa                                                                                                                          <  2:997025/139‑1 (MQ=255)                           tAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAAt                                                                                                           >  2:875668/1‑139 (MQ=255)                                tGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           <  1:1139623/118‑1 (MQ=255)                                tGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           >  2:1139623/1‑118 (MQ=255)                                  ccACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  2:664593/118‑1 (MQ=255)                                  ccACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  2:560418/118‑1 (MQ=255)                                  ccACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  1:560418/1‑118 (MQ=255)                                  ccACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  1:664593/1‑118 (MQ=255)                                        ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              >  2:634645/1‑139 (MQ=255)                                            gCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           >  1:330458/1‑106 (MQ=255)                                            gCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           <  2:330458/106‑1 (MQ=255)                                                  tAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGa                                                                                    >  1:430904/1‑139 (MQ=255)                                                          aaTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  2:806246/1‑94 (MQ=255)                                                          aaTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  1:806246/94‑1 (MQ=255)                                                          aaTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCa                                                                                                         >  2:401821/1‑110 (MQ=255)                                                          aaTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCa                                                                                                         <  1:401821/110‑1 (MQ=255)                                                            tGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  1:1613271/92‑1 (MQ=255)                                                            tGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  2:1613271/1‑92 (MQ=255)                                                                cATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACg                                                                      <  1:875668/139‑1 (MQ=255)                                                                 aTTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAAc                                                                                             <  2:220100/115‑1 (MQ=255)                                                                 aTTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAAc                                                                                             >  1:220100/1‑115 (MQ=255)                                                                        gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           <  1:450341/78‑1 (MQ=255)                                                                        gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           >  2:450341/1‑78 (MQ=255)                                                                        gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtg                                                                                                      >  2:1289001/1‑99 (MQ=255)                                                                        gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtg                                                                                                      <  1:1289001/99‑1 (MQ=255)                                                                        gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                              <  2:1639273/139‑1 (MQ=255)                                                                         cTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAtt                                                             <  1:634645/139‑1 (MQ=255)                                                                           ttAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTATGCAAACGAt                                                              >  1:615477/1‑136 (MQ=255)                                                                           ttAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                              <  2:615477/136‑1 (MQ=255)                                                                             aaCATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtg                                                                                                      <  1:1607594/94‑1 (MQ=255)                                                                             aaCATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCt                                                         >  1:99360/1‑139 (MQ=255)                                                                             aaCATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCt                                                         <  2:286554/139‑1 (MQ=255)                                                                             aaCATCCACATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtg                                                                                                      >  2:1607594/1‑94 (MQ=255)                                                                                  ccGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                    <  1:1416486/139‑1 (MQ=255)                                                                                    gCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  1:645555/1‑68 (MQ=255)                                                                                    gCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  2:645555/68‑1 (MQ=255)                                                                                    gCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTc                                                       <  2:998927/134‑1 (MQ=255)                                                                                    gCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTc                                                       >  1:998927/1‑134 (MQ=255)                                                                                    gCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCt                                                  >  2:1414183/1‑139 (MQ=255)                                                                                    gCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCt                                                  >  2:1254280/1‑139 (MQ=255)                                                                                        ctcCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                   <  1:727789/54‑1 (MQ=255)                                                                                        ctcCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                   >  2:727789/1‑54 (MQ=255)                                                                                         tcCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCa                                             <  1:1414183/139‑1 (MQ=255)                                                                                          ccTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAATAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCaa                                            >  2:1538935/1‑139 (MQ=255)                                                                                                tCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCttt                                      <  2:430904/139‑1 (MQ=255)                                                                                                tCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCttt                                      <  1:1254280/139‑1 (MQ=255)                                                                                                  aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCAtt                                                                                                       >  1:1315771/1‑72 (MQ=255)                                                                                                  aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCAtt                                                                                                       <  2:1315771/72‑1 (MQ=255)                                                                                                  aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              <  2:749062/81‑1 (MQ=255)                                                                                                  aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              >  1:749062/1‑81 (MQ=255)                                                                                                  aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCtt                                                                                      >  2:573135/1‑89 (MQ=255)                                                                                                  aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCtt                                                                                      <  1:573135/89‑1 (MQ=255)                                                                                                    cTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  2:397491/52‑1 (MQ=255)                                                                                                    cTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  1:397491/1‑52 (MQ=255)                                                                                                    cTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              >  1:1535835/1‑79 (MQ=255)                                                                                                    cTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              <  2:1535835/79‑1 (MQ=255)                                                                                                          aTAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACaa                                                                                 <  1:1517041/86‑1 (MQ=255)                                                                                                          aTAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACaa                                                                                 >  2:1517041/1‑86 (MQ=255)                                                                                                                     ggCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAAc                 >  2:1007615/1‑139 (MQ=255)                                                                                                                        cTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCgg              <  1:1007615/139‑1 (MQ=255)                                                                                                                                    tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCt  >  2:410290/1‑139 (MQ=255)                                                                                                                                    tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCt  >  2:1311439/1‑139 (MQ=255)                                                                                                                                    tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCt  <  2:99360/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        TTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCT  >  NC_000913/2685427‑2685696

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