Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,466,841 T→G 21.2% intergenic (‑16/+108) PP_3081 ← / ← PP_5560 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,466,8410TG21.2% 78.8 / 8.2 49intergenic (‑16/+108)PP_3081/PP_5560hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (12/26);  new base G (8/3);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 3.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTACCTGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGG  >  NC_002947/3466693‑3466986
                                                                                                                                                    |                                                                                                                                                 
gcgTACCTGAGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCt                                                                                                                                                   <  4:343064/149‑1 (MQ=255)
      cTGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGgcggcg                                                                                                                                             >  2:162657/1‑149 (MQ=255)
         cGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGcggcgg                                                                                                                                            <  6:347524/147‑1 (MQ=255)
           ggCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCCCCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGCAGGCGGCCaa                                                                                                                                        <  5:24925/149‑1 (MQ=255)
           ggCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCaa                                                                                                                                        <  7:352599/149‑1 (MQ=255)
           ggCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCaa                                                                                                                                        <  6:95020/149‑1 (MQ=255)
              aGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCGGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTTAAAGACGAGTAACTGACCGACATCGCGCCTTGGCCGCCGCAGGCGGCCAAGGc                                                                                                                                     <  6:375249/149‑1 (MQ=255)
                gTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGTAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCtg                                                                                                                                   >  6:292767/1‑149 (MQ=255)
                  tGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCtgtg                                                                                                                                 >  4:116946/1‑149 (MQ=255)
                      ccAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCCCCTGCCGCGGCCAACGCGGTGCAgg                                                                                                                             >  4:94570/1‑149 (MQ=255)
                       cAGATCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGCCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTCGCCGCCTGCCGCGCCCAACGCTGAGCAggg                                                                                                                            >  4:308706/1‑149 (MQ=255)
                         gCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGgca                                                                                                                          <  3:116946/149‑1 (MQ=255)
                          cTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGGCGGCGGCCCAAGCTGTGCAGGgccg                                                                                                                         >  6:9196/1‑147 (MQ=255)
                              tttGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCaaa                                                                                                                     <  3:94570/148‑1 (MQ=255)
                               ttGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAag                                                                                                                    <  5:292767/149‑1 (MQ=255)
                                 gCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGGCAACGCAGTGCAGGGCAGCAAagcc                                                                                                                  >  7:323713/1‑148 (MQ=255)
                                 ccGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGAGCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCggcggc                                                                                                                                           <  4:66900/123‑1 (MQ=255)
                                    gggTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGgca                                                                                                                          <  7:3/138‑1 (MQ=255)
                                    gggTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGGCGGCGGCCCCCGCGGTTCAgggga                                                                                                                          >  8:3/1‑136 (MQ=255)
                                      gTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGaa                                                                                                             <  7:66683/149‑1 (MQ=255)
                                            tctTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCCGCGAACCCGGTTCAGTGCACGCAACCCGCGAACAGc                                                                                                         >  5:329663/1‑147 (MQ=255)
                                              ttCACCGGCACCGGATGGCGGGTCTTGGAGGACCTGGTGCAGGGCATTGACGCCGACATCCGGCAAAGCCCAGGAACTGACCGCCATGGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCa                                                                                                       <  8:161983/147‑1 (MQ=255)
                                                cACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAgtgt                                                                                                   >  1:71309/1‑149 (MQ=255)
                                                           aaTGCGGGGTCCTGGCCGACCTGTGGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCAt                                                                                        <  3:237197/148‑1 (MQ=255)
                                                            aTGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATg                                                                                       <  2:71309/149‑1 (MQ=255)
                                                             tGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATCCCGCCGCCATCAGTCAACGACGAGCAACTCACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATgg                                                                                      <  2:191505/149‑1 (MQ=255)
                                                               ccGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGTAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGt                                                                                     <  3:308706/148‑1 (MQ=255)
                                                                 gggTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCCACGCTGTGCAGGGCAGCAAAGCACCGAACAGCCAGTGGCTGCGGTGCATGGTCa                                                                                   >  2:59185/1‑148 (MQ=255)
                                                                   gTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCCCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATc                                                                                <  6:220954/149‑1 (MQ=255)
                                                                         gCCGACCTGTTGCAGGGCAATGACGCCCCCATCAGCTATAGACGAGGAACTGACCGACTTCGCGCCTTGGCGGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTgg                                                                           <  2:283979/148‑1 (MQ=255)
                                                                             aCCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGa                                                                                               >  7:264359/1‑124 (MQ=255)
                                                                             aCCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGa                                                                                               <  8:264359/124‑1 (MQ=255)
                                                                             aCCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAggg                                                                      <  3:302633/149‑1 (MQ=255)
                                                                             aCCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGTCGGCGGCCAACCCTTTGCGGGGCAGCAAACCAGCGAACAGCCAGGGTCAGGGGGGCAAGGTGCAACCTTAGCAggg                                                                      >  6:329299/1‑149 (MQ=255)
                                                                                    gCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGGCGCCGGCGGGGGCCCACCCGGTGGAGGGCCGCCAAGCCGCGCAAACCCAGGGGGTGAGGTGGATGGTGCATCCTTGGCAGGGCGGCAcc                                                               >  4:173940/1‑148 (MQ=255)
                                                                                    gCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCAc                                                                <  7:93235/148‑1 (MQ=255)
                                                                                        ggCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTc                                                           <  1:162657/149‑1 (MQ=255)
                                                                                               aCGCCGCCATCAGTGAAAGACGAGCAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAg                                                    <  3:59627/149‑1 (MQ=255)
                                                                                                       aTCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCGTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTc                                            <  5:329299/149‑1 (MQ=255)
                                                                                                        tCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCt                                           >  8:288997/1‑149 (MQ=255)
                                                                                                              aaaGACGAGGAACTGACCGACATCGCGCCTTGGCCCCCGGCGGCCGGCAACCCTGTGCAGCGCAGGAAAGCAGCGCACAGCCAGAGGCGGAGGGGGATGGGGAATCCATGGGCGGGCCGCACCGCGGAGTCTAGCCCAGGTGTGCTTCa                                     >  2:66191/1‑149 (MQ=255)
                                                                                                                  aCGAGGATCTGACCGACATCGCGCCTTGGCCGCCTGCGGGGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAg                                                                                              <  8:206209/88‑1 (MQ=255)
                                                                                                                  aCGAGGAACTGACCGACATCGCGCCTTGGCCCCCTGCGGCGGCCCACGCAGTGCAGGGCCGCAAAGCACCCCACACCCAGTGTCTGAg                                                                                              >  7:206209/1‑88 (MQ=255)
                                                                                                                    gAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGGGTCTGAGGTGCATGGGGAATCCTGGGCAGGGCGGCCCGGCGCAGTCTAGCCCAGGGCGGCGTGAGCgggg                               >  8:272411/1‑147 (MQ=255)
                                                                                                                        aaCTGACCGACATCGCGCCTTGTCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATTGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCATGTCTGCTTCAGCGGCTGTcc                           <  8:323713/149‑1 (MQ=255)
                                                                                                                          cTGACCCACCTCGCGCCTTGGCCGCCGGCGGGGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTg                         <  2:178142/149‑1 (MQ=255)
                                                                                                                             aCCGACATCGCGCCTTGGCCGCCGGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGCACAGCCAGTGTGTGAGGTGGATGGTGAATACTTGGCAGGGGGTCCCGGGTCAGTCTAGCCCAGGTGTGCTTCAGCGGGGGTCCGGGtg                      >  2:214293/1‑148 (MQ=255)
                                                                                                                                      gcgcCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGt             <  6:329663/149‑1 (MQ=255)
                                                                                                                                           ttGGCCGCCGCAGGCGGCCAACGCTGTGCAGGGCGGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTc        >  4:188877/1‑149 (MQ=255)
                                                                                                                                                 gccTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTgg  <  7:288997/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCGTACCTGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGG  >  NC_002947/3466693‑3466986

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