Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,860,164 +A coding (405/657 nt) pphB → serine/threonine‑specific protein phosphatase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,860,1601.A100.0% 172.5 / NA 46I134K (ATA→AAA) pphBserine/threonine‑specific protein phosphatase 2
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (21/25);  total (21/25)

CGAAACTGGCGATGGCAATATGTGGCTTGCCAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGTTACAACAAATAAACGGCGCTG  >  NC_000913/2860017‑2860300
                                                                                                                                                |                                                                                                                                            
cGAAACTGGCGATGGCAATATGTGGCTTGCCAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATaaaaata                                                                                                                                       <  1:760261/151‑4 (MQ=255)
 gAAACTGGCGATGGCAATATGTGGCTTGCCAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATaaaaatat                                                                                                                                      <  2:572829/151‑5 (MQ=255)
 gAAACTGGCGATGGCAATATGTGGCTTGCCAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATaaaaatat                                                                                                                                      >  1:64794/1‑147 (MQ=255)
   aaCTGGCGATGGCAATATGTGGCTTGCCAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGc                                                                                                                                     <  2:528356/151‑1 (MQ=255)
    aCTGGCGATGGCAATATGTGGCTTGCCAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGcc                                                                                                                                    >  1:458585/1‑151 (MQ=255)
       ggCGATGGCAATATGTGGCTTGCCAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCa                                                                                                                                   <  1:745165/149‑1 (MQ=255)
                   atGTGGCTTGCCAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGc                                                                                                                         >  2:676330/1‑147 (MQ=255)
                             ccAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCgggg                                                                                                            >  1:243107/1‑150 (MQ=255)
                                gCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTg                                                                                                        <  2:506069/151‑1 (MQ=255)
                                   gtggtgACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAat                                                                                                     >  1:114513/1‑151 (MQ=255)
                                                tttCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTa                                                                                         >  2:421708/1‑150 (MQ=255)
                                                  tCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTaaa                                                                                       <  2:596788/150‑1 (MQ=255)
                                                    gATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGaa                                                                                    <  2:644417/151‑1 (MQ=255)
                                                    gATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGa                                                                                     >  1:744292/1‑150 (MQ=255)
                                                       ttAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAg                                                                                                          <  1:665059/126‑1 (MQ=255)
                                                       ttAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAg                                                                                                          >  2:665117/1‑126 (MQ=255)
                                                       ttAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAat                                                                                                     <  1:220777/131‑1 (MQ=255)
                                                       ttAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAat                                                                                                     >  2:220809/1‑131 (MQ=255)
                                                       ttAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAAt                                                                                  <  2:260011/150‑1 (MQ=255)
                                                        tAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATa                                                                                 <  2:593883/150‑1 (MQ=255)
                                                         aaaTGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGc                                                                               <  1:97718/151‑1 (MQ=255)
                                                          aaTGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCg                                                                              <  1:280611/151‑1 (MQ=255)
                                                           aTGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCgg                                                                             >  2:661236/1‑151 (MQ=255)
                                                              aTTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCgg                                                                             >  1:712420/1‑148 (MQ=255)
                                                                 cAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCg                                                                       <  2:495631/151‑1 (MQ=255)
                                                                     gCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAAt                                                                    <  1:491443/150‑1 (MQ=255)
                                                                       aacaaGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTAc                                                                 >  2:235564/1‑151 (MQ=255)
                                                                         caaGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTActc                                                               >  1:707043/1‑151 (MQ=255)
                                                                          aaGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTActct                                                              >  1:195513/1‑151 (MQ=255)
                                                                           agagGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTg                                                             <  2:458636/151‑1 (MQ=255)
                                                                           agagGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTg                                                             >  1:455727/1‑151 (MQ=255)
                                                                                cAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTg                                                        >  2:711846/1‑151 (MQ=255)
                                                                                   tAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTg                                                     >  2:774869/1‑151 (MQ=255)
                                                                                          gTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCgtg                                               <  2:290450/150‑1 (MQ=255)
                                                                                          gTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCgt                                                <  1:703712/149‑1 (MQ=255)
                                                                                              cTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGc                                            >  2:499961/1‑149 (MQ=255)
                                                                                                   aTTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATc                                     >  2:596971/1‑151 (MQ=255)
                                                                                                       cATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCtt                                 >  2:66021/1‑151 (MQ=255)
                                                                                                            ccTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATgg                            <  2:512950/151‑1 (MQ=255)
                                                                                                            ccTTCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATgg                            <  2:741722/151‑1 (MQ=255)
                                                                                                               tCCACATATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGc                           <  2:728846/149‑1 (MQ=255)
                                                                                                                  acaTATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGtt                      <  1:716849/151‑1 (MQ=255)
                                                                                                                   caTATTATTGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGTTa                     <  1:770881/151‑1 (MQ=255)
                                                                                                                        tattGAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGTTacaac                 >  1:742097/1‑150 (MQ=255)
                                                                                                                         attGAAAACACTAACGACAACATAAAAATATGCCATCGCACATGCAGAGTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGTTAcaacaa               <  1:292506/151‑1 (MQ=255)
                                                                                                                            gAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGTTACAACAAAt             <  1:581443/150‑1 (MQ=255)
                                                                                                                            gAAATCACTAACGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGTTACAACAAATa            <  2:37768/151‑1 (MQ=255)
                                                                                                                                     aaCGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGTTACAACAAATAAACGgcg     <  1:21608/149‑1 (MQ=255)
                                                                                                                                      aCGACAACATAAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGTTACAACAAATAAACGGCGCTg  >  1:764300/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                            
CGAAACTGGCGATGGCAATATGTGGCTTGCCAGCGGTGGTGACTGGTTTTTCGATTTAAATGATTCAGAGCAACAAGAGGCAATAGATCTGTTGCTGAAATTCCATCACCTTCCACATATTATTGAAATCACTAACGACAACATAAAATATGCCATCGCACATGCAGATTATCCGGGGAGTGAATATCTCTTTGGTAAAGAAATAGCGGAGAGCGAATTACTCTGGCCTGTTGATCGTGTGCAGAAATCGCTTAATGGCGAGTTACAACAAATAAACGGCGCTG  >  NC_000913/2860017‑2860300

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