Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,652,937 T→G 11.8% intergenic (‑38/+123) PP_4961 ← / ← PP_4962 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,652,9370TG11.8% 127.9 / 10.5 51intergenic (‑38/+123)PP_4961/PP_4962lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/21);  new base G (3/3);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GGTGTCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTG  >  NC_002947/5652789‑5653081
                                                                                                                                                    |                                                                                                                                                
ggTGTCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTtct                                                                                                                                                  <  7:128513/149‑1 (MQ=255)
  tgtCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCg                                                                                                                                                >  1:261801/1‑149 (MQ=255)
  tgtCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCg                                                                                                                                                <  2:261801/149‑1 (MQ=255)
     cGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGc                                                                                                                                             <  6:36210/149‑1 (MQ=255)
      gCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCt                                                                                                                                            >  6:260555/1‑149 (MQ=255)
      gCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCt                                                                                                                                            >  3:26523/1‑149 (MQ=255)
            ccAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCGGCCCTTCGCGGGCTCGCCCg                                                                                                                                      >  7:138440/1‑149 (MQ=255)
            ccAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCg                                                                                                                                      <  8:135185/149‑1 (MQ=255)
            ccAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCg                                                                                                                                      >  4:66077/1‑149 (MQ=255)
              aaCGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCt                                                                                                                                    <  6:181944/149‑1 (MQ=255)
              aaCGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCt                                                                                                                                    >  5:181944/1‑149 (MQ=255)
               aCGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTc                                                                                                                                   <  6:184115/149‑1 (MQ=255)
                 gggTCCAGTGGGCTGCAGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCGAGGGCCGTCACGGCCAAAAGCGCTTTCATTGAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCTCGCCCGCTccc                                                                                                                                 <  3:99363/149‑1 (MQ=39)
                   gTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCcac                                                                                                                               <  7:24575/149‑1 (MQ=255)
                              tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTg                                                                                                                     >  2:83830/1‑148 (MQ=255)
                               gCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCAGGGCCGTGACGGGCAAAAGCGCGTTGATTCAGTTTCCGCGATTTGGGGCTGCAGTGCAGCCCATCTCGGGCTAGCCCCCAGCCCCAAGGACAGTGCAg                                                                                                                  >  5:306678/1‑149 (MQ=255)
                                  ggggCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGtc                                                                                                                >  2:168411/1‑149 (MQ=255)
                                     gCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTctct                                                                                                             <  6:272568/149‑1 (MQ=255)
                                        ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCtgtg                                                                                                          <  7:26932/149‑1 (MQ=255)
                                         gTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                         <  2:104448/149‑1 (MQ=255)
                                             cgcTGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCCGCTGTGCGGCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGc                                                                                                     <  7:263089/149‑1 (MQ=21)
                                                cGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCGGCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCggg                                                                                                  <  8:138440/149‑1 (MQ=21)
                                                     gcagcaGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGc                                                                                             >  1:48849/1‑149 (MQ=255)
                                                      atgcagTGTGCATCCTGCCTGGGTCGCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGc                                                                                             <  6:306678/145‑1 (MQ=38)
                                                                cATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGc                                                                                  >  7:159144/1‑149 (MQ=255)
                                                                 aTCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGc                                                                                                                                             <  8:90661/89‑1 (MQ=255)
                                                                                ttttggCCGTGACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCctt                                                                  >  4:265091/5‑149 (MQ=255)
                                                                                 cAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttc                                                                 >  2:45683/1‑149 (MQ=255)
                                                                                    ggCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCTGTGCGGcccttcgcgg                                                              <  2:48849/149‑8 (MQ=38)
                                                                                      ccGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTctct                                                                                                             <  5:201919/100‑1 (MQ=255)
                                                                                      ccGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTctct                                                                                                             >  6:201919/1‑100 (MQ=255)
                                                                                          cACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGAc                                                         >  2:80428/1‑148 (MQ=255)
                                                                                                cAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAAt                                                  >  6:67487/1‑149 (MQ=255)
                                                                                                cAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCGTCTTCAAGACCACaaa                                                   >  3:229373/1‑148 (MQ=255)
                                                                                                  aaaGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                <  5:67487/149‑1 (MQ=255)
                                                                                                  aaaGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                <  1:80428/149‑1 (MQ=255)
                                                                                                              ttcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTATCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCttcttc                                                              >  5:237278/1‑123 (MQ=255)
                                                                                                                 aGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                <  1:15257/134‑1 (MQ=255)
                                                                                                                 aGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                >  2:15257/1‑134 (MQ=255)
                                                                                                                 aGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCACCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAAc                                              >  3:108093/1‑136 (MQ=255)
                                                                                                                    ttCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTAt                               >  7:181341/1‑148 (MQ=255)
                                                                                                                     tCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTa                             >  2:210809/1‑149 (MQ=255)
                                                                                                                       cTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGTGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAcc                           >  4:171618/1‑149 (MQ=255)
                                                                                                                        tCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAccc                          >  3:192139/1‑149 (MQ=255)
                                                                                                                               ggggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGTAGCGGGCTAGCCCGCTAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCaa                    <  1:45683/148‑1 (MQ=255)
                                                                                                                                gggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAgg                  <  8:159144/149‑1 (MQ=255)
                                                                                                                                    tgctgGGCCGCCCGTCTCGGGCTGGCCCGCTCCCACCGGTACAGTGCTGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCaa              <  4:229373/149‑1 (MQ=38)
                                                                                                                                     gctgTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCt                                 >  5:77047/1‑129 (MQ=255)
                                                                                                                                     gctgTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCt                                 <  6:77047/129‑1 (MQ=255)
                                                                                                                                        gtgCAGCCCTTCTCGGGCTCGCCCGCTCCCCCAGGTACAGTGCAGTCTCGGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGAGCACCACTCAACGTTCAGGTTGGCTATTACCCGCTCTAGGGGAACCCg          >  5:278295/1‑149 (MQ=38)
                                                                                                                                              cccTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCt    >  1:47424/1‑149 (MQ=255)
                                                                                                                                                cTTCTCGGGCTCGCCCGCACCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTg  <  1:168411/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGTGTCGCAGACCCAACGGGTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTG  >  NC_002947/5652789‑5653081

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