Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 19,009 C→T 19.4% S222F (TCC→TTC)  yybT → DHH subfamily phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073019,0090CT19.4% 119.8 / 20.8 57S222F (TCC→TTC) yybTDHH subfamily phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (18/27);  new base T (3/8);  total (21/36)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGAAGG  >  CP000730/18869‑19147
                                                                                                                                            |                                                                                                                                          
tGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTc                                                                                                                                            <  2:782540/141‑1 (MQ=255)
tGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTc                                                                                                                                            >  2:204720/1‑141 (MQ=255)
 gCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTcc                                                                                                                                           <  2:137193/141‑1 (MQ=255)
 gCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTcc                                                                                                                                           <  2:213985/141‑1 (MQ=255)
        tattTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAAt                                                                                                                                    >  1:275445/1‑141 (MQ=255)
              tAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAg                                                                                                                              <  2:343864/141‑1 (MQ=255)
               aGATAACTACGATGAGATTACGCAAAATATGAAGGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGc                                                                                                                             <  2:610143/141‑1 (MQ=255)
                     ttaagATGAGAATACGCAAAATATGAATGATAGGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTCttt                                                                                                                       <  2:445994/137‑1 (MQ=255)
                         gATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAAt                                                                                                                   <  1:234265/141‑1 (MQ=255)
                         gATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAAt                                                                                                                   <  1:418088/141‑1 (MQ=255)
                               aTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCaaaaa                                                                                                             >  2:338325/1‑141 (MQ=255)
                                 tACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAat                                                                                                           <  2:748606/141‑1 (MQ=255)
                                  aCGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGGCTATTTAAATCAAAAAata                                                                                                          <  2:36251/141‑1 (MQ=255)
                                       aaaTATGAATGATACGCAGCGTTCAGAAATCAACTTAATGGTAACGTGTGACAGTAGTTGATGGGCAACTGAGTATATTATATATTTCAAAAGATACTGTTACGATCAATTCGTAGTCAAAATAAAAAAAAAAATATTAGc                                                                                                     <  2:174216/141‑1 (MQ=255)
                                        aaTATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCt                                                                                                    >  2:63105/1‑141 (MQ=255)
                                        aaTATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCt                                                                                                    <  2:30837/141‑1 (MQ=255)
                                           atgaatgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGAc                                                                                                 <  1:556035/141‑1 (MQ=255)
                                               atgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTag                                                                                             >  1:190026/1‑141 (MQ=255)
                                               atgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTag                                                                                             >  2:594449/1‑141 (MQ=255)
                                               atgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTag                                                                                             >  2:681293/1‑141 (MQ=255)
                                               atgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGAGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTag                                                                                             >  2:674750/1‑141 (MQ=255)
                                                  aTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTagaag                                                                                          <  1:659789/141‑1 (MQ=255)
                                                  aTACGCAGCGTTCAGAAATAAACGCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTagaag                                                                                          <  2:548875/141‑1 (MQ=255)
                                                    aCGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAgaagaa                                                                                        <  1:127697/141‑1 (MQ=255)
                                                      gCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATc                                                                                      <  2:367346/141‑1 (MQ=255)
                                                            ttCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAAtt                                                                                <  1:332637/141‑1 (MQ=255)
                                                                gAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAACATATTAGCTGACTTAGACGAATCTAAAtttaaa                                                                            >  2:68660/1‑138 (MQ=255)
                                                                gAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGat                                                                            <  2:184301/141‑1 (MQ=255)
                                                                gAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGat                                                                            <  1:444239/141‑1 (MQ=255)
                                                                gAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCATATTAAAAAAAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGat                                                                            <  2:184376/141‑1 (MQ=255)
                                                                        cTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACATTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAg                                                                    <  1:782101/141‑1 (MQ=255)
                                                                        cTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAg                                                                    <  1:782182/141‑1 (MQ=255)
                                                                               gTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTTTTTTTTTTTTGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTa                                                             >  2:776273/1‑141 (MQ=255)
                                                                               gTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTa                                                             >  2:542892/1‑141 (MQ=255)
                                                                                  aCGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGt                                                          >  1:9847/1‑141 (MQ=255)
                                                                                    gcgTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGa                                                        >  1:82840/1‑141 (MQ=255)
                                                                                     cgTGTCATTAGTGGATGGGAAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATAAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGaa                                                       <  2:56190/141‑1 (MQ=255)
                                                                                     cgTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTaaaaaaaaaaaa                                                       <  2:224261/141‑10 (MQ=255)
                                                                                      gtgtCATTAGTCGATGGGCAACTGAGTATAATATATTTTACAAAAGATAAAGTTCCGATCAATTCGTAGCCTAATAAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAAtaaagaaaaa                                                      <  2:601902/141‑10 (MQ=255)
                                                                                           aTTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGaaaaagtg                                                 <  2:269629/141‑4 (MQ=255)
                                                                                              aGTCGAAGGGCAACTGAGTATAATAAATTTTTAAAAAGATAAAGTTCCGATCAATTCGTAGAATAATAAAAAAAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGtt                                              <  2:254530/141‑1 (MQ=255)
                                                                                                    tGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTAt                                        >  1:350033/1‑141 (MQ=255)
                                                                                                    tGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTAt                                        >  2:27076/1‑141 (MQ=255)
                                                                                                     gggCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATc                                       >  1:39311/1‑141 (MQ=255)
                                                                                                         aaCTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGc                                   <  1:782650/141‑1 (MQ=255)
                                                                                                         aaCTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGc                                   <  1:782565/141‑1 (MQ=255)
                                                                                                           cagAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATGCGAATCAAAAAATATTAGCTGACTTAGAAGAATCAAAATTTGATATCTTGATGCAATTACGTGAAAAAAGTGTTGGTTATCGTGcac                                 >  2:274067/3‑139 (MQ=255)
                                                                                                            tGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTNAAAAAAGTGTTGGTTATCGTGCCCa                                <  1:241323/141‑1 (MQ=255)
                                                                                                            tGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTAAAAAAAAAAAAAAAAGTGTTGGTTATCGTGCCCa                                <  2:275445/141‑1 (MQ=255)
                                                                                                                 ataataTATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTaa                           >  1:212233/1‑141 (MQ=255)
                                                                                                                   aataTATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACa                         >  1:347922/1‑141 (MQ=255)
                                                                                                                      atatTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTa                      <  2:97018/141‑1 (MQ=255)
                                                                                                                          ttttCAAAAGATACAGTTTCGATCAATTCGTAGACTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTAAATGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTa                  <  2:363396/141‑1 (MQ=255)
                                                                                                                           tttCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTaa                 >  2:418725/1‑140 (MQ=255)
                                                                                                                           tttCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTAt                 <  1:327014/141‑1 (MQ=255)
                                                                                                                             tCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCg               >  1:627855/1‑141 (MQ=255)
                                                                                                                                 tttttttatttCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGgtgt           >  2:50966/10‑141 (MQ=255)
                                                                                                                                    aTACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACATAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTgg        <  2:98187/141‑1 (MQ=255)
                                                                                                                                      aCAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTg      <  2:451255/141‑1 (MQ=255)
                                                                                                                                        aGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGaa    <  1:132036/141‑1 (MQ=255)
                                                                                                                                          ttCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGAAgg  <  2:632805/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
TGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGAAGG  >  CP000730/18869‑19147

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