Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,176,721 G→A T38I (ACC→ATC)  cpdA ← 3',5' cAMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,7210GA100.0% 183.0 / NA 57T38I (ACC→ATC) cpdA3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/29);  total (28/29)

CTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGAC  >  NC_000913/3176586‑3176856
                                                                                                                                       |                                                                                                                                       
cTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGatgt                                                                                                                                      <  1:1099455/139‑1 (MQ=255)
                 gAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCt                                                                                                                     >  1:492753/1‑139 (MQ=255)
                        tGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCt                                                                                                              >  1:449310/1‑139 (MQ=255)
                           tAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCtttt                                                                                                           <  1:500434/139‑1 (MQ=255)
                               ccGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGc                                                                                                       >  2:18029/1‑139 (MQ=255)
                               ccGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGc                                                                                                       <  2:449310/139‑1 (MQ=255)
                                  cAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCaaa                                                                                                    >  2:1223547/1‑139 (MQ=255)
                                        aTTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGatgt                                                                                                                                      <  2:263554/99‑1 (MQ=255)
                                        aTTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGatgt                                                                                                                                      >  1:263554/1‑99 (MQ=255)
                                        aTTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                              <  2:749003/139‑1 (MQ=255)
                                          tGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAg                                                                                            <  2:492753/139‑1 (MQ=255)
                                                 gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGatgt                                                                                                                                      >  2:326333/1‑90 (MQ=255)
                                                 gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGatgt                                                                                                                                      <  1:326333/90‑1 (MQ=255)
                                                 gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGttt                                                                                                                                    >  1:914742/1‑92 (MQ=255)
                                                 gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGttt                                                                                                                                    <  2:914742/92‑1 (MQ=255)
                                                 gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt                                                                                               >  2:1612312/1‑129 (MQ=255)
                                                 gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt                                                                                               <  1:1612312/129‑1 (MQ=255)
                                                 gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGt                                                                                     >  1:875309/1‑139 (MQ=255)
                                                       aaTCTCCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCtt                                                                                                             <  2:1014615/109‑1 (MQ=255)
                                                       aaTCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCtt                                                                                                             >  1:1014615/1‑109 (MQ=255)
                                                       aaTCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTg                                                                               >  1:619531/1‑139 (MQ=255)
                                                       aaTCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTg                                                                               >  1:796100/1‑139 (MQ=255)
                                                          cACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaa                                                                            <  2:1111877/139‑1 (MQ=255)
                                                            ccTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCg                                                                                       <  2:719630/126‑1 (MQ=255)
                                                            ccTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCg                                                                                       >  1:719630/1‑126 (MQ=255)
                                                                  gcgACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGa                                                                    >  2:925077/1‑139 (MQ=255)
                                                                          cAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGcc                                                            >  1:1490986/1‑139 (MQ=255)
                                                                              tCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAc                                                        >  1:499766/1‑139 (MQ=255)
                                                                              tCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAc                                                        <  1:1223547/139‑1 (MQ=255)
                                                                              tCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAc                                                        <  1:1279865/139‑1 (MQ=255)
                                                                               cGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCa                                                                                                      <  1:1251173/92‑1 (MQ=255)
                                                                               cGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCa                                                                                                      >  2:1251173/1‑92 (MQ=255)
                                                                               cGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAcc                                                       <  2:619531/139‑1 (MQ=255)
                                                                                 aaTTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGatgt                                                                                                                                      >  1:821693/1‑58 (MQ=255)
                                                                                 aaTTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGatgt                                                                                                                                      <  2:821693/58‑1 (MQ=255)
                                                                                  aTTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAgg                                                                                                <  2:868891/95‑1 (MQ=255)
                                                                                  aTTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAgg                                                                                                >  1:868891/1‑95 (MQ=255)
                                                                                        gCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAgg                                              <  2:875309/139‑1 (MQ=255)
                                                                                         cTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTGTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGa                                             <  1:711625/139‑1 (MQ=255)
                                                                                               tgGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCGAGATGCTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTa                                                                             >  2:1079509/1‑101 (MQ=255)
                                                                                               tgGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTa                                                                             <  1:1079509/101‑1 (MQ=255)
                                                                                                   cGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACa                                   <  2:499766/139‑1 (MQ=255)
                                                                                                      aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGttt                                                                                                                                    <  2:1122973/39‑1 (MQ=255)
                                                                                                      aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGttt                                                                                                                                    >  1:1122973/1‑39 (MQ=255)
                                                                                                      aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAATGGTTAACAGGc                                >  2:1047782/1‑139 (MQ=255)
                                                                                                      aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGc                                >  2:755124/1‑139 (MQ=255)
                                                                                                      aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGc                                <  1:1147047/139‑1 (MQ=255)
                                                                                                      aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGc                                <  1:925077/139‑1 (MQ=255)
                                                                                                          ccTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAAAAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTc                            <  1:1047782/139‑1 (MQ=255)
                                                                                                           cTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTcc                           >  1:603985/1‑139 (MQ=255)
                                                                                                             ccAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCaa                         <  2:603985/139‑1 (MQ=255)
                                                                                                                gCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATg                      <  1:18029/139‑1 (MQ=255)
                                                                                                                  aCCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGgt                    <  1:755124/139‑1 (MQ=255)
                                                                                                                    cGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCtt                                                                                                             <  1:1553028/48‑1 (MQ=255)
                                                                                                                    cGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCtt                                                                                                             >  2:1553028/1‑48 (MQ=255)
                                                                                                                    cGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCtt                                                                                                             <  1:87026/48‑1 (MQ=255)
                                                                                                                    cGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCtt                                                                                                             >  2:87026/1‑48 (MQ=255)
                                                                                                                     gCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTc                 >  1:1279477/1‑139 (MQ=255)
                                                                                                                       cTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCt               >  2:281353/1‑139 (MQ=255)
                                                                                                                        tGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCtt              >  1:866671/1‑139 (MQ=255)
                                                                                                                                ctcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCa      <  2:866671/139‑1 (MQ=255)
                                                                                                                                   ccAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGt                                                                              <  1:33246/64‑1 (MQ=255)
                                                                                                                                   ccAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGt                                                                              >  2:33246/1‑64 (MQ=255)
                                                                                                                                    cAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGACGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacgac  >  1:172098/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGAC  >  NC_000913/3176586‑3176856

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: ‑