Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 871,871 A→G 9.4% intergenic (+17/+19) PP_0755 → / ← PP_0756 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947871,8710AG 9.4% 140.3 / 5.6 53intergenic (+17/+19)PP_0755/PP_0756hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base G (2/3);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

TACATGGCAGTAGACGAGGGGGTAGTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTT  >  NC_002947/871729‑872018
                                                                                                                                              |                                                                                                                                                   
tACATGGCAGTAGACGAGGGGGTAGTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAg                                                                                                                                               <  3:17428/149‑1 (MQ=255)
     ggCAGTAGACGAGGGGGTAGTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGcccc                                                                                                                                           >  1:193256/1‑148 (MQ=255)
      gCAGTAGACGAGGGGGTAGTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGccc                                                                                                                                            >  1:214930/1‑146 (MQ=255)
      gCAGTAGACGAGGGGGTAGTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGccc                                                                                                                                            <  2:214930/146‑1 (MQ=255)
       cAGTAGACGAGGGGGTAGTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATTGGGGCTGCGTTGCAGccccgat                                                                                                                                        <  4:31464/149‑4 (MQ=255)
             aCGAGGGGGTAGTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTGTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGAt                                                                                                                                  <  4:341592/149‑1 (MQ=255)
                      tAGTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTa                                                                                                                           >  3:267731/1‑147 (MQ=255)
                        gTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAgg                                                                                                                       >  3:63093/1‑149 (MQ=255)
                           gACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACAATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGc                                                                                                                    <  7:54312/149‑1 (MQ=255)
                               cTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCctgc                                                                                                                <  4:63093/149‑1 (MQ=255)
                                 ttCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCAGCAGGCACctgtag                                                                                                              >  1:306930/1‑146 (MQ=255)
                                       cTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGcctgac                                                                                                        >  5:300962/1‑149 (MQ=255)
                                           gggTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCAGCAGGCACCTGTAGcctgatcag                                                                                                     <  2:306930/148‑5 (MQ=255)
                                            ggTGTCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGcctgacctga                                                                                                    >  2:156616/1‑148 (MQ=255)
                                                  gAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACGTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGt                                                                                             >  8:5757/1‑149 (MQ=255)
                                                  gAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCGTTGCAGccccgat                                                                                                                                        >  6:291391/1‑103 (MQ=255)
                                                  gAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCGTTGCAGccccgat                                                                                                                                        <  5:291391/106‑4 (MQ=255)
                                                    aaTACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTg                                                                                           <  6:300962/149‑1 (MQ=255)
                                                    aaTACTTCTACATCCTGGGCTGTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTg                                                                                           <  4:24211/149‑1 (MQ=255)
                                                      gacgtCTACATGAAGGGCTTTGCGTGCGTGGTGATGGCGCTGCGGCTGTCCTGGTTCAGCAGCTGGGTCTGAGAACCTTCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGTCCTGATACTGGTTgcg                                                                                         <  8:161808/145‑1 (MQ=255)
                                                       aCTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTg                                                                                               <  8:319246/142‑1 (MQ=255)
                                                       aCTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTg                                                                                               >  7:319246/1‑142 (MQ=255)
                                                        cTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGccc                                                                                                                                            <  5:298198/96‑1 (MQ=255)
                                                        cTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGccc                                                                                                                                            >  6:298198/1‑96 (MQ=255)
                                                                  cTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCg                                                                              >  1:13167/1‑148 (MQ=255)
                                                                   tGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAAc                                                                                                                                       <  3:313485/90‑1 (MQ=255)
                                                                   tGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAAc                                                                                                                                       >  4:313485/1‑90 (MQ=255)
                                                                    gggCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCa                                                                            >  3:103742/1‑148 (MQ=255)
                                                                            cGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGcctg                                                                                                          <  7:236255/110‑1 (MQ=255)
                                                                            cGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGcctg                                                                                                          >  8:236255/1‑110 (MQ=255)
                                                                               gCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAacac                                                                 >  2:332939/1‑148 (MQ=255)
                                                                                 cTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCaa                                                                                                                                                    >  4:316525/1‑63 (MQ=255)
                                                                                 cTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCaa                                                                                                                                                    <  3:316525/63‑1 (MQ=255)
                                                                                    gCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGt                                                            <  1:332939/148‑1 (MQ=255)
                                                                                      gATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACg                                                         <  2:193256/149‑1 (MQ=255)
                                                                                        tGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGc                                                       <  3:346229/149‑1 (MQ=255)
                                                                                         ggCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGcc                                                      >  7:333673/1‑149 (MQ=255)
                                                                                            gctgctgctGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCgct                                                   <  4:267731/149‑1 (MQ=255)
                                                                                                ctgctgTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGctc                                               <  1:156616/149‑1 (MQ=255)
                                                                                                  gctgTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCtt                                             >  4:121265/1‑149 (MQ=255)
                                                                                                   ctgTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGt                                                                                             <  7:75965/100‑1 (MQ=255)
                                                                                                   ctgTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGt                                                                                             >  8:75965/1‑100 (MQ=255)
                                                                                                               ccccgatGGTTCTCAGACCCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGAAACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATc                                >  4:286030/7‑149 (MQ=255)
                                                                                                                    cTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGgc                           <  3:121265/149‑1 (MQ=255)
                                                                                                                        gTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCg                       >  7:199115/1‑149 (MQ=255)
                                                                                                                        gTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCACGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCg                       >  4:86873/1‑149 (MQ=255)
                                                                                                                         tCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGt                      <  4:103742/149‑1 (MQ=255)
                                                                                                                                  cATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCAttct             >  8:148633/1‑149 (MQ=255)
                                                                                                                                  cATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCAttct             >  3:273563/1‑149 (MQ=255)
                                                                                                                                   aTCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAg                                                              <  2:204772/99‑1 (MQ=255)
                                                                                                                                    ttggggCTGCGTTGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGggca                                      <  4:298659/120‑1 (MQ=255)
                                                                                                                                    ttggggCTGCGTTGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGggca                                      >  3:298659/3‑122 (MQ=255)
                                                                                                                                      ggggCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTcc          <  7:5757/148‑1 (MQ=255)
                                                                                                                                             cAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTccttggcctt  >  6:181474/1‑149 (MQ=255)
                                                                                                                                             cAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTccttggcc    <  7:173958/147‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TACATGGCAGTAGACGAGGGGGTAGTGGACGCTTTCGAACTTCGGGTGCCGAAATACTTCTACATCCTGGGCTTTTCGTGCCTGGCGATGGCGCTGCTGCTGTCCTGGTTCAGCACCTGGGTCTGAGAACCATCGGGGCTGCAACGCAGCCCCAAACCTGATCAGGCTACAGGTGCCTGCTGCCTGACCTGATACTGGTTGCGCACCGGCGTCGCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTT  >  NC_002947/871729‑872018

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: