Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,155,242 C→G 19.7% intergenic (‑52/+154) 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,2420CG19.7% 78.7 / 11.1 35intergenic (‑52/+154)PP_5398/PP_5399transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (13/15);  new base G (5/2);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTACTGCAGAATTTCCTCAAAGTTTCCAA  >  NC_002947/6155148‑6155385
                                                                                              |                                                                                                                                               
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                                                                           <  5:319722/149‑1 (MQ=255)
                tttCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACgg                                                                           <  3:142354/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)
                                                       attAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTg                                    >  4:201372/3‑149 (MQ=25)
                                                            gggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCagag                                                                   <  6:293957/113‑1 (MQ=255)
                                                            gggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCagag                                                                   >  5:293957/1‑113 (MQ=255)
                                                                 aTCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAg                                                <  2:314517/127‑1 (MQ=255)
                                                                 aTCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAg                                                >  1:314517/1‑127 (MQ=255)
                                                                    gcgTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTACTGCa                       <  5:20176/149‑1 (MQ=255)
                                                                     cgTCGGCTTATAAATACGCGGACACGATCGCGCTTAATGCTGGAGTTCGGAAGGTGAGTTCCGCGGACGGTTTCTCGCTGTTTGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGTCGCAGGGAAATCAGGTGCGGAAGTACTGCAg                      <  1:188534/149‑1 (MQ=25)
                                                                            ttAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAgg                                               >  6:41012/1‑117 (MQ=255)
                                                                            ttAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAgg                                               <  5:41012/117‑1 (MQ=255)
                                                                                aaaTACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTACTGCAGAATTTCCTCa            >  6:42939/1‑148 (MQ=255)
                                                                                         gACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTACTGCAGAATTTCCTCAAAGTTTCCaa  <  5:42939/149‑1 (MQ=255)
                                                                                              |                                                                                                                                               
CCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTACTGCAGAATTTCCTCAAAGTTTCCAA  >  NC_002947/6155148‑6155385

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: