Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,155,204 A→T 20.4% intergenic (‑14/+192) PP_5398 ← / ← PP_5399 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,155,2040AT20.4% 41.3 / 5.6 24intergenic (‑14/+192)PP_5398/PP_5399transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base T (3/2);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGG  >  NC_002947/6155148‑6155349
                                                        |                                                                                                                                                 
ccTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCAGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGGTCGGAAGGGGAGTTCCCCGGACGCTTACTCGCTg                                                       <  3:184923/149‑1 (MQ=255)
 cTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGTAGTTCGGAAGGTGAGTTCCACGGACGCTTACTCGCTg                                                       >  7:288532/1‑148 (MQ=21)
 cTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCAGCTTAAAAATACGCGTACACCATCGCCCTTAATGATGGAGTTCGGAAGGTGAGTTCCCAGGACGCTTACTCGCTg                                                       >  7:292555/1‑148 (MQ=255)
  tGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTa                                                     >  8:250958/1‑149 (MQ=35)
          gAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTc                                              >  2:154212/1‑148 (MQ=255)
           aaCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTc                                              >  6:20176/1‑147 (MQ=255)
            aCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTc                                              >  8:25173/1‑146 (MQ=255)
                tttCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACgg                                       <  3:142354/149‑1 (MQ=255)
                tttCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACgg                                       <  5:319722/149‑1 (MQ=255)
                   cGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGa                                     <  4:34661/148‑1 (MQ=255)
                       tAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCaga                                <  7:250958/149‑1 (MQ=255)
                        aaGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCagag                               <  4:61446/149‑1 (MQ=255)
                          gagCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGtt                             <  6:322615/149‑1 (MQ=255)
                          gagCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCAGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGtt                             >  1:93971/1‑149 (MQ=255)
                          gagCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGGCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGtt                             <  6:324792/149‑1 (MQ=255)
                           agCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTg                            >  4:64758/1‑149 (MQ=255)
                                          tGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCa             <  7:25173/149‑1 (MQ=255)
                                             aaaTCCTGGCATTAAGGGTGATCGCTTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCATAGTTGTACTTTCTGGACGCAgg           <  4:1351/148‑1 (MQ=255)
                                              aaTCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAggg          >  3:270123/1‑148 (MQ=255)
                                                 ccTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACGATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACGCGGTGTATGCGTTGGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAAt      >  5:250067/1‑149 (MQ=255)
                                                     gCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCaga                                >  6:203274/1‑119 (MQ=17)
                                                     gCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCaga                                <  5:203274/119‑1 (MQ=17)
                                                     gCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAgg  >  8:110724/1‑149 (MQ=25)
                                                     gCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAgg  >  8:115021/1‑149 (MQ=25)
                                                        |                                                                                                                                                 
CCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGG  >  NC_002947/6155148‑6155349

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: