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%	223.5 / NA	61	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 (31/30);  total (31/30)

TGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTC  >  NC_000913/2685431‑2685695                                                                                                                                    |                                                                                                                                       tGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACa                                                                                                                                <  2:1397193/139‑1 (MQ=255)  gACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACac                                                                                                                               >  2:290164/1‑138 (MQ=255)    cGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACg                                                                                                                             <  1:1376563/139‑1 (MQ=255)       aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          >  1:757770/1‑139 (MQ=255)       aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          <  1:290164/139‑1 (MQ=255)       aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          <  2:973628/139‑1 (MQ=255)                       tAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGaa                                                                                                           <  1:28071/139‑1 (MQ=255)                       tAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAAt                                                                                                          >  1:1126157/1‑139 (MQ=255)                             gCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgt                                                                                                    >  2:1330919/1‑139 (MQ=255)                               cacaGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGc                                                                                                  >  1:543133/1‑139 (MQ=255)                               cacaGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGc                                                                                                  >  2:200887/1‑139 (MQ=255)                               cacaGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTAACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGc                                                                                                  >  1:1349524/1‑139 (MQ=255)                                    ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          <  1:1140193/110‑1 (MQ=255)                                    ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          >  2:1140193/1‑110 (MQ=255)                                    ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             >  2:881289/1‑139 (MQ=255)                                    ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             <  1:1330919/139‑1 (MQ=255)                                     tCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAAc                                                                                            <  2:543133/139‑1 (MQ=255)                                         catcatAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGt                                                                                        >  2:284952/1‑139 (MQ=255)                                            catAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCtt                                                                                     <  1:962421/139‑1 (MQ=255)                                               aaTCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGAc                                                                                  <  2:757770/139‑1 (MQ=255)                                                        tGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGa                                                                                                                                  >  1:1090915/1‑82 (MQ=255)                                                        tGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGa                                                                                                                                  <  2:1090915/82‑1 (MQ=255)                                                                  aCGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACg                                                                                                                             >  2:944849/1‑77 (MQ=255)                                                                  aCGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACg                                                                                                                             <  1:944849/77‑1 (MQ=255)                                                                  aCGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACg                                                               >  1:913391/1‑139 (MQ=255)                                                                    gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCg                                                                                                 <  2:1382418/103‑1 (MQ=255)                                                                    gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCg                                                                                                 >  1:1382418/1‑103 (MQ=255)                                                                         aaCATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          >  2:1241802/1‑73 (MQ=255)                                                                         aaCATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          <  1:1241802/73‑1 (MQ=255)                                                                           cATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTc                                                      <  1:881289/139‑1 (MQ=255)                                                                              ccGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   <  1:200887/139‑1 (MQ=255)                                                                                      ccTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             <  1:749718/89‑1 (MQ=255)                                                                                      ccTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             >  2:749718/1‑89 (MQ=255)                                                                                        tGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGa                                                                                                                         <  2:288198/59‑1 (MQ=255)                                                                                        tGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGa                                                                                                                         >  1:288198/1‑59 (MQ=255)                                                                                              aGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                   >  1:572134/1‑139 (MQ=255)                                                                                              aGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                   >  2:258887/1‑139 (MQ=255)                                                                                                cTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          >  1:188201/1‑50 (MQ=255)                                                                                                cTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          <  2:188201/50‑1 (MQ=255)                                                                                                cTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACg                                                                     <  1:722634/103‑1 (MQ=255)                                                                                                cTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACg                                                                     >  2:722634/1‑103 (MQ=255)                                                                                                cTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCttt                                     >  2:1041549/1‑135 (MQ=255)                                                                                                cTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCttt                                     <  1:1041549/135‑1 (MQ=255)                                                                                                cTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCt                                       >  2:572139/1‑133 (MQ=255)                                                                                                cTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCt                                       <  1:572139/133‑1 (MQ=255)                                                                                                  aaCAATAAAATTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGaa                               <  2:1349524/139‑1 (MQ=255)                                                                                                          aaTTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgtg                                                                                                   >  1:446469/1‑63 (MQ=255)                                                                                                          aaTTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgtg                                                                                                   <  2:446469/63‑1 (MQ=255)                                                                                                          aaTTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             <  2:934350/69‑1 (MQ=255)                                                                                                          aaTTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             >  1:934350/1‑69 (MQ=255)                                                                                                           aTTTTTGGCC‑TTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCt                      <  1:258887/139‑1 (MQ=255)                                                                                                                                gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTATGCAAACGAt                                                             >  2:1124902/1‑80 (MQ=255)                                                                                                                                gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTATGCAAACGAt                                                             <  1:1124902/80‑1 (MQ=255)                                                                                                                                gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCaaa                                                                 <  2:1442308/76‑1 (MQ=255)                                                                                                                                gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCaaa                                                                 >  1:1442308/1‑76 (MQ=255)                                                                                                                                gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATaa                            <  1:1089661/113‑1 (MQ=255)                                                                                                                                gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATaa                            >  2:1089661/1‑113 (MQ=255)                                                                                                                                gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTc  >  2:1417935/1‑139 (MQ=255)                                                                                                                                gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTc  >  1:647841/1‑139 (MQ=255)                                                                                                                                 tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgtg                                                                                                   <  2:756364/41‑1 (MQ=255)                                                                                                                                 tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   >  1:284172/1‑89 (MQ=255)                                                                                                                                    |                                                                                                                                       TGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCC‑TTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTC  >  NC_000913/2685431‑2685695

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