Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,053,629 T→A 11.8% intergenic (‑29/+37) 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,6290TA11.8% 86.5 / 4.0 34intergenic (‑29/+37)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 = 9.69e-01

AGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTGG  >  NC_002947/2053488‑2053775
                                                                                                                                             |                                                                                                                                                  
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                              <  8:280849/138‑1 (MQ=255)
                                                                                                                          tCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAg                              >  7:280849/1‑138 (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)
                                                                                                                                             |                                                                                                                                                  
AGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTGG  >  NC_002947/2053488‑2053775

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: