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% 216.5 / NA 66T143I (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 (31/35);  total (31/35)

GCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTT  >  NC_000913/3587244‑3587504
                                                                                                                              |                                                                                                                                      
gcggGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            <  2:921529/139‑1 (MQ=255)
gcggGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            <  1:798434/139‑1 (MQ=255)
     aTGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGtt                                                                                                                       >  1:1260311/1‑139 (MQ=255)
             cAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCaa                                                                                                               >  2:56586/1‑139 (MQ=255)
                       aaTCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                     <  1:56586/139‑1 (MQ=255)
                       aaTCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                     <  2:652350/139‑1 (MQ=255)
                              aTAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGtt                                                                                                                       <  2:596901/114‑1 (MQ=255)
                              aTAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGtt                                                                                                                       >  1:596901/1‑114 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGATAATTATTGGGCCGTCAAGTTGGGACCAATCACGTCAATCGCATCGGTAc                                                                                          <  1:1331247/139‑1 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                     >  2:361454/1‑128 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                     <  1:361454/128‑1 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCg                                                                                                   >  1:591142/1‑130 (MQ=255)
                                  ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCg                                                                                                   <  2:591142/130‑1 (MQ=255)
                                    cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAg                                                                                        <  2:525386/139‑1 (MQ=255)
                                    cGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAg                                                                                        >  1:481864/1‑139 (MQ=255)
                                         aaTCAGCAGCATATCCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCTCATGCTGACATTCCGTTGAAAATTCTTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                                   <  2:1260311/139‑1 (MQ=255)
                                              gcagcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATg                                                                                     >  1:1367166/1‑132 (MQ=255)
                                              gcagcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATg                                                                                     <  2:1367166/132‑1 (MQ=255)
                                              gcagcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCca                                                                              <  1:1498840/139‑1 (MQ=255)
                                                    tatGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAcccc                                                                        <  1:1913288/139‑1 (MQ=255)
                                                    tatGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAcccc                                                                        <  1:986357/139‑1 (MQ=255)
                                                         tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCgg                                                                                             >  1:1613713/1‑113 (MQ=255)
                                                         tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCgg                                                                                             <  2:1613713/113‑1 (MQ=255)
                                                         tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                                   >  2:1944176/1‑139 (MQ=255)
                                                            cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcagc                                                                >  2:1749995/1‑139 (MQ=255)
                                                              gagaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCaa                                                              <  2:556564/139‑1 (MQ=255)
                                                                      aaaCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            <  1:2055301/69‑1 (MQ=255)
                                                                      aaaCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            >  2:2055301/1‑69 (MQ=255)
                                                                          ggCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCCTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                                  <  1:1749995/139‑1 (MQ=255)
                                                                           gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcgc                                                 >  1:1045497/1‑139 (MQ=255)
                                                                             gACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                                   >  1:1412819/1‑71 (MQ=255)
                                                                             gACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGa                                                                                                                   <  2:1412819/71‑1 (MQ=255)
                                                                                  gCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGtt                                          >  2:146059/1‑139 (MQ=255)
                                                                                    tGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCAc                                                                                                           >  2:1856106/1‑72 (MQ=255)
                                                                                    tGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCAc                                                                                                           <  1:1856106/72‑1 (MQ=255)
                                                                                      aaTCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATc                                                                                               <  1:637341/82‑1 (MQ=255)
                                                                                      aaTCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATc                                                                                               >  2:637341/1‑82 (MQ=255)
                                                                                        tCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                    <  1:1944176/139‑1 (MQ=255)
                                                                                        tCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                    <  2:481864/139‑1 (MQ=255)
                                                                                           gCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCca                                                                              >  2:1170877/1‑94 (MQ=255)
                                                                                           gCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCca                                                                              <  1:1170877/94‑1 (MQ=255)
                                                                                           gCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                    <  2:836799/136‑1 (MQ=255)
                                                                                           gCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                    >  1:836799/1‑136 (MQ=255)
                                                                                            cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            >  2:1272833/1‑47 (MQ=255)
                                                                                            cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            <  1:1272833/47‑1 (MQ=255)
                                                                                            cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGt                                >  1:740540/1‑139 (MQ=255)
                                                                                             ccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCaa                                                              <  1:1718019/108‑1 (MQ=255)
                                                                                             ccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCaa                                                              >  2:1718019/1‑108 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTggg                                                                                                                                 >  2:1243577/1‑38 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTggg                                                                                                                                 <  1:1243577/38‑1 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            >  2:1770639/1‑43 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTa                                                                                                                            <  1:1770639/43‑1 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATg                                                                                     <  2:1557131/82‑1 (MQ=255)
                                                                                                gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATg                                                                                     >  1:1557131/1‑82 (MQ=255)
                                                                                                      cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                    <  2:1753718/125‑1 (MQ=255)
                                                                                                      cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                    >  1:1753718/1‑125 (MQ=255)
                                                                                                      cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGa                      >  2:265756/1‑139 (MQ=255)
                                                                                                      cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGa                      >  2:783939/1‑139 (MQ=255)
                                                                                                                cATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                     >  1:387620/1‑50 (MQ=255)
                                                                                                                cATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                     <  2:387620/50‑1 (MQ=255)
                                                                                                                   tCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg         <  2:1045497/139‑1 (MQ=255)
                                                                                                                          aaaaTTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAcccc                                                                        >  2:1805677/1‑69 (MQ=255)
                                                                                                                          aaaaTTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAcccc                                                                        <  1:1805677/69‑1 (MQ=255)
                                                                                                                          aaaaTTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTgg                                              >  1:194012/1‑95 (MQ=255)
                                                                                                                          aaaaTTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTgg                                              <  2:194012/95‑1 (MQ=255)
                                                                                                                          aaaaTTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCttt  >  2:187640/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                      
GCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTT  >  NC_000913/3587244‑3587504

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