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% 109.0 / NA 35T38I (ACC→ATC) cpdA3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/16);  total (19/16)

GCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATG  >  NC_000913/3176589‑3176836
                                                                                                                                    |                                                                                                                   
gCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTa                                                                                                               >  1:675009/1‑139 (MQ=255)
                     tGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCt                                                                                          <  1:391510/139‑1 (MQ=255)
                             cGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCt                                                                                          >  2:52636/1‑131 (MQ=255)
                             cGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCt                                                                                          <  1:52636/131‑1 (MQ=255)
                                                gcTAAATCACCTGTAGGGACAATCAGGTCGAATTCCTGCTGGGGTGGCTGAATCGCCTCCAGCACCGCCAGGTAACTCTCCTAGATGTTTACCCCTAACAGGCCTTCGTGCTTTTGTGCAAAAAGATGAGGGTTGGTat                                                               >  2:1037465/1‑138 (MQ=255)
                                                    aaTCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCaa                                                                                 >  2:91826/1‑117 (MQ=255)
                                                    aaTCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCaa                                                                                 <  1:91826/117‑1 (MQ=255)
                                                    aaTCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCg                                                                   <  2:821900/131‑1 (MQ=255)
                                                    aaTCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCg                                                                   >  1:821900/1‑131 (MQ=255)
                                                     aTCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGt                                                          >  1:1055345/1‑139 (MQ=255)
                                                            gTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATcc                                                   >  2:1074120/1‑139 (MQ=255)
                                                                cgACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGAc                                               >  2:296337/1‑139 (MQ=255)
                                                                cgACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGAc                                               >  1:900965/1‑139 (MQ=255)
                                                                 gACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGAct                                              <  1:1145267/139‑1 (MQ=255)
                                                                  gcAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGa                                                                         >  1:282140/2‑111 (MQ=255)
                                                                  gcAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGa                                                                         <  2:282140/110‑1 (MQ=255)
                                                                                         gtgtGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAgg                      >  2:1389375/1‑139 (MQ=255)
                                                                                              gCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaa                                                        >  1:1016971/1‑100 (MQ=255)
                                                                                              gCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaa                                                        <  2:1016971/100‑1 (MQ=255)
                                                                                                cGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACa               <  1:296337/139‑1 (MQ=255)
                                                                                                   aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCaaa                                                                                <  2:221021/71‑1 (MQ=255)
                                                                                                   aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCaaa                                                                                >  1:221021/1‑71 (MQ=255)
                                                                                                   aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                          >  1:830750/1‑77 (MQ=255)
                                                                                                   aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGAGCAAActgatg                                                                          <  2:830750/77‑6 (MQ=255)
                                                                                                             gCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATg  <  1:1389375/139‑1 (MQ=255)
                                                                                                                 cGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAggg                                                                                                   <  2:620079/38‑1 (MQ=255)
                                                                                                                 cGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAggg                                                                                                   >  1:620079/1‑38 (MQ=255)
                                                                                                                  gCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTa                                                         <  1:117057/79‑1 (MQ=255)
                                                                                                                  gCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTa                                                         >  2:117057/1‑79 (MQ=255)
                                                                                                                   ccTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATg  >  1:794597/1‑133 (MQ=255)
                                                                                                                   ccTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGCTAACAGGCTTTCCAAATg  <  2:794597/133‑1 (MQ=255)
                                                                                                                              tcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt                                                                           <  2:495468/49‑1 (MQ=255)
                                                                                                                              tcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt                                                                           >  1:495468/1‑49 (MQ=255)
                                                                                                                              tcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGt                                                                 >  2:114840/1‑59 (MQ=255)
                                                                                                                              tcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGt                                                                 <  1:114840/59‑1 (MQ=255)
                                                                                                                                    |                                                                                                                   
GCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATG  >  NC_000913/3176589‑3176836

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