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% 110.4 / NA 35T143I (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 (19/16);  total (19/16)

GCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCG  >  NC_000913/3587259‑3587499
                                                                                                               |                                                                                                                                 
gCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTcaga                                                                                                              <  1:173238/133‑4 (MQ=255)
gCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTcaga                                                                                                              >  2:173238/1‑130 (MQ=255)
 cAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCa                                                                                                       >  1:849641/1‑139 (MQ=255)
        aaTCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAAt                                                                                                <  1:746953/139‑1 (MQ=255)
        aaTCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAAt                                                                                                <  1:961182/139‑1 (MQ=255)
             ggATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCAt                                                                                           <  2:849641/139‑1 (MQ=255)
              gATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATc                                                                                          >  1:376005/1‑139 (MQ=255)
                   ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAAt                                                                                                         >  2:282841/1‑119 (MQ=255)
                   ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAAt                                                                                                         <  1:282841/119‑1 (MQ=255)
                        tGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTa                                                                                                                                 >  2:918955/1‑90 (MQ=255)
                        tGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATg                                                                                <  2:376005/139‑1 (MQ=255)
                                    atatGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAccc                                                                    >  1:761479/1‑139 (MQ=255)
                                     tatGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                       >  2:734264/1‑87 (MQ=255)
                                     tatGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                       <  1:734264/87‑1 (MQ=255)
                                        gctgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgc                                                                >  1:172622/1‑139 (MQ=255)
                                          tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGc                                                                                             <  1:505352/108‑1 (MQ=255)
                                          tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGc                                                                                             >  2:505352/1‑108 (MQ=255)
                                          tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                              >  2:588191/1‑139 (MQ=255)
                                            ccGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                            <  2:147087/139‑1 (MQ=255)
                                                            gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCcaca                                                                       <  1:506607/112‑1 (MQ=255)
                                                            gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCcaca                                                                       >  2:506607/1‑112 (MQ=255)
                                                             cGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCt                                           <  2:172622/139‑1 (MQ=255)
                                                                         tCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                               <  2:761479/139‑1 (MQ=255)
                                                                             cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCCACTCTGCTGCGCGCTGGGGTTTATTGACGGt                           >  1:918892/1‑139 (MQ=255)
                                                                             cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                             <  2:199893/121‑1 (MQ=255)
                                                                             cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                             >  1:199893/1‑121 (MQ=255)
                                                                                 gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATaa                       >  2:995609/1‑139 (MQ=255)
                                                                                       cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                              >  2:750279/1‑78 (MQ=255)
                                                                                       cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                              <  1:750279/78‑1 (MQ=255)
                                                                                            ccTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCg            <  2:735695/139‑1 (MQ=255)
                                                                                                    tCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg    >  2:546399/1‑139 (MQ=255)
                                                                                                    tCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg    <  1:588191/139‑1 (MQ=255)
                                                                                                      cGTTGAAAATTTTTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCg  >  1:360788/1‑139 (MQ=255)
                                                                                                      cGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCg  >  2:140139/1‑139 (MQ=255)
                                                                                                      cGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCg  >  1:1013672/1‑139 (MQ=255)
                                                                                                               |                                                                                                                                 
GCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCG  >  NC_000913/3587259‑3587499

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