Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,995,396 A→T 7.1% intergenic (‑117/+33) PP_3521 ← / ← PP_3522 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,995,3960AT 7.1% 166.2 / 4.2 56intergenic (‑117/+33)PP_3521/PP_3522membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base T (2/2);  total (29/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 = 8.39e-01

GTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCAGGCC  >  NC_002947/3995250‑3995543
                                                                                                                                                  |                                                                                                                                                   
gTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAAc                                                                                                                                                  <  1:1808130/150‑1 (MQ=255)
       atcttGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTTTACCcgcgaa                                                                                                                                           <  1:369299/146‑3 (MQ=255)
            ggTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACgg                                                                                                                                      <  1:2099253/150‑1 (MQ=255)
              tCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGTGGGTAAACCCGCGCCCACGGGc                                                                                                                                    <  1:900777/150‑1 (MQ=255)
                              ggTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCtttt                                                                                                                    <  1:1966864/150‑1 (MQ=255)
                                  tgGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTCAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAgg                                                                                                                <  1:1847469/150‑1 (MQ=255)
                                         ctctGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAg                                                                                                         >  2:2174717/1‑150 (MQ=255)
                                                  ggATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGAcggc                                                                                                >  2:2158871/1‑150 (MQ=255)
                                                     ttgttgTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACggcggc                                                                                             <  1:858958/150‑1 (MQ=255)
                                                           ttATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACAAACGCCCACGAACAACCGAGCCGGCTTTTCAGGATTTGAAAAAAACGAAGGCTAATga                                                                                       >  2:10079/1‑149 (MQ=255)
                                                             ataGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCtt                                                                                     <  1:115854/150‑1 (MQ=255)
                                                                          gTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAg                                                                        >  2:591709/1‑150 (MQ=255)
                                                                           tCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGc                                                                       >  2:2078292/1‑150 (MQ=255)
                                                                             aGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTc                                                                                    <  2:2195757/135‑1 (MQ=255)
                                                                             aGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTc                                                                                    >  1:2195757/1‑135 (MQ=255)
                                                                              ggACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATg                                                                    >  1:2079258/1‑150 (MQ=255)
                                                                                   ggCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTc                                                               <  2:2079258/150‑1 (MQ=255)
                                                                                    gCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGcc                                                                                                                           <  2:1087058/89‑1 (MQ=255)
                                                                                    gCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGcc                                                                                                                           >  1:1087058/1‑89 (MQ=255)
                                                                                    gCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTcc                                                              >  1:1356536/1‑150 (MQ=255)
                                                                                       ccGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGt                                                           <  1:10079/150‑1 (MQ=255)
                                                                                       ccGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGt                                                           <  1:2078292/150‑1 (MQ=255)
                                                                                        cGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTg                                                          >  1:2180228/1‑150 (MQ=255)
                                                                                         gATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGAc                                                                                                                                 <  1:2315805/78‑1 (MQ=255)
                                                                                         gATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGAc                                                                                                                                 >  2:2315805/1‑78 (MQ=255)
                                                                                          atatACAAGGAAGAACCACCGTGAGGATGGCGCCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGcc                                                        <  1:1196805/150‑1 (MQ=255)
                                                                                          atatACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGcc                                                        >  1:1390685/1‑150 (MQ=255)
                                                                                           tataCAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTTTTCGCGGGTAAACCcgc                                                                                                                                              >  1:1805962/1‑63 (MQ=255)
                                                                                           tataCAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCcgc                                                                                                                                              <  2:1805962/63‑1 (MQ=255)
                                                                                              aCAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTTTACCcgcg                                                                                                                                             >  2:2019188/1‑61 (MQ=25)
                                                                                              aCAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTTTACCcgcg                                                                                                                                             <  1:2019188/61‑1 (MQ=25)
                                                                                              aCAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGcc                                                    >  2:2111167/1‑150 (MQ=255)
                                                                                                aaGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTccacca                                                                         >  2:1947365/1‑127 (MQ=255)
                                                                                                aaGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTccacca                                                                         <  1:1947365/127‑1 (MQ=255)
                                                                                                aaGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATTCCGTCCAGTGACAGCCgg                                                  >  2:448645/1‑150 (MQ=255)
                                                                                                 aGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGCCCGAGCCGGCTTTTCTGGTTTTGAGCAAGACGGCGGCTCCTACTTGGATCTCCACCAGCATGCCGTCAATTGCCAGCCggg                                                 >  2:378926/1‑150 (MQ=38)
                                                                                                 aGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCggg                                                 <  1:880599/150‑1 (MQ=255)
                                                                                                          caccGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCAACAGCATGCCCTCCAGTGCCATCCGGGGCACCGGTa                                        >  1:401619/1‑150 (MQ=255)
                                                                                                           accGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTAt                                       <  1:2158871/150‑1 (MQ=255)
                                                                                                                 aGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGt                                 <  1:2174717/150‑1 (MQ=255)
                                                                                                                 aGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGt                                 <  1:591709/150‑1 (MQ=255)
                                                                                                                  ggATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTg                                >  2:160951/1‑150 (MQ=255)
                                                                                                                   gATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGc                               >  1:1886862/1‑150 (MQ=255)
                                                                                                                     tGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAg                             >  2:1684542/1‑150 (MQ=255)
                                                                                                                                aCGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCAt                  <  2:2180228/150‑1 (MQ=255)
                                                                                                                                  gCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGt                       >  2:2296784/1‑143 (MQ=255)
                                                                                                                                  gCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGt                       <  1:2296784/143‑1 (MQ=255)
                                                                                                                                  gCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCAtgt                >  2:1586265/1‑150 (MQ=255)
                                                                                                                                    cGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCa                                                                                                                  <  2:2024970/50‑1 (MQ=255)
                                                                                                                                    cGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCa                                                                                                                  >  1:2024970/1‑50 (MQ=255)
                                                                                                                                    cGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTgc              >  2:192667/1‑150 (MQ=255)
                                                                                                                                      gtgggcgcgGGTTTACCCGCGAACACCGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTgcgc            >  1:290385/6‑150 (MQ=255)
                                                                                                                                            gcgGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCa      <  1:2111167/150‑1 (MQ=255)
                                                                                                                                                gTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCAGGcc  >  1:663940/1‑150 (MQ=255)
                                                                                                                                                gTAAACCCGCGCCCACGGGCGACCGAGACGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCa                                                             >  2:2167788/1‑91 (MQ=255)
                                                                                                                                                gTAAACCCGCGCCCACGGGCGACCGAGACGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCa                                                             <  1:2167788/91‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCAGGCC  >  NC_002947/3995250‑3995543

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