Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 104,876 G→C 10.5% P345P (CCG→CCC)  PP_0099 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947104,8760GC10.5% 166.4 / 8.5 57P345P (CCG→CCC) PP_0099hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/28);  new base C (4/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GTATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCT  >  NC_002947/104732‑105020
                                                                                                                                                |                                                                                                                                                
gTACGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCTGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTg                                                                                                                                              >  5:37385/1‑149 (MQ=255)
 tATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTgg                                                                                                                                             <  6:285729/149‑1 (MQ=255)
   tGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTggg                                                                                                                                            <  3:279909/148‑1 (MQ=255)
    gACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGTGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTgg                                                                                                                                             <  8:254040/146‑1 (MQ=255)
    gACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTgg                                                                                                                                             >  7:254040/1‑146 (MQ=255)
          cGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGgc                                                                                                                                    <  1:41912/149‑1 (MQ=255)
          cGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGgc                                                                                                                                    <  2:127015/149‑1 (MQ=255)
            gTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCCGGTGGGCaaaacgc                                                                                                                                    <  3:241667/147‑7 (MQ=255)
            gTTCGGCGTGCGCATGCGTGGTCAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGCCCGCGAGGCGCACGAGCTGGATTGCTCGGCGATTTGACCACCCGGTGGGCaaaacgcc                                                                                                                                   >  4:241667/1‑141 (MQ=255)
                    tGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTCCCCACCGGGTGGGCAGATGGCGCTGTTCTGa                                                                                                                          >  7:109461/1‑149 (MQ=255)
                         aTGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGa                                                                                                                     <  3:191543/149‑1 (MQ=255)
                          tGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGaa                                                                                                                    >  1:338803/1‑149 (MQ=255)
                            cGTGGTGAAGGAGTATTTGCTGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAac                                                                                                                   <  6:37385/148‑1 (MQ=255)
                             gtggtgAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAaccc                                                                                                                 <  3:40452/149‑1 (MQ=255)
                             gtggtgAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAaccc                                                                                                                 >  4:40452/1‑149 (MQ=255)
                                  gAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAaccccaccc                                                                                                            >  3:290620/1‑149 (MQ=255)
                                      gAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGCACCCCACCCCCca                                                                                                        >  1:123231/1‑149 (MQ=255)
                                          aTTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATgg                                                                                                                                     <  6:88445/116‑1 (MQ=255)
                                          aTTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATgg                                                                                                                                     >  5:88445/1‑116 (MQ=255)
                                               ccGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCt                                                                                               <  7:256538/149‑1 (MQ=255)
                                                    ttGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCa                                                                                          <  4:290620/149‑1 (MQ=255)
                                                       cTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCAtt                                                                                        <  5:214144/148‑1 (MQ=255)
                                                           cgcAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTca                                                                                   <  3:327394/149‑1 (MQ=255)
                                                            gcAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGTGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTcac                                                                                  <  2:57740/149‑1 (MQ=255)
                                                              aGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTCTGCCCACCCGGTGGGCaaaac                                                                                                                                      >  7:200713/1‑91 (MQ=255)
                                                              aGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTCTGCCCACCCGGTGGGCaaaac                                                                                                                                      <  8:200713/95‑5 (MQ=255)
                                                               gcgcTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAAc                                                                               >  4:238504/1‑149 (MQ=255)
                                                                       ccAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGAc                                                                                                                         >  3:143198/1‑99 (MQ=255)
                                                                       ccAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGAc                                                                                                                         <  4:143198/99‑1 (MQ=255)
                                                                       ccAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACaa                                                                                >  6:186161/1‑140 (MQ=255)
                                                                       ccAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACaa                                                                                <  5:186161/140‑1 (MQ=255)
                                                                        cAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCTGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaa                                                                       >  7:87187/1‑148 (MQ=255)
                                                                        cAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCACCCACAGGGGGCTCCGCATTCCACACAACGCAACtaa                                                                       >  6:55711/1‑148 (MQ=255)
                                                                         aaGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACACCACCCACCACAGGGGGCTGCGCATTGCTCACAACGCAACtaat                                                                      >  7:259967/1‑148 (MQ=255)
                                                                          aGGCAATGAACCGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaataa                                                                    <  8:305466/149‑1 (MQ=255)
                                                                           ggCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATaac                                                                   <  8:109461/149‑1 (MQ=255)
                                                                             cAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATaacaa                                                                 <  2:241181/149‑1 (MQ=255)
                                                                                       gCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTg                                                       <  2:338803/149‑1 (MQ=255)
                                                                                        cTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgc                                                      <  3:238504/149‑1 (MQ=255)
                                                                                         tGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgcg                                                     >  4:251023/1‑149 (MQ=255)
                                                                                             tttCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTg                                                                                       <  7:10354/111‑1 (MQ=255)
                                                                                             tttCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTg                                                                                       >  8:10354/1‑111 (MQ=255)
                                                                                               tCGAAGGGCGCTATGCGCAGGAGCTGGATTGCTCGGCGTTTTTCCCACCGGGTTGTCAGATGGCGCTGTTCTGACATGAACACCACCCCCCACAGGGGTCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAAt                                                <  8:259967/148‑1 (MQ=255)
                                                                                                 gAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACa                                             >  1:141372/1‑149 (MQ=255)
                                                                                                      gcgcGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACt                                                                         <  8:133359/116‑1 (MQ=255)
                                                                                                      gcgcGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACt                                                                         >  7:133359/1‑116 (MQ=255)
                                                                                                          gAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCaa                                    <  3:251023/149‑1 (MQ=255)
                                                                                                            ggCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGgcg                                                                                                                                   >  5:136241/1‑52 (MQ=255)
                                                                                                            ggCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGgcg                                                                                                                                   <  6:136241/52‑1 (MQ=255)
                                                                                                               gcAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGAt                               >  3:154782/1‑149 (MQ=255)
                                                                                                                      cTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGtttttt                        <  2:123231/149‑1 (MQ=255)
                                                                                                                         ttttGCTCGGCGTTTTGCCCTCCGGGTTGTCAGATTGAGATTTTATGAAATTAACCCCACCCCCCACAGGGGGCTGCGAATTTCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATAATGGCAACTGATGtttttt                        <  7:120682/144‑1 (MQ=255)
                                                                                                                                 ggCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCCCAGGGGGCTGCGCATTGCTGACAACGCAACTAATAACAACCAATAATTGCGGTATTTCATATTGGCTACTGATGTTTTTTTGCTTTTATcg             >  4:232113/1‑148 (MQ=255)
                                                                                                                                 ggCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCtattat               >  1:128391/1‑147 (MQ=255)
                                                                                                                                   ccatctGCCCACCCGGTGGGCAAATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAg           >  7:85690/6‑149 (MQ=255)
                                                                                                                                   ccatctGCCCACCCGGTGGGCAAATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAg           >  1:213230/6‑149 (MQ=255)
                                                                                                                                            cACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCt  >  4:251175/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GTATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCT  >  NC_002947/104732‑105020

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