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%	218.1 / NA	60	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 (32/28);  total (32/28)

CGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCT  >  NC_000913/2685434‑2685696                                                                                                                                |                                                                                                                                        cGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACg                                                                                                                              <  1:952738/139‑1 (MQ=255)  gTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACgg                                                                                                                             >  1:544330/1‑139 (MQ=255)    aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           >  2:538206/1‑139 (MQ=255)           ccTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  1:1137342/134‑1 (MQ=255)           ccTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  2:1137342/1‑134 (MQ=255)           ccTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGta                                                                                                                    >  2:1293409/1‑139 (MQ=255)               cTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAAc                                                                                                                >  2:1300567/1‑139 (MQ=255)                         tGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtg                                                                                                      >  2:234848/1‑139 (MQ=255)                                 ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              <  2:43679/139‑1 (MQ=255)                                 ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              <  2:544330/139‑1 (MQ=255)                                      catcatAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGt                                                                                         <  1:234848/139‑1 (MQ=255)                                           tAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGa                                                                                    >  2:1563284/1‑139 (MQ=255)                                           tAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGa                                                                                    >  2:825883/1‑139 (MQ=255)                                              tCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACaa                                                                                 <  1:1563284/139‑1 (MQ=255)                                                   aaTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  1:3567/1‑94 (MQ=255)                                                   aaTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  2:3567/94‑1 (MQ=255)                                                          aTTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGc                                                                     <  1:748588/139‑1 (MQ=255)                                                                cGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGa                                                               >  1:223844/1‑139 (MQ=255)                                                                cGCTTAAACATCCGAACCTCTGGACTCAGCT‑ACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                              <  1:1167428/139‑1 (MQ=255)                                                                 gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                   <  2:653157/70‑1 (MQ=255)                                                                 gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                   >  1:653157/1‑70 (MQ=255)                                                                 gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                              <  1:825883/139‑1 (MQ=255)                                                                     tAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGc                                                                     >  1:317403/1‑128 (MQ=255)                                                                     tAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGc                                                                     <  2:317403/128‑1 (MQ=255)                                                                      aaCATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACa                                                                                                                                 <  1:530047/67‑1 (MQ=255)                                                                      aaCATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTAGCCTTTATAGGCGGTCCCGTTGGACa                                                                                                                                 >  2:530047/1‑67 (MQ=255)                                                                      aaCATCCGCATCTCCTGACTCAGCT‑ACAATAAAATTTTTGGCCTTTTTAGGCGGTCCCGTTGGACAACCGCGGACAGTATAACCGAAACATTGTGCGATAACAGGTCTTGACCAAGGGATTTAAGCAAACGGTTTccct                                                        >  2:1167428/1‑137 (MQ=255)                                                                           ccGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGt                                                                                                                     >  2:873755/1‑74 (MQ=255)                                                                           ccGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGt                                                                                                                     <  1:873755/74‑1 (MQ=255)                                                                             gCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCt                                                  >  2:1573396/1‑139 (MQ=255)                                                                               aTCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              <  2:1574267/93‑1 (MQ=255)                                                                               aTCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              >  1:1574267/1‑93 (MQ=255)                                                                                     tGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAgg                                          <  1:1224066/139‑1 (MQ=255)                                                                                       aCTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCt                                        <  1:1300567/139‑1 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  1:965684/54‑1 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  2:965684/1‑54 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGa                                                                                                                          >  2:403301/1‑53 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGa                                                                                                                          <  1:403301/53‑1 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGa                                                                                                                          >  1:192135/1‑53 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGa                                                                                                                          <  2:192135/53‑1 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    >  2:1147886/1‑139 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    >  1:1150478/1‑139 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    >  1:545445/1‑139 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTg                                     >  1:1284281/1‑138 (MQ=255)                                                                                           aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTg                                     <  2:1284281/138‑1 (MQ=255)                                                                                             cTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGgag                                  >  2:1439476/1‑139 (MQ=255)                                                                                             cTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGgag                                  >  1:1432044/1‑139 (MQ=255)                                                                                                  aaTAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGACGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATaa                             >  1:1132779/1‑139 (MQ=255)                                                                                                           tttGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATGATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGc                    <  1:1147886/139‑1 (MQ=255)                                                                                                                  tttATAGGCGGTCCCGTTGTACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAAATTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGa             <  2:1132779/139‑1 (MQ=255)                                                                                                                      taGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACg         <  1:538206/139‑1 (MQ=255)                                                                                                                         gCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGAt      <  1:1640477/139‑1 (MQ=255)                                                                                                                            gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTa                                                                        >  2:276696/1‑70 (MQ=255)                                                                                                                            gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTa                                                                        <  1:276696/70‑1 (MQ=255)                                                                                                                            gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    <  2:1493182/106‑1 (MQ=255)                                                                                                                            gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    >  1:1493182/1‑106 (MQ=255)                                                                                                                            gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTc   >  2:101880/1‑139 (MQ=255)                                                                                                                            gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTc   >  1:331481/1‑139 (MQ=255)                                                                                                                             tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACg                                               <  2:718753/94‑1 (MQ=255)                                                                                                                             tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCt  >  1:687761/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                        CGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCT  >  NC_000913/2685434‑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: ‑