Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,624,187 G→C 14.0% intergenic (‑92/+81) PP_3194 ← / ← PP_3195 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,624,1870GC14.0% 115.5 / 7.7 42intergenic (‑92/+81)PP_3194/PP_3195hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/16);  new base C (3/3);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTT  >  NC_002947/3624052‑3624336
                                                                                                                                       |                                                                                                                                                     
gTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGgc                                                                                                                                         <  1:99133/150‑1 (MQ=255)
gTGCCGTCGAACAGGGTGGAGCAGGGCAGGTGCAGGGCACTCATTGGTGAGTTCTTTCATCTTGGGGGTGGTGCAAGCCTAGGCTGTGGCGACTGCACTGTGCCTGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTGgc                                                                                                                                         <  1:322931/150‑1 (MQ=255)
    cGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGctc                                                                                                                                     >  1:539786/1‑150 (MQ=255)
        gAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTg                                                                                                                                 >  2:368228/1‑150 (MQ=255)
              ggTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCt                                                                                                                                            >  1:145022/1‑133 (MQ=255)
              ggTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCt                                                                                                                                            <  2:145022/133‑1 (MQ=255)
                 ggAGCGGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACGAGCAGATGCTATTGATGGTCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGaaa                                                                                                                        <  2:129332/150‑1 (MQ=255)
                          caggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACgg                                                                                                               >  1:52221/1‑150 (MQ=255)
                                   ggCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGTTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCg                                                                                                      <  2:218039/150‑1 (MQ=255)
                                      aCTCTTTGTTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGgcg                                                                                                   >  2:468874/1‑150 (MQ=255)
                                          aTTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGcgccgc                                                                                               >  2:270326/1‑150 (MQ=255)
                                                   ctctctCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAg                                                                                      <  1:405790/150‑1 (MQ=255)
                                                   ctctctCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAg                                                                                      <  1:48075/150‑1 (MQ=255)
                                                     ctctcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGTGACTCGGAGCCAACccagccagca                                                                                    <  2:330525/150‑3 (MQ=255)
                                                      tctcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGc                                                                                   >  1:365218/1‑150 (MQ=255)
                                                      tctcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGc                                                                                   >  1:547618/1‑150 (MQ=255)
                                                      tctcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGc                                                                                   >  1:423251/1‑150 (MQ=255)
                                                      tctcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGc                                                                                   >  1:448056/1‑150 (MQ=255)
                                                                        gCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCaa                                                                                           <  1:219051/124‑1 (MQ=255)
                                                                        gCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCaa                                                                                           >  2:219051/1‑124 (MQ=255)
                                                                            ggctaggctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAAAGAGAACGGCGCCAATCGGCGCAGCCCAAGCCAGGGCCGGCAGCGCCGATGAGGCCggg                                                             >  2:269262/1‑148 (MQ=255)
                                                                             gctaggctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCgtgt                                                            <  1:269262/150‑1 (MQ=255)
                                                                                        gCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAgg                                                 <  2:448056/150‑1 (MQ=255)
                                                                                             tGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTc                                            <  1:85944/150‑1 (MQ=255)
                                                                                             tGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTc                                            >  1:379806/1‑150 (MQ=255)
                                                                                             tGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTGgc                                                                                                                                         >  1:80499/1‑57 (MQ=255)
                                                                                             tGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTGgc                                                                                                                                         <  2:80499/57‑1 (MQ=255)
                                                                                               cccGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCgg                                          >  2:199687/1‑150 (MQ=255)
                                                                                                ccGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCAGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGc                                         >  1:161865/1‑150 (MQ=255)
                                                                                                 cGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGgc                                                                                                    <  1:334405/90‑1 (MQ=255)
                                                                                                 cGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGgc                                                                                                    >  2:334405/1‑90 (MQ=255)
                                                                                                         gcagATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTc                                <  2:365218/150‑1 (MQ=255)
                                                                                                               gCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACga                                                                                                                     >  1:336414/1‑59 (MQ=255)
                                                                                                               gCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACga                                                                                                                     <  2:336414/59‑1 (MQ=255)
                                                                                                                    tGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGt                     <  1:270326/150‑1 (MQ=255)
                                                                                                                      aTGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGc                   >  2:205308/1‑150 (MQ=255)
                                                                                                                          ccAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCgg               >  1:118491/1‑150 (MQ=255)
                                                                                                                           cAGGTATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGc              >  2:302945/1‑150 (MQ=255)
                                                                                                                              gTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAtt           >  1:189226/1‑150 (MQ=255)
                                                                                                                               tATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTa                       >  2:409948/1‑137 (MQ=255)
                                                                                                                               tATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTa                       <  1:409948/137‑1 (MQ=255)
                                                                                                                                 tCGATCGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGAt        <  2:342584/150‑1 (MQ=255)
                                                                                                                                      ggCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGc                                         >  2:364239/1‑112 (MQ=255)
                                                                                                                                      ggCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGc                                         <  1:364239/112‑1 (MQ=255)
                                                                                                                                       gCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCCATACGGTCTGCGACTGGTTGCCCGGCATTGATGCTGtt  >  2:162005/1‑150 (MQ=255)
                                                                                                                                       |                                                                                                                                                     
GTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTGTT  >  NC_002947/3624052‑3624336

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: