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_0009132,685,5610TC100.0% 104.1 / NA 29intergenic (‑54/‑274)glyA/hmpserine hydroxymethyltransferase/fused nitric oxide dioxygenase/dihydropteridine reductase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/12);  total (17/12)

TCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAG  >  NC_000913/2685425‑2685664
                                                                                                                                        |                                                                                                       
tctTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGt                                                                                                       >  1:392696/1‑139 (MQ=255)
            aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                           >  2:714163/1‑139 (MQ=255)
            aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                           <  1:714163/139‑1 (MQ=255)
              tttttCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                         <  1:635544/139‑1 (MQ=255)
                                  gCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgt                                                                     >  2:644237/1‑139 (MQ=255)
                                  gCCACAGTTCGGCATCATAATCGGCAATGTCCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgt                                                                     >  2:638009/1‑139 (MQ=255)
                                         ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                              >  2:449950/1‑139 (MQ=255)
                                           cGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                         >  2:24222/1‑110 (MQ=255)
                                           cGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                         <  1:24222/110‑1 (MQ=255)
                                                   tAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                           >  1:464540/1‑100 (MQ=255)
                                                   tAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                           <  2:464540/100‑1 (MQ=255)
                                                        ggCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACg                                                                                              >  1:173431/1‑92 (MQ=255)
                                                        ggCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACg                                                                                              <  2:173431/92‑1 (MQ=255)
                                                            aTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGa                                                                                          <  1:581879/92‑1 (MQ=255)
                                                            aTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGa                                                                                          >  2:581879/1‑92 (MQ=255)
                                                                                   ccGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                    >  2:649259/1‑139 (MQ=255)
                                                                                          tcCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCa             <  1:649259/139‑1 (MQ=255)
                                                                                                   aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGtat                                                                                   <  1:145255/60‑1 (MQ=255)
                                                                                                   aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGtat                                                                                   >  2:145255/1‑60 (MQ=255)
                                                                                                   aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg    >  2:543151/1‑139 (MQ=255)
                                                                                                     cTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGgag  >  1:254334/1‑139 (MQ=255)
                                                                                                             aaaaTTTTTGGCCTTTATAGGCGGTCCCGTTGGACaa                                                                                                >  2:147064/1‑37 (MQ=255)
                                                                                                             aaaaTTTTTGGCCTTTATAGGCGGTCCCGTTGGACaa                                                                                                <  1:147064/37‑1 (MQ=255)
                                                                                                                                     tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                              >  2:381891/1‑47 (MQ=255)
                                                                                                                                     tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                              <  1:381891/47‑1 (MQ=255)
                                                                                                                                     tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                              >  2:132472/1‑79 (MQ=255)
                                                                                                                                     tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                              <  1:132472/79‑1 (MQ=255)
                                                                                                                                     tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCa             <  2:182225/96‑1 (MQ=255)
                                                                                                                                     tCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCa             >  1:182225/1‑96 (MQ=255)
                                                                                                                                        |                                                                                                       
TCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAG  >  NC_000913/2685425‑2685664

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