Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,041 Δ1 bp 10.0% intergenic (+130/+59) PP_5440 → / ← PP_t72 hypothetical protein/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,0410A.10.0% 149.1 / 11.4 50intergenic (+130/+59)PP_5440/PP_t72hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base . (2/3);  total (20/30)
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.04e-01

CATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGA  >  NC_002947/343896‑344185
                                                                                                                                                 |                                                                                                                                                 
cATGTCGGACGATTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCAc                                                                                                                                                  >  3:216866/1‑147 (MQ=255)
  tGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAg‑g                                                                                                                                             <  6:125558/149‑1 (MQ=255)
      ggACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCggg                                                                                                                                         >  6:330383/1‑149 (MQ=255)
        aCGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGtt                                                                                                                                       >  1:88327/1‑149 (MQ=255)
            gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCtt                                                                                                                                   <  1:300796/149‑1 (MQ=255)
            gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCtt                                                                                                                                   <  3:140183/149‑1 (MQ=255)
               aGGTTTGTCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGTTTTttgtt                                                                                                                                <  7:122659/149‑1 (MQ=39)
               aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttgtt                                                                                                                                <  6:216942/149‑1 (MQ=255)
               aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttgt                                                                                                                                 <  2:63831/148‑1 (MQ=255)
                ggTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGAGTGCATGGCTAGCCTTCAGGTGGCGTTGTGAGGGCTCGTGGTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttgttg                                                                                                                               <  7:243341/149‑1 (MQ=255)
                 gTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTtgttgt                                                                                                                              <  4:216866/149‑1 (MQ=255)
                  tttGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCtt                                                                                                                                   <  7:293413/143‑1 (MQ=255)
                  tttGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCtt                                                                                                                                   >  8:293413/1‑143 (MQ=255)
                      gCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCt                                                                                                                          <  7:231341/148‑1 (MQ=255)
                       cccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTc                                                                                                                        <  8:289380/149‑1 (MQ=255)
                         cGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCgg                                                                                                                      <  3:180576/149‑1 (MQ=255)
                                     aTGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCaa                                                                                                          >  3:89239/1‑149 (MQ=255)
                                         cTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTTTTCTTCGACGCCGACTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑TCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                                      <  4:237876/149‑1 (MQ=255)
                                         cTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                                      <  5:330383/149‑1 (MQ=255)
                                             gggTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAg                                                                                                   >  4:170166/1‑148 (MQ=255)
                                                   ctcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCtt                                                                                             >  7:17387/1‑148 (MQ=255)
                                                      gATTGCATGGCTAGCCTTCAGGGTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGACGCTTGaa                                                                                          <  6:184688/148‑1 (MQ=255)
                                                               gCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGccc                                                                                <  2:88327/149‑1 (MQ=255)
                                                                cTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCa                                                                               <  2:169913/149‑1 (MQ=255)
                                                                  aGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCa                                                                               >  4:235204/1‑147 (MQ=255)
                                                                    ccTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                                  <  8:55379/94‑1 (MQ=255)
                                                                    ccTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                                  >  7:55379/1‑94 (MQ=255)
                                                                          ggTTGCGTTGGGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGt                                                                     <  7:213193/149‑1 (MQ=255)
                                                                                   gtgAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGTAGTCGAAc                                                                >  3:220191/1‑145 (MQ=255)
                                                                                   gtgAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAc                                                                <  4:220191/145‑1 (MQ=255)
                                                                                      aGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGa                                                          <  8:17387/148‑1 (MQ=255)
                                                                                           tcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAg                                                    >  7:250550/1‑149 (MQ=255)
                                                                                           tcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGAACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAg                                                    >  1:91123/1‑149 (MQ=255)
                                                                                            cGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAg                                                    >  2:167197/1‑148 (MQ=255)
                                                                                             gtgtTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGa                                                  >  8:27706/1‑149 (MQ=255)
                                                                                                   ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTc                                             >  3:320772/1‑148 (MQ=255)
                                                                                                     cGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                              <  1:137027/145‑1 (MQ=255)
                                                                                                      gACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATc                                         <  4:89239/149‑1 (MQ=255)
                                                                                                         gccgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCctc                                      <  3:189555/149‑1 (MQ=255)
                                                                                                         gccgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCctc                                      <  3:235204/149‑1 (MQ=255)
                                                                                                              cTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGctc                                 >  8:215926/1‑149 (MQ=255)
                                                                                                                              caGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATc                 <  8:250550/149‑1 (MQ=255)
                                                                                                                                gACAACAAGAAACCCGC‑CTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTg               >  3:230559/1‑149 (MQ=39)
                                                                                                                                 acaacaAGAAACCCGC‑CTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATcc                                        <  3:207019/123‑1 (MQ=38)
                                                                                                                                 acaacaAGAAACCCGC‑CTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATcc                                        >  4:207019/1‑123 (MQ=38)
                                                                                                                                 acaacaAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTgg              >  6:89322/1‑149 (MQ=255)
                                                                                                                                   aacaaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCGCTACCGACTGAGCTATCTGGGc            <  7:215926/149‑1 (MQ=255)
                                                                                                                                        aaaaCCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac         >  2:311511/1‑147 (MQ=255)
                                                                                                                                         aaaCCCGC‑CTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAg      <  4:230559/149‑1 (MQ=255)
                                                                                                                                             ccGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGa  <  7:27706/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGA  >  NC_002947/343896‑344185

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