Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,472,112 A→G 7.5% A10A (GCA→GCG)  PP_3962 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,472,1120AG 7.5% 148.8 / 3.0 53A10A (GCA→GCG) PP_3962hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/23);  new base G (2/2);  total (28/25)
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.47e-01

CTAGCTAGAGCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGATGTGCAGATGGAAGCGCATCACTACTGCTCCATC  >  NC_002947/4471981‑4472253
                                                                                                                                   |                                                                                                                                             
ctagctagAGCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTgg                                                                                                                               >  1:32646/1‑148 (MQ=255)
    ctagAGCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGACGcc                                                                                                                          <  2:72848/149‑1 (MQ=255)
     tagAGCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCaa                                                                                                                                              >  1:171312/1‑128 (MQ=255)
     tagAGCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCaa                                                                                                                                              <  2:171312/128‑1 (MQ=255)
     tagAGCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCg                                                                                                                         >  5:215728/1‑149 (MQ=255)
       gagCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCg                                                                                                                         <  5:71615/147‑1 (MQ=255)
       gagCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCACGAAACCGCAAAGACTCAGTTGATGGAAGCCGgc                                                                                                                       >  4:24709/1‑149 (MQ=255)
         gCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGgcgc                                                                                                                     >  7:126604/1‑149 (MQ=255)
            aCTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGcgccgc                                                                                                                  >  8:75222/1‑149 (MQ=255)
                tcATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCTGTTGCTGGAAGCCGGCGCCGCGAcc                                                                                                              >  3:95668/1‑149 (MQ=255)
                  aTGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCtt                                                                                                            <  6:215728/149‑1 (MQ=255)
                              atcaGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTg                                                                                                           >  3:248571/1‑138 (MQ=255)
                              atcaGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTg                                                                                                           <  4:248571/138‑1 (MQ=255)
                               tcaGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGATCTGCTCCTGAGAAGAAGTCGCATGTGGAGACGCATGGCAGGCCAGAAACCGCAAAGACTCAGTTGATGGAAGCCGGCGCCGCGACATTGCAGGCCAAgccg                                                                                               >  3:246260/1‑149 (MQ=255)
                                  gATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAgccgcc                                                                                             >  2:23062/1‑148 (MQ=255)
                                      tataACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCAtt                                                                                         <  8:70205/148‑1 (MQ=255)
                                       ataACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCGAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCAtt                                                                                         >  2:257620/1‑147 (MQ=255)
                                       ataACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCGAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCAtt                                                                                         >  2:260930/1‑147 (MQ=255)
                                          aCCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCGAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAgg                                                                                     <  1:260930/148‑1 (MQ=255)
                                          aCCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCGAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAgg                                                                                     <  1:257620/148‑1 (MQ=255)
                                            ccTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCa                                                                                  >  1:297090/1‑149 (MQ=255)
                                                ccccTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATc                                                                                <  8:111154/147‑1 (MQ=255)
                                                      gCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGcc                                                                          <  8:126604/147‑1 (MQ=255)
                                                        cTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTAcc                                                                      >  4:157343/1‑149 (MQ=255)
                                                         tGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCt                                                                     >  6:280424/1‑149 (MQ=255)
                                                           tAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCt                                                                     <  2:32646/147‑1 (MQ=255)
                                                             gTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGAt                                                                 >  8:35326/1‑149 (MQ=255)
                                                                ttctGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGc                                                              >  8:168651/1‑149 (MQ=255)
                                                                     tgGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCTGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTc                                                           <  4:95668/147‑1 (MQ=255)
                                                                     tgGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTc                                                           <  3:258230/147‑1 (MQ=255)
                                                                     tgGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTc                                                           <  1:23062/147‑1 (MQ=255)
                                                                     tgGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCa                                                         >  5:156834/1‑149 (MQ=255)
                                                                         cTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACttt                                                     <  7:168651/149‑1 (MQ=255)
                                                                            gctctgctcCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTa                                                   >  6:107765/1‑148 (MQ=255)
                                                                              tctgctcCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACaa                                                >  5:1432/1‑149 (MQ=255)
                                                                                     cTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACAt                                         >  2:742/1‑149 (MQ=255)
                                                                                      tGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAAACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAACCCCCCCATTGAGGCCAACAACCCCAACCTGGATGGATTCCACTTTTACAACAGACaa                                         >  8:46787/1‑147 (MQ=255)
                                                                                        agaCGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATcc                                       <  8:106418/148‑1 (MQ=255)
                                                                                        agaCGAAGTCGCCTGAGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATcc                                       <  8:102096/148‑1 (MQ=255)
                                                                                          aCGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGa                                    >  2:212390/1‑149 (MQ=255)
                                                                                          aCGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGa                                    >  3:193092/1‑149 (MQ=255)
                                                                                           cGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGa                                    <  3:24709/148‑1 (MQ=255)
                                                                                            gAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGAtg                                  >  7:50303/1‑149 (MQ=255)
                                                                                            gAAGTCGCCTGTGGAGACGCCTGGCAGGACCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGAtg                                  <  7:46787/149‑1 (MQ=255)
                                                                                                tCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGAtgtg                                >  2:290770/1‑147 (MQ=255)
                                                                                                   ccTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGATGTGCAGAt                           >  4:74435/1‑149 (MQ=255)
                                                                                                    cTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCAc                                                        <  1:9934/119‑1 (MQ=255)
                                                                                                    cTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCAc                                                        >  2:9934/1‑119 (MQ=255)
                                                                                                    cTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGATGTGCAGATg                          >  5:114837/1‑149 (MQ=255)
                                                                                                         gagaCGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGTTGGCTTCCACTTTTACAACGGACATCCTGATGTGCAGATGGAAgc                     <  2:167648/149‑1 (MQ=255)
                                                                                                               ccTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGATGTGCAGATGGAAGCGCATCa               <  1:212390/149‑1 (MQ=255)
                                                                                                                   gCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGATGTGCAGATGGAAGCGCATCactac           <  7:75222/149‑1 (MQ=255)
                                                                                                                    cAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCTCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGTCCATCAACGCCTAACTGGATGGCTTCCTCTTTTACAACGGACATCCTGATGTGCAGATGGAAGCGCATCactact          <  4:246260/149‑1 (MQ=255)
                                                                                                                         ccGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGATGTGCAGATGGAAGCGCATCACTACTGCTcc     <  1:742/149‑1 (MQ=255)
                                                                                                                            aaaCCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGATGTGCAGATGGAAGCGCATCACTACTGCTCCATc  <  2:297090/149‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                             
CTAGCTAGAGCCACTCTCATGCCACAACTCATCAGATTTATAACCCTTCCCCTGGCCTGTAGTCTTCTGTGGGCTTGCTCTGCTCCTGAGACGAAGTCGCCTGTGGAGACGCCTGGCAGGCCCGAAACCGCAAAGACTCAGTTGCTGGAAGCCGGCGCCGCGACCTTGCAGGCCAAGCCGCCCATTGAGGCCATCAACGCCTACCTGGATGGCTTCCACTTTTACAACGGACATCCTGATGTGCAGATGGAAGCGCATCACTACTGCTCCATC  >  NC_002947/4471981‑4472253

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: ‑