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% 123.9 / NA 40T143I (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 (18/22);  total (18/22)

AACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCT  >  NC_000913/3587234‑3587502
                                                                                                                                        |                                                                                                                                    
aaCACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTa                                                                                                                                    <  1:1293501/139‑1 (MQ=255)
          gcggGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          <  2:797333/139‑1 (MQ=255)
          gcggGATGTTAAACAGCAACCGTAATCCGGATAGTGCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          <  2:537035/139‑1 (MQ=255)
                          cAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCa                                                                                                          >  2:950723/1‑139 (MQ=255)
                            aCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACg                                                                                                        >  2:1270998/1‑139 (MQ=255)
                                    ccGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGc                                                                                                <  2:1218790/139‑1 (MQ=255)
                                            ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTAc                                                                                        >  1:709940/1‑139 (MQ=255)
                                              cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAg                                                                                      >  1:549735/1‑139 (MQ=255)
                                                     tCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAAt                                                                               >  2:1023074/1‑139 (MQ=255)
                                                              tatGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAcccc                                                                      >  2:414528/1‑139 (MQ=255)
                                                                     ccGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                               <  1:950723/139‑1 (MQ=255)
                                                                      cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAattatt                                                                                                                                  >  2:512214/1‑71 (MQ=255)
                                                                      cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAattatt                                                                                                                                  <  1:512214/71‑1 (MQ=255)
                                                                          gaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAAct                                                          <  2:549735/139‑1 (MQ=255)
                                                                                   cGGCGACATGCTGAATCAGCCCAGCACTCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          >  2:1035085/1‑66 (MQ=255)
                                                                                   cGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          <  1:1035085/66‑1 (MQ=255)
                                                                                    ggCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTggg                                                                                                                               >  1:371959/1‑60 (MQ=255)
                                                                                    ggCGACATGCTGAAGCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTggg                                                                                                                               <  2:371959/60‑1 (MQ=255)
                                                                                     gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcgc                                               <  1:414528/139‑1 (MQ=255)
                                                                                     gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcgc                                               >  1:1591524/1‑139 (MQ=255)
                                                                                      cGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCt                                              <  2:1234241/139‑1 (MQ=255)
                                                                                           tGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGt                                         >  1:1698685/1‑139 (MQ=255)
                                                                                                  tCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  1:1270998/139‑1 (MQ=255)
                                                                                                  tCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  2:1591524/139‑1 (MQ=255)
                                                                                                    aGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACg                                <  1:1625270/139‑1 (MQ=255)
                                                                                                      cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          >  1:1396389/1‑47 (MQ=255)
                                                                                                      cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                          <  2:1396389/47‑1 (MQ=255)
                                                                                                      cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTcaga                                                                                                                 <  2:442062/56‑4 (MQ=255)
                                                                                                      cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGa                                                                                                                 >  1:442062/1‑56 (MQ=255)
                                                                                                      cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCAt                                                                                              >  1:1510207/1‑75 (MQ=255)
                                                                                                      cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCAt                                                                                              <  2:1510207/75‑1 (MQ=255)
                                                                                                       ccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGta                             <  2:686751/139‑1 (MQ=255)
                                                                                                           cACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATaaa                         >  1:370028/1‑139 (MQ=255)
                                                                                                                cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  1:1245300/125‑1 (MQ=255)
                                                                                                                cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  >  2:1245300/1‑125 (MQ=255)
                                                                                                                        gACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAg            >  2:1276973/1‑139 (MQ=255)
                                                                                                                             tCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg       <  1:1276973/139‑1 (MQ=255)
                                                                                                                               cGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                               <  2:627824/81‑1 (MQ=255)
                                                                                                                               cGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                               >  1:627824/1‑81 (MQ=255)
                                                                                                                                 ttGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  >  1:805746/1‑108 (MQ=255)
                                                                                                                                 ttGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  2:805746/108‑1 (MQ=255)
                                                                                                                                  tGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCt  <  2:679313/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
AACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCT  >  NC_000913/3587234‑3587502

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