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% 185.2 / NA 58T143I (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 (30/28);  total (30/28)

GCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTT  >  NC_000913/3587244‑3587504
                                                                                                                              |                                                                                                                                      
gcggGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            <  2:1528779/139‑1 (MQ=255)
   ggATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAg                                                                                                                         >  1:1312248/1‑139 (MQ=255)
                cAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCCGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCa                                                                                                            >  1:1013326/1‑139 (MQ=255)
                cAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                              >  1:703412/1‑137 (MQ=255)
                cAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                              <  2:703412/137‑1 (MQ=255)
                 aaCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                                   >  1:802406/1‑131 (MQ=255)
                 aaCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                                   <  2:802406/131‑1 (MQ=255)
                  aCCGTAATCCGGATAGTTCGTTTAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGAc                                                                                                                  <  1:451888/131‑1 (MQ=255)
                  aCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGATGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGAc                                                                                                                  >  2:451888/1‑131 (MQ=255)
                  aCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCaa                                                                                                               >  2:975984/1‑134 (MQ=255)
                  aCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCaa                                                                                                               <  1:975984/134‑1 (MQ=255)
                  aCCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACg                                                                                                          >  1:1415742/1‑139 (MQ=255)
                      tAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                                   >  2:64546/1‑126 (MQ=255)
                      tAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                                   <  1:64546/126‑1 (MQ=255)
                            ggATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCAt                                                                                                <  2:1415742/139‑1 (MQ=255)
                              aTAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTggg                                                                                                                                 >  2:1652377/1‑104 (MQ=255)
                              aTAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTggg                                                                                                                                 <  1:1652377/104‑1 (MQ=255)
                                 gTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTa                                                                                           <  2:1312248/139‑1 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            <  2:179771/105‑1 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            >  1:657369/1‑105 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            >  1:179771/1‑105 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            <  2:657369/105‑1 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                              <  2:92393/119‑1 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAAt                                                                                                              >  1:92393/1‑119 (MQ=255)
                                    cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATc                                                                                               <  1:649191/132‑1 (MQ=255)
                                    cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATc                                                                                               >  2:649191/1‑132 (MQ=255)
                                    cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAg                                                                                        >  2:1351563/1‑139 (MQ=255)
                                       tGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATg                                                                                     >  2:971714/1‑139 (MQ=255)
                                           tCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAAt                                                                                 >  2:453715/1‑139 (MQ=255)
                                                agcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGTAATCcaca                                                                            <  2:1013326/139‑1 (MQ=255)
                                                agcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                                   >  2:686859/1‑132 (MQ=255)
                                                agcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                                   <  1:686859/132‑1 (MQ=255)
                                                   atatGCTGCCGGAGAGGAGAAACGGCGACATGCTAAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGACGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAccc                                                                         <  1:453715/139‑1 (MQ=255)
                                                       gctgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgc                                                                     >  2:312097/1‑139 (MQ=255)
                                                           ccGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                                 <  1:326868/139‑1 (MQ=255)
                                                            cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTggg                                                                                                                                 >  1:281972/1‑74 (MQ=255)
                                                            cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTggg                                                                                                                                 <  2:281972/74‑1 (MQ=255)
                                                            cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            <  1:256268/79‑1 (MQ=255)
                                                            cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            >  2:256268/1‑79 (MQ=255)
                                                                     gaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTTAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTctgct                                                       <  1:971714/139‑1 (MQ=255)
                                                                     gaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                                   <  1:1556683/79‑1 (MQ=255)
                                                                     gaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGATCGGa                                                                                                                   >  2:1556683/1‑79 (MQ=255)
                                                                             gACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTg                                               >  2:384745/1‑139 (MQ=255)
                                                                              aCATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTgg                                              <  1:312097/139‑1 (MQ=255)
                                                                                         cAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGAc                                   >  2:904076/1‑139 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCgg                                                                                                                    >  2:216540/1‑51 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCgg                                                                                                                    <  1:216540/51‑1 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                                  >  2:232674/1‑117 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                                  <  1:232674/117‑1 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATaa                            >  2:606265/1‑139 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATaa                            >  1:826950/1‑139 (MQ=255)
                                                                                                              gACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCaa                                                              <  1:769148/91‑1 (MQ=255)
                                                                                                              gACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCaa                                                              >  2:769148/1‑91 (MQ=255)
                                                                                                                   tCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg         <  1:904076/139‑1 (MQ=255)
                                                                                                                        tGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCt    <  1:243735/139‑1 (MQ=255)
                                                                                                                         gAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCca                                                                              <  2:841048/64‑1 (MQ=255)
                                                                                                                         gAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCca                                                                              >  1:841048/1‑64 (MQ=255)
                                                                                                                          aaaaTTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCttt  >  1:1055568/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                      
GCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTT  >  NC_000913/3587244‑3587504

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