Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,841,897 G→A *136* (TAG→TAA
T80I (ACT→ATT) 
nudG →
ynjH ←
CTP pyrophosphohydrolase, also hydrolyzes 2‑hydroxy‑dATP, 8‑hydroxy‑dGTP, 5‑hydroxy‑CTP, dCTP and 5‑methyl‑dCTP
DUF1496 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,841,8970GA100.0% 150.3 / NA 48*136* (TAG→TAA
T80I (ACT→ATT) 
nudG
ynjH
CTP pyrophosphohydrolase, also hydrolyzes 2‑hydroxy‑dATP, 8‑hydroxy‑dGTP, 5‑hydroxy‑CTP, dCTP and 5‑methyl‑dCTP
DUF1496 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/23);  total (25/23)

ACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGC  >  NC_000913/1841774‑1842023
                                                                                                                           |                                                                                                                              
acaTGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGt                                                                                                                 >  2:1155292/1‑139 (MQ=255)
   tGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtcat                                                                                                              <  1:1096388/139‑1 (MQ=255)
   tGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtcat                                                                                                              >  2:1096388/1‑139 (MQ=255)
        aTCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGc                                                                                                         >  2:999983/1‑139 (MQ=255)
            ggCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGc                                                                                                     >  2:19738/1‑139 (MQ=255)
            ggCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGc                                                                                                     >  1:948847/1‑139 (MQ=255)
                   gtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                              <  2:948847/139‑1 (MQ=255)
                       ggtGCTCACCTGACGAGGCGCTGCAATATCCGCTGGCACCTGCTGACATTCCATTATTAGAGGCGTTTATGGCCTTCCGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGt                                                                                          >  2:893751/1‑139 (MQ=255)
                           cTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTc                                                                                      >  2:1302648/1‑139 (MQ=255)
                                        gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAAt                                                                         >  1:706629/1‑139 (MQ=255)
                                          gcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTa                                                                       >  2:557729/1‑139 (MQ=255)
                                          gcTGCAATATCCGCTGGCCCCCGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTACGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTa                                                                       >  2:465513/1‑139 (MQ=255)
                                           cTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                                 >  2:927710/1‑112 (MQ=255)
                                           cTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                                 <  1:927710/112‑1 (MQ=255)
                                             gCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                                   >  1:1109112/1‑108 (MQ=255)
                                             gCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                                   <  2:1109112/108‑1 (MQ=255)
                                                   tCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCtt                                                              <  1:1155292/139‑1 (MQ=255)
                                                      gCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGc                                                                                                         >  1:1163421/1‑93 (MQ=255)
                                                      gCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGc                                                                                                         <  2:1163421/93‑1 (MQ=255)
                                                        tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGaa                                                         <  2:1226337/139‑1 (MQ=255)
                                                        tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGaa                                                         >  1:1226337/1‑139 (MQ=255)
                                                          gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCt                                                                              <  1:495589/116‑1 (MQ=255)
                                                          gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCt                                                                              >  2:495589/1‑116 (MQ=255)
                                                                      aTTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                           <  1:548997/59‑1 (MQ=255)
                                                                      aTTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                           >  2:548997/1‑59 (MQ=255)
                                                                      aTTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgc                                                                                                <  1:1282162/86‑1 (MQ=255)
                                                                      aTTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgc                                                                                                >  2:1282162/1‑86 (MQ=255)
                                                                        tCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGAt                                         <  2:105967/139‑1 (MQ=255)
                                                                        tCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGAt                                         <  1:893751/139‑1 (MQ=255)
                                                                              attaGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgc                                   >  2:679809/1‑139 (MQ=255)
                                                                               ttaGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                  <  1:182743/139‑1 (MQ=255)
                                                                                  gagGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgcagca                               <  1:1302648/139‑1 (MQ=255)
                                                                                    ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCttt                                                                            >  1:1119586/1‑92 (MQ=255)
                                                                                    ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCttt                                                                            <  2:1119586/92‑1 (MQ=255)
                                                                                    ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                             <  1:465513/139‑1 (MQ=255)
                                                                                      cGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                  <  2:173625/132‑1 (MQ=255)
                                                                                      cGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                  >  1:173625/1‑132 (MQ=255)
                                                                                            tGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                                   >  1:1023357/1‑77 (MQ=255)
                                                                                            tGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                                   <  2:1023357/77‑1 (MQ=255)
                                                                                                    cgcgCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                                   <  2:118158/69‑1 (MQ=255)
                                                                                                    cgcgCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                                   >  1:118158/1‑69 (MQ=255)
                                                                                                               aCCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGc  <  1:19738/139‑1 (MQ=255)
                                                                                                                   gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAtt                                                                                                  <  2:928050/39‑1 (MQ=38)
                                                                                                                   gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAtt                                                                                                  >  1:928050/1‑39 (MQ=38)
                                                                                                                   gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                                   <  2:616621/38‑1 (MQ=255)
                                                                                                                   gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                                   >  1:616621/1‑38 (MQ=38)
                                                                                                                   gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                  <  2:73142/103‑1 (MQ=255)
                                                                                                                   gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                  >  1:73142/1‑103 (MQ=255)
                                                                                                                           |                                                                                                                              
ACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGC  >  NC_000913/1841774‑1842023

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