Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 874,140 T→C S388P (TCT→CCT)  pdeI → putative c‑di‑GMP phosphodiesterase PdeI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913874,1400TC100.0% 105.3 / NA 31S388P (TCT→CCT) pdeIputative c‑di‑GMP phosphodiesterase PdeI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/16);  total (15/16)

TATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGA  >  NC_000913/874016‑874266
                                                                                                                            |                                                                                                                              
tATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACaat                                                                                                                        >  1:364091/1‑133 (MQ=255)
 aTCAGCATCACAAATGCGATCGCGTAAACGAACGATGCCGACATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACaata                                                                                                                       <  2:129752/133‑1 (MQ=255)
   cAGCATCACAAATTCGATAGCGGAAATGAATGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACaataat                                                                                                                     <  1:244187/133‑1 (MQ=255)
     gCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCa                                                                                                                   <  1:333573/133‑1 (MQ=255)
     gCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGCTTTATTAAGAATGTGGCGATTATTGCGCCTAACAATAATCa                                                                                                                   >  2:176075/1‑133 (MQ=255)
        tCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTa                                                                                                                >  2:79954/1‑133 (MQ=255)
        tCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTa                                                                                                                >  2:447375/1‑133 (MQ=255)
        tCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATAAGGCGATTCTTGTGCATAACAAGAATCAGTa                                                                                                                >  1:152460/1‑133 (MQ=255)
                    tCGCGGAAATGAAAGTTGCCGAAATTGAACCGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTggg                                                                                                    <  1:176075/133‑1 (MQ=255)
                    tCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTggg                                                                                                    <  1:369559/133‑1 (MQ=255)
                      gcgGAAATGAAAGATTCCGAAATTGAACGGTCCAATAAATTTCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGa                                                                                                  <  2:167707/133‑1 (MQ=255)
                      gcgGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGa                                                                                                  <  1:106672/133‑1 (MQ=255)
                                   aTGCCGAAATTGAACCGTCAAATATATTACTCGCACTGATCCCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACtgtg                                                                                     <  2:174021/133‑1 (MQ=255)
                                      ccGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTgg                                                                                  >  1:480918/1‑133 (MQ=255)
                                        gAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTc                                                                                >  2:510871/1‑133 (MQ=255)
                                                  ggTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGt                                                                      <  1:510871/133‑1 (MQ=255)
                                                       aataaatTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAg                                                                 >  1:878/1‑133 (MQ=255)
                                                         taaatTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCg                                                               <  2:286988/133‑1 (MQ=255)
                                                              tGCTCTCACTGATCTCTTACGATCATGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGtt                                                          <  2:289095/133‑1 (MQ=255)
                                                              tGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGtt                                                          <  1:104102/133‑1 (MQ=255)
                                                              tGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGtt                                                          <  1:268936/133‑1 (MQ=255)
                                                              tGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACCGTGTGGTCTGGCAGCCGTGGAGGCGGTGtt                                                          >  1:251922/1‑133 (MQ=255)
                                                                  ctcACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGgt                                                      >  1:357859/1‑133 (MQ=255)
                                                                     aCTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGtt                                                   >  1:363494/1‑133 (MQ=255)
                                                                          tctcTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTtgatga                                              >  1:122242/1‑133 (MQ=255)
                                                                               tACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAg                                             <  1:57556/129‑1 (MQ=255)
                                                                               tACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAg                                             >  2:57556/1‑129 (MQ=255)
                                                                                        gAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGAGTTCATAAGCGCACTAAAATCg                                >  1:193158/1‑133 (MQ=255)
                                                                                            gTGGTTTTATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAg                            <  1:447375/133‑1 (MQ=255)
                                                                                                ttttATTAAAAATATGGCGATTATTGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGaaa                        >  2:348552/1‑133 (MQ=255)
                                                                                                                      attGAGCCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCgaga  <  2:363494/133‑1 (MQ=255)
                                                                                                                            |                                                                                                                              
TATCAGCATCAAAAATTCGATCGCGGAAATGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGA  >  NC_000913/874016‑874266

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