Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,053,631 T→A 11.8% intergenic (‑31/+35) PP_1828 ← / ← PP_1829 hypothetical protein/alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,053,6310TA11.8% 88.3 / 4.2 34intergenic (‑31/+35)PP_1828/PP_1829hypothetical protein/alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (19/11);  new base A (2/2);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

AGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTGGGGCT  >  NC_002947/2053488‑2053779
                                                                                                                                               |                                                                                                                                                    
agggcggcCTGGCGGAGCTGCTGGAACTGTGCGGCGGGCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCGGGGGTTCTTCGGCCGGCTCGCCGGCGGGGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCggca                                                                                                                                                 <  6:323005/143‑1 (MQ=255)
                     tGGAATTGTTGGGCGGTCCGCTCGGCCCAGGGCCGGGCGGGGGTTTCCTGGGGGTCTTCGGCCGGCTCGCCGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCCATGCCTGCCGATACAGCAGGCAtag                                                                                                                                          <  4:114255/135‑3 (MQ=255)
                     tGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCCATGCCTGCCGATACAGCAGGCAtag                                                                                                                                          >  3:114255/1‑133 (MQ=255)
                                ggcggTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTa                                                                                                                 >  7:72691/1‑149 (MQ=255)
                                       ttgcgggCCCAGGGCTGGGGGGGTGTTTCGGGGGGTTGTTCGGCCGCCGCGCCGGGGGGGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGaa                                                                                                          <  6:30492/144‑1 (MQ=255)
                                         cTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGaa                                                                                                          >  5:30492/1‑147 (MQ=255)
                                               cccgggccgggcGGGGGTTTCGTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGAc                                                                                                                                   <  6:203796/113‑1 (MQ=255)
                                               ccAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGAc                                                                                                                                   >  5:203796/1‑116 (MQ=255)
                                                                   ccTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATtc                                                                              >  1:301323/1‑149 (MQ=255)
                                                                       gggTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTAAGCCGAAGGCAAATGCTGGATGAACGTCACAATtcgctc                                                                          >  7:202442/1‑149 (MQ=255)
                                                                                        cGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCa                                                                     <  2:189161/137‑1 (MQ=255)
                                                                                        cGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCa                                                                     >  1:189161/1‑137 (MQ=255)
                                                                                           cggcggTGGGGGAGGTGGTCATCAGGTCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTAGCCAAGCTCACCCGCCAGCg                                                      >  6:80418/1‑149 (MQ=255)
                                                                                                     ggAGGTTGTCATCAGGTCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGTTGAACGTCACTATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCggg                                            <  5:80418/149‑1 (MQ=255)
                                                                                                            gtCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAg                                      <  2:301323/148‑1 (MQ=255)
                                                                                                               atcaGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTg                                        <  1:338375/143‑1 (MQ=255)
                                                                                                               atcaGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTg                                        >  2:338375/1‑143 (MQ=255)
                                                                                                                   gggCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGt                              >  3:295089/1‑149 (MQ=255)
                                                                                                                          tCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAg                                  >  7:280849/1‑138 (MQ=255)
                                                                                                                          tCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAg                                  <  8:280849/138‑1 (MQ=255)
                                                                                                                          tCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctg                       >  6:164117/1‑149 (MQ=255)
                                                                                                                             tCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGccc                                                                      <  5:65051/99‑1 (MQ=255)
                                                                                                                             tCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGccc                                                                      >  6:65051/1‑99 (MQ=255)
                                                                                                                             tCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGctgctg                    >  5:261949/1‑149 (MQ=255)
                                                                                                                             tCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGctgctg                    >  3:5884/1‑149 (MQ=255)
                                                                                                                             tCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGctgctg                    >  3:84073/1‑149 (MQ=255)
                                                                                                                               aGCCATGCCTGCCGATACAGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctg                       <  7:309374/144‑1 (MQ=255)
                                                                                                                               aGCCATGCCTGCCGATACAGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctg                       >  8:309374/1‑144 (MQ=255)
                                                                                                                                    tGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGAAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATg             >  5:294503/1‑149 (MQ=255)
                                                                                                                                     gCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGc            >  3:76040/1‑149 (MQ=255)
                                                                                                                                      ccTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGcc           >  6:271679/1‑149 (MQ=255)
                                                                                                                                      ccTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGcc           >  3:78388/1‑149 (MQ=255)
                                                                                                                                      ccTGATGTTTCGGGAGTCATGGCTGAGAGCTTTGCCTGTTTGACTTACGCCGAAGTCAAATGCTGGATGAACGTCACAATTCGCTCGGCCAGCTCACTCGCCATCGGCAGCTCTTGGTTGAGGTAGCTGTCGAGCTGCTGGCTCATGc            <  6:294503/148‑1 (MQ=255)
                                                                                                                                           tgTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTgg      <  4:78388/149‑1 (MQ=255)
                                                                                                                                               tCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTGGGGCt  >  7:259055/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTGGGGCT  >  NC_002947/2053488‑2053779

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