Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,231,470 A→C 43.9% M39L (ATG→CTG)  prpC → possible PP2C protein phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,231,4700AC43.9% 15.7 / 66.3 59M39L (ATG→CTG) prpCpossible PP2C protein phosphatase
Reads supporting (aligned to +/- strand):  ref base A (18/15);  new base C (16/10);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATATAAATTTTC  >  CP000730/1231332‑1231608
                                                                                                                                           |                                                                                                                                          
cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATg                                                                                                                                          >  2:538330/1‑141 (MQ=255)
    tAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCCAACTAATCAACAACTTTTAGTTGTGAGTG‑ATGGTCTGGTTg                                                                                                                                      >  2:419225/1‑141 (MQ=255)
       aacgggcgtaattaCAAATGCTCGAAGTAGAATTGATGACGGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGcc                                                                                                                                   <  1:288494/128‑1 (MQ=255)
       aaaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGcc                                                                                                                                   >  2:213863/1‑141 (MQ=255)
       aaaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGcc                                                                                                                                   <  2:257942/141‑1 (MQ=255)
        aaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGTACTACATATGTATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCa                                                                                                                                  <  2:208390/141‑1 (MQ=255)
        aaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTAATGATACGGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCCAACTAATCAACACCTTTTAGTTCTGTGTG‑ATGGTCTGGGTgtcaa                                                                                                                                  >  2:251376/1‑137 (MQ=255)
         aCGAGGAAAAGAAAAATGCAAGAGGGACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCAt                                                                                                                                 <  2:560114/141‑1 (MQ=255)
                aaaaACAAATGCTAGAGGCACAATTTTTTACTGAGACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTGTTTATAATCAACCTAATCAACGACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAg                                                                                                                          <  1:564085/141‑1 (MQ=255)
                 aaGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACCACCTTTAATTACATTTA‑ATGGTATGGGTGGCCATAAAGCacg                                                                                                                         >  2:126067/1‑139 (MQ=255)
                  aGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGga                                                                                                                        >  1:489215/1‑141 (MQ=255)
                        aTGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑AGGGTATGGGTGGCCATAAAGCAGGAGAAGtt                                                                                                                  <  2:602739/141‑1 (MQ=255)
                            tAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGCAGTTGCaa                                                                                                              >  1:11398/1‑141 (MQ=255)
                            tAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAAGCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCaa                                                                                                              >  1:201616/1‑141 (MQ=255)
                            tAGAGGCACAATTTTTTACGGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCaa                                                                                                              >  1:335907/1‑141 (MQ=255)
                            tAGAGGCACAATTTGTTACAGATACTGGACAACATAGAGATAAGCATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCaa                                                                                                              >  1:335721/1‑141 (MQ=255)
                               agGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTa                                                                                                           <  2:381393/141‑1 (MQ=255)
                                      aTTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTgg                                                                                                    >  1:508322/1‑140 (MQ=255)
                                          tttACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACa                                                                                                >  1:540952/1‑141 (MQ=255)
                                          tttACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTCTTGTGAGCATGG‑AGTGGCGGAAATAAACCAGGAGAAGTTGCAAGTAAATTTGTTACa                                                                                                >  1:138541/1‑141 (MQ=255)
                                          tttACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTGTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGCGAAGTTGCAAGTAAATTTGTTACa                                                                                                >  1:387365/1‑141 (MQ=255)
                                           ttACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAg                                                                                               <  1:343061/141‑1 (MQ=255)
                                             aCTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGAt                                                                                             <  2:284696/141‑1 (MQ=255)
                                                    cTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCTAGTAAATTTCTTACAGATGAGTTGa                                                                                      <  1:307788/141‑1 (MQ=255)
                                                    cTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                                                                                      <  1:482276/141‑1 (MQ=255)
                                                          aacaTAGAGGTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTGGTTCTGTGTG‑ATGGTATCGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATccg                                                                                >  2:403740/1‑140 (MQ=255)
                                                          aacaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATccc                                                                                <  2:219975/141‑1 (MQ=255)
                                                          aacaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATccc                                                                                >  1:436638/1‑141 (MQ=255)
                                                            caTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGt                                                                              >  1:218369/1‑141 (MQ=255)
                                                             aTAGAGATAAGAATGAGGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGtt                                                                             >  1:304925/1‑141 (MQ=255)
                                                                    tAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGc                                                                      <  1:213863/141‑1 (MQ=255)
                                                                    tAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAGAtacgctggaggacg                                                                      >  2:426243/1‑128 (MQ=255)
                                                                      aGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCgg                                                                    >  2:50281/1‑141 (MQ=255)
                                                                      aGAATGAAGATGCGGGTGGTAGTTTTTATAATCAAACTAATCAACAACTTTTAGTTATGTGTG‑ATGGTCTGGGTCGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGCAGCgg                                                                    <  1:450231/141‑1 (MQ=255)
                                                                         aTGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGGGTTGAAATCCCGTTTTGACGCGGaaa                                                                 >  1:260771/1‑141 (MQ=255)
                                                                         aTGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGaaa                                                                 >  1:504331/1‑141 (MQ=255)
                                                                          tggagagtcGGGTTGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGaaaa                                                                <  2:471251/133‑1 (MQ=255)
                                                                           gAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCACCAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAAt                                                               <  2:64568/141‑1 (MQ=255)
                                                                              gATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCtt                                                            >  1:479998/1‑141 (MQ=255)
                                                                              gATGCGGGTGATATTTTTTATAATCATACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTAGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAATTTTGTTACAGATGAGTTGAAGTCCGGTTTTGAAACGGAAAATcat                                                            >  2:374415/1‑139 (MQ=255)
                                                                               aTGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTta                                                           <  2:99368/141‑1 (MQ=255)
                                                                                 gCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTtata                                                         >  2:244276/1‑141 (MQ=255)
                                                                                   gggTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGGTGAGTTGAAATCCCGTGTTGAAGCGGAAAATCTTATAGa                                                       >  2:420483/1‑141 (MQ=255)
                                                                                       ggtATTTTTTGTAATCCACCTAATCAACACCTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaacaa                                                   <  2:201616/141‑1 (MQ=255)
                                                                                        gtATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaacaac                                                  >  1:154749/1‑141 (MQ=255)
                                                                                            tttttATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCa                                              <  1:280796/141‑1 (MQ=255)
                                                                                               ttATAATCAAACTAATCAACAACTTTTAGTTATGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGc                                           >  2:325933/1‑141 (MQ=255)
                                                                                                 acaaccacaccaaTCAACAACTGTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAACCACCAGAAGCTCCAAGTAAATTTTTTCCAGATGAGTTGAAATCCCGCTTTGAATCGGAAAATCTTATAGAACAACATCCAGCtc                                         <  2:185600/130‑2 (MQ=255)
                                                                                                  taATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGa                                        >  1:89913/1‑141 (MQ=255)
                                                                                                  taATCAAACTAATCAAAAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGa                                        <  2:129013/141‑1 (MQ=255)
                                                                                                      cAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAAt                                    >  1:258064/1‑141 (MQ=255)
                                                                                                           tAATCAACAACTTTTAGTTTCGTGTG‑ATGGTATCGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAACTGGtt                               <  2:342559/141‑1 (MQ=255)
                                                                                                                  caacTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGtaat                        <  1:50281/141‑1 (MQ=255)
                                                                                                                    acTTTTAGTTCTGTGTG‑ATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGtaataa                      <  1:217259/141‑1 (MQ=255)
                                                                                                                     cTTTTAGATCTGTGTG‑ATGGCCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTaataat                     >  1:385094/1‑141 (MQ=255)
                                                                                                                       tttAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAatat                   >  1:540520/1‑141 (MQ=255)
                                                                                                                                gtgtg‑ATGGTATGGGTGGCCATAAAGCAGGAGATGTTGCAAGTAAATTTGATACAGATGAGTTGAACTCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAatatat          <  1:416106/141‑1 (MQ=255)
                                                                                                                                    g‑ATGGTCTGGGTGGCCATAAAGAAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATATAAAt      >  2:492974/1‑141 (MQ=255)
                                                                                                                                         gTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATAAAAATTTTc  >  2:104620/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTG‑ATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATATAAATTTTC  >  CP000730/1231332‑1231608

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