Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,311,740 T→A 8.9% intergenic (‑202/+76) PP_1144 ← / ← rapA membrane protein/RNA polymerase‑associated protein RapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,311,7400TA 8.9% 120.3 / 4.4 45intergenic (‑202/+76)PP_1144/rapAmembrane protein/RNA polymerase‑associated protein RapA
Reads supporting (aligned to +/- strand):  ref base T (25/16);  new base A (2/2);  total (27/18)
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.92e-01

AGGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCG  >  NC_002947/1311599‑1311888
                                                                                                                                             |                                                                                                                                                    
agGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCg                                                                                                                                              >  2:821664/1‑150 (MQ=255)
  gCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGt                                                                                                                                            >  2:135520/1‑150 (MQ=255)
      tataGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGAt                                                                                                                                        >  1:182836/1‑150 (MQ=255)
                  tACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCttt                                                                                                                            <  2:989482/150‑1 (MQ=255)
                    cTTAAGCCTTCCTGACTAATGGCACAACGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGc                                                                                                                          <  2:299188/150‑1 (MQ=255)
                    cTTAAGCCTTCCTGACAAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCACTAATGCCTTCCTGTGGGAGCCGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGc                                                                                                                          <  1:105060/150‑1 (MQ=255)
                          ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtg                                                                                                                     >  1:1260619/1‑149 (MQ=255)
                          ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtg                                                                                                                     <  2:1260619/149‑1 (MQ=255)
                          ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGT                                                                                                                    >  2:813688‑M1/1‑150 (MQ=255)
                               cTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCG                                                                                                               >  1:1168395‑M1/1‑150 (MQ=255)
                                      aTGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGc                                                                                                        >  1:530368/1‑150 (MQ=255)
                                          cacaATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCACCATCCACTCCGGGGGATGCCCCCGGCTTTTCCGGTGTTTGCGGGCGGGCCCGc                                                                                                    >  2:959501/1‑150 (MQ=255)
                                           acaATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTTTTCGCGGGCGTGCCCGCt                                                                                                   >  2:209780/1‑150 (MQ=255)
                                            caATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTc                                                                                                  <  2:953773/150‑1 (MQ=255)
                                                             tGGTTATAACTCCAATCGCAGTGAGGAATAATTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgggg                                                                                 >  2:789530/1‑150 (MQ=255)
                                                                         cAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcc                                                                     >  1:172502/1‑150 (MQ=255)
                                                                          aaTCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcca                                                                    <  2:182836/150‑1 (MQ=255)
                                                                          aaTCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcca                                                                    >  2:511853/1‑150 (MQ=255)
                                                                               cAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGccac                                                               >  1:569326/1‑150 (MQ=255)
                                                                                aGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGccacc                                                              <  2:1168395/150‑1 (MQ=255)
                                                                                  tGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAg                                                            >  2:1114511/1‑150 (MQ=255)
                                                                                       aaTAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGGCTGCGCTAGCCAGCCACCAGCAcgc                                                       >  1:1284811/1‑150 (MQ=255)
                                                                                         tAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCcac                                                                                              <  2:302419/109‑1 (MQ=255)
                                                                                         tAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCcac                                                                                              >  1:302419/1‑109 (MQ=255)
                                                                                                     ccTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAgg                                         <  1:821664/150‑1 (MQ=255)
                                                                                                     ccTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAgg                                         >  2:221226/1‑150 (MQ=255)
                                                                                                           gCCTTCCTGTGGGAGCGGGCGAGCCGGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAgg                                   <  1:209780/150‑1 (MQ=255)
                                                                                                           gCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAgg                                   >  1:25794/1‑150 (MQ=255)
                                                                                                              ttCCTGTGGGAGCGGGCGAGCCCGACCACCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                >  1:359420/1‑150 (MQ=255)
                                                                                                                ccTGTGGGAGGGGGCGAGCTCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                              <  1:959501/150‑1 (MQ=255)
                                                                                                                ccTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                              <  1:135520/150‑1 (MQ=255)
                                                                                                                      ggAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTtctc                        <  2:98760/150‑1 (MQ=255)
                                                                                                                            ggCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCAt                  >  1:353720/1‑150 (MQ=255)
                                                                                                                             gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCTGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCAtt                 <  2:359420/150‑1 (MQ=255)
                                                                                                                               gAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGc               <  1:400163/150‑1 (MQ=255)
                                                                                                                                   ccGCGAAGCAAGCACTGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCaagcccgctcccacagggtaaggtgctgcggtttagactttgtggtacagctggctgccttcctggcggaac                                    <  2:662931‑M1/125‑73 (MQ=255)
                                                                                                                                   ccGCGAAGCAAGCACTGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCaagcccgctcccacagggtaaggtgctgcggtttagactttgtggtacagctggctgccttcctggcggaac                                    >  1:662931‑M1/1‑53 (MQ=255)
                                                                                                                                   ccGCGAAGCAAGCACTGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGCGCaagcccgctcccacagggtgaggtgctgcggtttagactttgtggtacagctggctgccttcctggcggaaccgctcggcctgctcacgcatgcctt           >  1:733249‑M1/1‑53 (MQ=255)
                                                                                                                                     tcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGcc         >  2:882004/2‑150 (MQ=255)
                                                                                                                                     gcgAAGCAAGCACTGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCaagcccgctcccacagggtgaggtgctgcggtttagactttgtggtacagctggctgccttcctggcggaaccgctcggcctgctcacgcatgccttgc         <  2:733249‑M1/150‑100 (MQ=255)
                                                                                                                                        tggCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATc      >  2:859016/3‑150 (MQ=255)
                                                                                                                                         ggCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCt     >  1:270109/2‑150 (MQ=255)
                                                                                                                                            aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGcc              <  2:314221/138‑1 (MQ=255)
                                                                                                                                            aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGcc              >  1:314221/1‑138 (MQ=255)
                                                                                                                                            aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcg  >  1:761240/1‑150 (MQ=255)
                                                                                                                                             tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagc                                                                      <  1:839141/81‑1 (MQ=255)
                                                                                                                                             tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagc                                                                      >  2:839141/1‑81 (MQ=255)
                                                                                                                                             tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTag                                                                       >  2:1080943/1‑80 (MQ=255)
                                                                                                                                             tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTag                                                                       <  1:1080943/80‑1 (MQ=255)
                                                                                                                                             tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGgcgc                                      >  2:959353/1‑113 (MQ=255)
                                                                                                                                             tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGgcgc                                      <  1:959353/113‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
AGGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCG  >  NC_002947/1311599‑1311888

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