Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,506,031 G→T 14.8% intergenic (‑258/‑228) PP_5566 ← / → PP_5567 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,506,0310GT14.8% 81.5 / 2.0 27intergenic (‑258/‑228)PP_5566/PP_5567hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/7);  new base T (2/2);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

TGTGACGGGCCATTCGGCGATTGCCGGGTCGCCTCCGGGCGCCAGTGAGGAAGATGAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTGCAGGCAGTCGCATGTTTCCGCAAGC  >  NC_002947/3505904‑3506161
                                                                                                                               |                                                                                                                                  
tgtgACGGGCCATTCGGCGATTGCCGGGTCGCCTCCGGGCGCCAGTGAGGAAGATGAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCtt                                                                                                                            >  1:252661/1‑136 (MQ=255)
tgtgACGGGCCATTCGGCGATTGCCGGGTCGCCTCCGGGCGCCAGTGAGGAAGATGAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCtt                                                                                                                            >  1:338318/1‑136 (MQ=255)
       ggCCATTCGGCGATTGCCGGGTCGCCTCCGGGCGCCAGTGAGGAAGATGAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCgg                                                                                                                     <  2:338318/136‑1 (MQ=255)
          cATTCGGCGATTGCCGGGTCGCCTCCGGGCGCCAGTGAGGAAGATGAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACATGGCAGCTTTTGCCGgtgt                                                                                                                  >  2:47839/1‑136 (MQ=255)
            ttCGGCGATTGCCGGGTCGCCTCCGGGCGCCAGTGAGGAAGATGAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGttt                                                                                                                <  1:39031/136‑1 (MQ=255)
                                   cGGGCGCCAGTGAGGAAGATGAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCctggct                                                                                         <  1:304721/136‑1 (MQ=255)
                                             tGAGGAAGATGAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTgg                                                                               >  2:440701/1‑136 (MQ=255)
                                                       gAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAAcc                                                                     >  2:437859/1‑136 (MQ=255)
                                                       gAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCATGGCTGGCATTCTGGAATGTGAAcc                                                                     >  2:33194/1‑136 (MQ=255)
                                                            aTGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTc                                                                >  1:446560/1‑136 (MQ=255)
                                                                        cATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCa                                                    <  2:447706/136‑1 (MQ=255)
                                                                             ccTTGCAGGGCGGTATTGCCCCGTTTTCCTATGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATTCa                                               >  2:64269/1‑136 (MQ=255)
                                                                              cTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCat                                              >  1:508405/1‑136 (MQ=255)
                                                                                     ggCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCa                                       >  1:215422/1‑136 (MQ=255)
                                                                                          tATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCActtta                                  >  1:409087/1‑135 (MQ=255)
                                                                                               ccccGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCActttcctttc                             >  2:536054/1‑136 (MQ=255)
                                                                                                       tCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCGTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTgcaggca                     >  1:503836/1‑136 (MQ=255)
                                                                                                       tCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTAGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTgcaggca                     >  1:545251/1‑136 (MQ=255)
                                                                                                         cTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTGCAGGCAGt                   <  2:215422/136‑1 (MQ=255)
                                                                                                             gCCTGTAGTGAGTTGACATGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTGCAGGCAGTCGCa               <  1:47839/136‑1 (MQ=255)
                                                                                                                  tAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTACTTTCTGCAGGCAGTCGCATGttt          >  2:327148/1‑136 (MQ=255)
                                                                                                                     tGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTGCAGGCAGTCGCATGTTTCCg       >  1:294346/1‑136 (MQ=255)
                                                                                                                        gTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTGCAGGCAGTCGCATGTTTCCGCaa    <  1:437859/136‑1 (MQ=255)
                                                                                                                        gTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTGCAGGCAGTCGCATGTTTCCGCaa    <  1:33194/136‑1 (MQ=255)
                                                                                                                          tGACATGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTATTTCAATGCATGATGGCACTTTCCTTTCTGCAGGCAGTCGCATGTTTCCGCAAGc  >  1:87871/1‑135 (MQ=255)
                                                                                                                          tGACATGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTATTTCAATGCATGATGGCACTTTCCTTTCTGCAGGCAGTCGCATGTTTCCGCAAGc  <  2:87871/135‑1 (MQ=255)
                                                                                                                          tGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTGCACGCAGTCGCATGTTTCCGCAAGc  >  1:333367/1‑136 (MQ=255)
                                                                                                                               |                                                                                                                                  
TGTGACGGGCCATTCGGCGATTGCCGGGTCGCCTCCGGGCGCCAGTGAGGAAGATGAAGGATGGCGATATGACATTGCCTTGCAGGGCGGTATTGCCCCGTTTTCCTAGGCCTGTAGTGAGTTGACAGGGCAGCTTTTGCCGGTGTTTCCTTTGGGCCTTTCTGCCTGGCTGGCATTCTGGAATGTGAACCCTGTCGCATTTATTTCAATGCATGATGGCACTTTCCTTTCTGCAGGCAGTCGCATGTTTCCGCAAGC  >  NC_002947/3505904‑3506161

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: