Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,587,370 C→T T143I (ACT→ATT) 
*248* (TAG→TAA
yhhA →
ugpQ ←
DUF2756 family protein
glycerophosphodiester phosphodiesterase, cytosolic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,587,3700CT100.0% 179.3 / NA 55T143I (ACT→ATT) 
*248* (TAG→TAA
yhhA
ugpQ
DUF2756 family protein
glycerophosphodiester phosphodiesterase, cytosolic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/27);  total (28/27)

CAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGG  >  NC_000913/3587239‑3587497
                                                                                                                                   |                                                                                                                               
cAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTggg                                                                                                                          <  1:171997/139‑1 (MQ=255)
  aCGGCGGGATGTTAAACAGCAACCATAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGcc                                                                                                                        >  1:395244/1‑139 (MQ=255)
            gTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCg                                                                                                              >  1:841443/1‑139 (MQ=255)
                 aCAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCa                                                                                                         <  2:841443/139‑1 (MQ=255)
                       aCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                     <  2:356691/121‑1 (MQ=255)
                       aCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                     >  1:356691/1‑121 (MQ=255)
                                            tGAATCAACAGCATATGCTGCCGGAGAGGAGAAACGGCGACTTGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTCAAAATTATTgg                                                                                                                           <  2:874294/94‑1 (MQ=255)
                                            tGAATCAACAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTgg                                                                                                                           >  1:874294/1‑94 (MQ=255)
                                                              tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                            >  1:992493/1‑139 (MQ=255)
                                                                ccGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                          <  2:164817/139‑1 (MQ=255)
                                                                ccGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                          >  2:1336105/1‑139 (MQ=255)
                                                                 cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                     <  2:769289/79‑1 (MQ=255)
                                                                 cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                     >  1:769289/1‑79 (MQ=255)
                                                                 cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcagc                                                         >  2:1185149/1‑139 (MQ=255)
                                                                     gaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                            >  1:400340/1‑84 (MQ=255)
                                                                     gaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                            <  2:400340/84‑1 (MQ=255)
                                                                     gaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAAct                                                     >  2:1010844/1‑139 (MQ=255)
                                                                        gagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                     <  2:1452510/72‑1 (MQ=255)
                                                                        gagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                     >  1:1452510/1‑72 (MQ=255)
                                                                        gagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCa                                                                                                         <  1:604017/84‑1 (MQ=255)
                                                                        gagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCa                                                                                                         >  2:604017/1‑84 (MQ=255)
                                                                           aaaCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTctgctg                                               <  1:694441/139‑1 (MQ=255)
                                                                              cGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                       <  2:1438170/80‑1 (MQ=255)
                                                                              cGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                       >  1:1438170/1‑80 (MQ=255)
                                                                                gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTctg                                                  >  2:826464/1‑131 (MQ=255)
                                                                                gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTctg                                                  <  1:826464/131‑1 (MQ=255)
                                                                                gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcgc                                          <  1:1010844/139‑1 (MQ=255)
                                                                                  gACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                           <  1:93464/136‑1 (MQ=255)
                                                                                  gACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                           >  2:93464/1‑136 (MQ=255)
                                                                                  gACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTg                                        >  2:1052982/1‑139 (MQ=255)
                                                                                     aTGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTgggg                                     >  2:1037769/1‑139 (MQ=255)
                                                                                         tGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTa                                 >  1:847586/1‑139 (MQ=255)
                                                                                               aGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACg                           <  1:1052982/139‑1 (MQ=255)
                                                                                                gCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGt                                                                                                                      >  1:635859/1‑47 (MQ=255)
                                                                                                gCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGt                                                                                                                      <  2:635859/47‑1 (MQ=255)
                                                                                                gCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                     <  2:919875/48‑1 (MQ=255)
                                                                                                gCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                          <  1:1301503/107‑1 (MQ=255)
                                                                                                gCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                          >  2:1301503/1‑107 (MQ=255)
                                                                                                gCCCAGCACGCCGCAGCCTGACATGCCGTTGAAAATTATTGGGCCGTa                                                                                                                     >  1:919875/1‑48 (MQ=38)
                                                                                                 cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTa                                                                                                                               <  1:488232/37‑1 (MQ=255)
                                                                                                     gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATaa                     >  2:62045/1‑139 (MQ=255)
                                                                                                     gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATaaa                    >  1:442071/1‑139 (MQ=255)
                                                                                                           cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                       <  2:931165/51‑1 (MQ=255)
                                                                                                           cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                       >  1:931165/1‑51 (MQ=255)
                                                                                                           cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                            >  2:864117/1‑78 (MQ=255)
                                                                                                           cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                            <  1:864117/78‑1 (MQ=255)
                                                                                                           cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGa               <  2:992493/139‑1 (MQ=255)
                                                                                                                     cATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTc     <  1:1185149/139‑1 (MQ=255)
                                                                                                                        tCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCaa                                                       <  1:989081/86‑1 (MQ=255)
                                                                                                                        tCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCaa                                                       >  2:989081/1‑86 (MQ=255)
                                                                                                                        tCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg  <  1:1336105/139‑1 (MQ=255)
                                                                                                                          cGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                       <  2:39965/36‑1 (MQ=255)
                                                                                                                          cGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                       >  1:39965/1‑36 (MQ=255)
                                                                                                                              gAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                            >  2:275862/1‑59 (MQ=255)
                                                                                                                              gAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                            <  1:275862/59‑1 (MQ=255)
                                                                                                                                   |                                                                                                                               
CAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGG  >  NC_000913/3587239‑3587497

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: