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%	185.1 / NA	51	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 (29/22);  total (29/22)

TCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCT  >  NC_000913/2685425‑2685696                                                                                                                                         |                                                                                                                                        tctTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGt                                                                                                                                       >  1:640518/1‑139 (MQ=255)       tGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACa                                                                                                                                 <  1:279323/139‑1 (MQ=255)       tGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCAGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACa                                                                                                                                 <  2:449592/139‑1 (MQ=255)         aCGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCACATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAAc                                                                                                                               >  1:161867/1‑139 (MQ=255)             aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           >  2:674447/1‑139 (MQ=255)             aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           >  2:152814/1‑139 (MQ=255)             aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           <  2:663759/139‑1 (MQ=255)             aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           <  2:222108/139‑1 (MQ=255)             aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           <  2:228707/139‑1 (MQ=255)               tttttCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  1:798853/139‑1 (MQ=255)                  ttCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAg                                                                                                                      <  2:534848/139‑1 (MQ=255)                     cTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGtat                                                                                                                   <  1:595086/139‑1 (MQ=255)                             tAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                   >  1:547177/1‑115 (MQ=255)                             tAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                   <  2:547177/115‑1 (MQ=255)                             tAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           <  2:166259/123‑1 (MQ=255)                             tAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                           >  1:166259/1‑123 (MQ=255)                                ccTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCACCTCCTGACTCAGCTAACAATAGAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCAt                                                                                                        <  2:475207/139‑1 (MQ=255)                                  tGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtg                                                                                                      <  1:49134/139‑1 (MQ=255)                                    ccACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgtg                                                                                                    >  1:67008/1‑139 (MQ=255)                                     cacaGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGc                                                                                                   >  1:309379/1‑139 (MQ=255)                                         gTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATa                                                                                               >  2:610910/1‑139 (MQ=255)                                          ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                              <  1:152814/139‑1 (MQ=255)                                              gCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAgg                                                                                          >  1:139973/1‑139 (MQ=255)                                                    tAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  2:7710/1‑102 (MQ=255)                                                    tAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  1:7710/102‑1 (MQ=255)                                                         ggCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAg                                                                               >  2:641232/1‑139 (MQ=255)                                                              tGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAAtt                                                                          >  1:231607/1‑139 (MQ=255)                                                               gTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAAttt                                                                         >  2:718351/1‑139 (MQ=255)                                                                   aTTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGc                                                                     <  2:67008/139‑1 (MQ=255)                                                                     ttCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCaa                                                                   >  2:199786/1‑139 (MQ=255)                                                                      tCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCaaa                                                                  <  1:70927/139‑1 (MQ=255)                                                                        aCGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCaaagg                                                                >  1:308469/1‑137 (MQ=255)                                                                          gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         <  2:768475/80‑1 (MQ=255)                                                                          gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                         >  1:768475/1‑80 (MQ=255)                                                                          gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                              >  2:582439/1‑139 (MQ=255)                                                                                aCATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCtt                                                        <  1:199786/139‑1 (MQ=255)                                                                                               gACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGACTACGCAAGGc                                         >  2:519398/1‑139 (MQ=255)                                                                                                    aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGATAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    >  1:605429/1‑139 (MQ=255)                                                                                                    aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGTCTTTgg                                    >  2:189952/1‑139 (MQ=255)                                                                                                    aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    <  2:640518/139‑1 (MQ=255)                                                                                                    aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    >  1:630525/1‑139 (MQ=255)                                                                                                    aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    >  2:549679/1‑139 (MQ=255)                                                                                                    aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    >  1:262168/1‑139 (MQ=255)                                                                                                    aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                    <  1:519398/139‑1 (MQ=255)                                                                                                     gCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGga                                   <  2:320774/139‑1 (MQ=255)                                                                                                         aCAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAAt                               >  2:724348/1‑139 (MQ=255)                                                                                                                             tataGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGaaa           >  1:95719/1‑139 (MQ=255)                                                                                                                                    ggTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGAttt    >  2:420564/1‑139 (MQ=255)                                                                                                                                     gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCttt                                      >  1:75563/1‑104 (MQ=255)                                                                                                                                     gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCttt                                      <  2:75563/104‑1 (MQ=255)                                                                                                                                      tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCt  >  2:111925/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                        TCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCT  >  NC_000913/2685425‑2685696

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