Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,224,278 G→C 7.9% G104G (GGC→GGG)  PP_1962 ← site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,224,2780GC 7.9% 148.8 / 4.4 52G104G (GGC→GGG) PP_1962site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base C (2/2);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

AGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGA  >  NC_002947/2224129‑2224425
                                                                                                                                                     |                                                                                                                                                   
aGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTgggg                                                                                                                                                     <  1:1487655/150‑1 (MQ=255)
    aTAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCa                                                                                                                                                 >  1:1751552/1‑150 (MQ=255)
      aaGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAAt                                                                                                                                               <  2:851631/150‑1 (MQ=255)
        ggAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAg                                                                                                                                             >  1:1793307/1‑150 (MQ=255)
        ggAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATATAAATCTAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAg                                                                                                                                             >  1:116674/1‑150 (MQ=255)
         gAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAgg                                                                                                                                            >  1:843985/1‑150 (MQ=255)
         gAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGACCCAATAgg                                                                                                                                            >  2:1735322/1‑150 (MQ=255)
         gAAAGCCAAAATCGCAACATGAAAAACAGGCATGAACAATGACTCGTTTTGCCACTTGCTTTTGCGTCCTTGAGGGCCAGCTACAAGGGCGGAGTGGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCACATTCGGGCCCAATAgg                                                                                                                                            >  1:611494/1‑150 (MQ=255)
          aaaGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGt                                                                                                                                           <  2:824866/150‑1 (MQ=255)
                        aaCATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAAc                                                                                                                             <  2:843985/150‑1 (MQ=255)
                        aaCATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAAc                                                                                                                             >  1:1923425/1‑150 (MQ=255)
                            tGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGAt                                                                                                                         >  1:463557/1‑150 (MQ=255)
                                  aCAGTCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGt                                                                                                                   <  1:26022/150‑1 (MQ=255)
                                        cTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGAt                                                                                                                         <  2:362662/138‑1 (MQ=255)
                                        cTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGAt                                                                                                                         >  1:362662/1‑138 (MQ=255)
                                        cTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGAtt                                                                                                             >  1:1209016/1‑150 (MQ=255)
                                         tGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGAttt                                                                                                            >  2:1905424/1‑150 (MQ=255)
                                               atGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCttt                                                                                                      <  2:1789974/150‑1 (MQ=255)
                                               atGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCttt                                                                                                      >  2:785384/1‑150 (MQ=255)
                                                   cTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGgat                                                                                                  <  2:615738/150‑1 (MQ=255)
                                                      ggTTGTCCACTTTCTTTTGCGTCGTCGTTCCCATCATTACTCTTCCGATTTCTTATCAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTTAGATTTCGCTTTGgatgat                                                                                               <  2:116674/150‑1 (MQ=255)
                                                                    ttttGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAAt                                                                                 >  2:1458716/1‑150 (MQ=255)
                                                                      ttGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATAc                                                                               <  2:1923425/150‑1 (MQ=255)
                                                                              tcgAGGGCCCGATTAGGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGGACGAATACGATTGAAg                                                                       >  2:80963/1‑150 (MQ=255)
                                                                                  gggCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAg                                                                   <  2:482168/150‑1 (MQ=255)
                                                                                      ccGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTaaa                                                               <  2:1793307/150‑1 (MQ=255)
                                                                                       cGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGATGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAAt                                                              <  1:115020/150‑1 (MQ=255)
                                                                                           tAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAg                                                          <  2:463557/150‑1 (MQ=255)
                                                                                           tAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAg                                                          <  1:1905424/150‑1 (MQ=255)
                                                                                                 tCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGaa                                                                                          >  2:740627/1‑112 (MQ=255)
                                                                                                 tCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGaa                                                                                          <  1:740627/112‑1 (MQ=255)
                                                                                                          ttATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTaaa                                           >  1:1346844/1‑150 (MQ=255)
                                                                                                          ttATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTaaa                                           >  2:1810945/1‑150 (MQ=255)
                                                                                                             tGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACg                                                                                            <  1:1607613/98‑1 (MQ=255)
                                                                                                             tGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACg                                                                                            >  2:1607613/1‑98 (MQ=255)
                                                                                                              gagaTCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAAt                                       >  1:257270/1‑150 (MQ=255)
                                                                                                                gaTCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTg                                     >  1:1098300/1‑150 (MQ=255)
                                                                                                                  tCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGcc                                   <  1:1458716/150‑1 (MQ=255)
                                                                                                                          cTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAg                           >  1:1466154/1‑150 (MQ=255)
                                                                                                                          cTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAg                           <  2:1346844/150‑1 (MQ=255)
                                                                                                                           tCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAgg                          >  2:1693221/1‑150 (MQ=255)
                                                                                                                                aaTCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGtt                     <  1:785384/150‑1 (MQ=255)
                                                                                                                                          aGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGt           <  2:1466154/150‑1 (MQ=255)
                                                                                                                                           gCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTa          >  1:1771205/1‑150 (MQ=255)
                                                                                                                                            ccTATTGGGCCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTaaa                                           <  2:22125/116‑1 (MQ=255)
                                                                                                                                            ccTATTGGGCCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTaaa                                           >  1:22125/1‑116 (MQ=255)
                                                                                                                                             cTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTaaa                                           >  1:403431/1‑115 (MQ=255)
                                                                                                                                             cTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTaaa                                           <  2:403431/115‑1 (MQ=255)
                                                                                                                                             cTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTAtc        >  2:1879426/1‑150 (MQ=255)
                                                                                                                                              tATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAAGGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTAtct       >  2:271633/1‑150 (MQ=255)
                                                                                                                                               aTTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTAtctc      >  1:1428086/1‑150 (MQ=255)
                                                                                                                                               aTTGGGCCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTc                    <  2:518934/136‑1 (MQ=255)
                                                                                                                                               aTTGGGCCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTc                    >  1:518934/1‑136 (MQ=255)
                                                                                                                                                   gggCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGa  <  2:320867/150‑1 (MQ=255)
                                                                                                                                                   gggCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGa  <  2:1098300/150‑1 (MQ=255)
                                                                                                                                                    ggCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGa                                                       >  2:1421409/1‑96 (MQ=255)
                                                                                                                                                    ggCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGa                                                       <  1:1421409/96‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGA  >  NC_002947/2224129‑2224425

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