Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,546,869 G→A 10.9% intergenic (+25/+15) PP_2238 → / ← PP_5504 metallopeptidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,546,8690GA10.9% 137.8 / 7.8 55intergenic (+25/+15)PP_2238/PP_5504metallopeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base A (2/4);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAG  >  NC_002947/2546727‑2547010
                                                                                                                                              |                                                                                                                                             
gACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGcc                                                                                                                                        <  1:645122/150‑1 (MQ=255)
   cAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtt                                                                                                                                     >  1:1232463/1‑150 (MQ=255)
    aaGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCttt                                                                                                                                    <  1:1921864/150‑1 (MQ=255)
      ggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGGTGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtc                                                                                                                                  >  1:1648600/1‑150 (MQ=255)
      ggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtc                                                                                                                                  >  1:454487/1‑150 (MQ=255)
       gaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGACCTGGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtttttc                                                                                                                                 <  2:1232463/150‑3 (MQ=255)
         ggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtttttccg                                                                                                                               <  2:540162/150‑5 (MQ=255)
          gaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTa                                                                                                                              >  2:56529/1‑150 (MQ=255)
                   tGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGCTTGCGAACTGCCCTTTTCTCTATTTACAgat                                                                                                                     <  1:1601882/150‑1 (MQ=255)
                         aaCGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCCGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatgat                                                                                                               <  2:468201/150‑1 (MQ=255)
                          aCGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                                                                                              <  2:1285986/150‑1 (MQ=255)
                              cTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGc                                                                                                          >  1:1740291/1‑150 (MQ=255)
                                 ccACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTcc                                                                                                       >  2:1307820/1‑150 (MQ=255)
                                  cACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                                                                                              <  1:578336/142‑1 (MQ=255)
                                  cACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                                                                                              >  2:578336/1‑142 (MQ=255)
                                   aCGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATATGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCgg                                                                                                     >  1:452050/1‑150 (MQ=255)
                                                 gAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCccgc                                                                                       >  1:672791/1‑150 (MQ=255)
                                                    cgcCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCgg                                                                                    <  1:828577/150‑1 (MQ=255)
                                                       cTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGtt                                                                                 >  1:412771/1‑150 (MQ=255)
                                                          gcgcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGTCCTTTTCTCTATTTTCAGATTATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGag                                                                              <  1:1307820/150‑2 (MQ=255)
                                                          gcgcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGAt                                                                              <  1:56529/150‑1 (MQ=255)
                                                           cgcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATg                                                                             <  1:1117903/150‑1 (MQ=255)
                                                            gcCCTTGCTCCAAGGGGATGCGCTAGCATGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGc                                                                            >  1:1843191/1‑150 (MQ=255)
                                                             cccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTAGTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGcc                                                                           <  1:1321147/150‑1 (MQ=255)
                                                              ccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCg                                                                          >  1:275622/1‑150 (MQ=255)
                                                                ttGCTCCAACGGGATGCGCTAGCCTGGCTGGAGGTGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACNGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGaa                                                                        <  1:1317833/150‑1 (MQ=255)
                                                                 tGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAAc                                                                       <  2:869150/150‑1 (MQ=255)
                                                                   cTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACaa                                                                     >  1:791464/1‑150 (MQ=255)
                                                                      cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACTGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGGTGAACAAATa                                                                  >  1:1751053/1‑150 (MQ=255)
                                                                      cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATa                                                                  <  1:221405/150‑1 (MQ=255)
                                                                         ggggATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTc                                                               >  2:1760065/1‑150 (MQ=255)
                                                                           ggATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCTCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCAc                                                             <  2:1707793/150‑1 (MQ=255)
                                                                                        tggctggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATc                                                <  2:412771/150‑1 (MQ=255)
                                                                                        tggctggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATc                                                >  1:486211/1‑150 (MQ=255)
                                                                                            tggAGATGCGTACTGCACCGCTGTCACGCAGCGGAAAAAGGGCAGTTTGCGAAATGCGCTTTTCTCTATTTTTAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTa                                            >  2:813240/1‑150 (MQ=255)
                                                                                            tggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTa                                            >  1:1275382/1‑150 (MQ=255)
                                                                                                gaTGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGc                                        >  2:45094/1‑150 (MQ=255)
                                                                                                    cGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGa                                    <  2:486211/150‑1 (MQ=255)
                                                                                                      tACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGc                                                         >  1:699217/1‑127 (MQ=255)
                                                                                                      tACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGc                                                         <  2:699217/127‑1 (MQ=255)
                                                                                                      tACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAgc                                  >  1:1224436/1‑150 (MQ=255)
                                                                                                           cACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtct                                                                                                                                 <  1:900817/50‑1 (MQ=255)
                                                                                                           cACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtct                                                                                                                                 >  2:900817/1‑50 (MQ=255)
                                                                                                           cACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatgat                                                                                                               >  2:372575/1‑68 (MQ=255)
                                                                                                           cACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatgat                                                                                                               <  1:372575/68‑1 (MQ=255)
                                                                                                              cGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATg                          >  2:944351/1‑150 (MQ=255)
                                                                                                               gCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGa                         <  2:672791/150‑1 (MQ=255)
                                                                                                                         aGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGCGCGCGCAATGACCCGCTAGcc               >  2:43295/1‑150 (MQ=255)
                                                                                                                           ctgAAAAGGGGCCTTTTGCGAACTGCCCTTTTCTCTCTTTACAGATGATGATCATGCTCCGCCTGGTGTCGCCCGCCGGGTTGATGCCGAACAAATAGTCGCCGGCTTCTTCATCGCTATCGTCCGAGCGCGCAATGACCCGGTAGCcgc             <  1:813240/148‑1 (MQ=255)
                                                                                                                             gAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCcgcgc           <  2:1224436/150‑1 (MQ=255)
                                                                                                                             gAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCcgcgc           >  1:454574/1‑150 (MQ=255)
                                                                                                                                    ggCAGTTCGCAAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAgc                                  <  2:322487/119‑1 (MQ=255)
                                                                                                                                    ggCAGTTCGCAAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAgc                                  >  1:322487/1‑119 (MQ=255)
                                                                                                                                     gCAGTTTGCGAACTGCCCTTTTATACCCGTGGAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCa   >  1:1507427/1‑150 (MQ=255)
                                                                                                                                      cAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAg  <  1:1760065/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAG  >  NC_002947/2546727‑2547010

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