Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,624,144 G→C 18.1% 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,1440GC18.1% 99.4 / 17.2 47intergenic (‑92/+81)PP_3194/PP_3195hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/14);  new base C (5/3);  total (30/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

GGTGTGCAGGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCT  >  NC_002947/3623996‑3624290
                                                                                                                                                    |                                                                                                                                                  
ggTGTGCAGGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATgg                                                                                                                                                    <  5:183265/149‑1 (MQ=255)
              tGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTGgc                                                                                                                                      <  8:171607/149‑1 (MQ=255)
                 cGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCt                                                                                                                                         >  1:19718/1‑143 (MQ=255)
                 cGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCt                                                                                                                                         <  2:19718/143‑1 (MQ=255)
                 cGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTg                                                                                                                                        <  6:192771/144‑1 (MQ=255)
                 cGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTg                                                                                                                                        >  5:192771/1‑144 (MQ=255)
                   tCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGAGctc                                                                                                                                  >  8:82321/1‑148 (MQ=255)
                     gAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATAGATACATGGCGctctct                                                                                                                               >  6:224244/1‑149 (MQ=255)
                        cAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCAGCGATACCTGGCGCTCTCTGgc                                                                                                                            >  1:277467/1‑149 (MQ=255)
                        cAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATGGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGgc                                                                                                                            <  7:207100/149‑1 (MQ=255)
                                aGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAAc                                                                                                                    <  2:47994/149‑1 (MQ=255)
                                aGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGGGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAAc                                                                                                                    <  2:43233/149‑1 (MQ=255)
                                 gCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTgg                                                                                                                                       <  1:136638/129‑1 (MQ=255)
                                 gCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTgg                                                                                                                                       >  2:136638/1‑129 (MQ=255)
                                   aGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAAc                                                                                                                    >  3:38025/1‑146 (MQ=255)
                                     tgcaggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGaa                                                                                                               >  6:167019/1‑149 (MQ=255)
                                     tgcaggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGaa                                                                                                               >  4:2111/1‑149 (MQ=255)
                                     tgcaggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGaa                                                                                                               >  7:232879/1‑149 (MQ=255)
                                           tgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGgcgc                                                                                                         >  5:291751/1‑149 (MQ=255)
                                                   aCTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGAGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGgc                                                                                                 <  7:82321/149‑1 (MQ=255)
                                                            tGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGtcc                                                                                          >  2:7609/1‑147 (MQ=255)
                                                            tGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCaa                                                                                        <  7:123318/149‑1 (MQ=255)
                                                               gCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTAATGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGc                                                                                      >  1:98280/1‑148 (MQ=255)
                                                                ctctctCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGc                                                                                      >  1:223678/1‑147 (MQ=255)
                                                                     tcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAgg                                                                                  <  5:167019/146‑1 (MQ=255)
                                                                     tcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGcc                                                                               >  7:83812/1‑149 (MQ=255)
                                                                       aTCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCaa                                                                                        <  7:87599/138‑1 (MQ=255)
                                                                       aTCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCaa                                                                                        >  8:87599/1‑138 (MQ=255)
                                                                               ggtggtGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGcc                                                                      >  4:190604/1‑148 (MQ=255)
                                                                                    tgCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAg                                                                >  5:42069/1‑149 (MQ=255)
                                                                                     gCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGTCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAgg                                                               >  7:247324/1‑149 (MQ=255)
                                                                                     gCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGTCCAAGCCAGGGCCGGCAGCGCCGATCAgg                                                               <  1:7609/149‑1 (MQ=255)
                                                                                     gCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAgg                                                               >  3:287626/1‑149 (MQ=255)
                                                                                      cAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGc                                                              >  5:284500/1‑149 (MQ=255)
                                                                                      cAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGc                                                              >  5:9019/1‑149 (MQ=255)
                                                                                               gctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGgcc                                                     <  2:223678/149‑1 (MQ=255)
                                                                                               gctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGgcc                                                     >  8:111963/1‑149 (MQ=255)
                                                                                                                       cagATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGAGCCAATCGGCGCCGCGCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTc                             <  5:224244/149‑1 (MQ=255)
                                                                                                                           ttcttttGCTGGCCCGGTATCGAGGGCTTCGTTACCTGGCGCTCTCTGGCGCATAAACGTGACCGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGc                          <  4:165846/143‑1 (MQ=255)
                                                                                                                            gCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGc                          >  3:165846/1‑147 (MQ=255)
                                                                                                                                        cAGTTATCTATGGCATCTATACCTGGCGCTCTCTGGCGCAGAAACTATAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCgg            <  8:232879/149‑1 (MQ=255)
                                                                                                                                        cAGGTATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCgg            >  6:246149/1‑149 (MQ=255)
                                                                                                                                        cAGGTATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCgg            >  7:67786/1‑149 (MQ=255)
                                                                                                                                          ggTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCa          >  7:252056/1‑149 (MQ=255)
                                                                                                                                           gTATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAt         >  4:227339/1‑149 (MQ=255)
                                                                                                                                             aTCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAtt        >  6:245754/1‑148 (MQ=255)
                                                                                                                                             aTCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAtt        >  1:51057/1‑148 (MQ=255)
                                                                                                                                                  tGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCt  <  4:287626/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTGTGCAGGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCT  >  NC_002947/3623996‑3624290

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