Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,155,227 A→T 17.7% intergenic (‑37/+169) 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,2270AT17.7% 74.6 / 7.8 33intergenic (‑37/+169)PP_5398/PP_5399transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/14);  new base T (4/2);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

CCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTACTGCAG  >  NC_002947/6155148‑6155365
                                                                               |                                                                                                                                          
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                                                <  5:203274/119‑1 (MQ=17)
                                                     gCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCaga                                                >  6:203274/1‑119 (MQ=17)
                                                     gCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAgg                  >  8:115021/1‑149 (MQ=25)
                                                     gCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAgg                  >  8:110724/1‑149 (MQ=25)
                                                       attAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTg                >  4:201372/3‑149 (MQ=25)
                                                            gggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCagag                                               >  5:293957/1‑113 (MQ=255)
                                                            gggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCagag                                               <  6:293957/113‑1 (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)
                                                                               |                                                                                                                                          
CCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTACTGCAG  >  NC_002947/6155148‑6155365

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: