Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,939,267 C→T 34.3% intergenic (‑213/+227) PP_2570 ← / ← PP_2572 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,939,2670CT34.3% 11.4 / 13.9 22intergenic (‑213/+227)PP_2570/PP_2572transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (8/6);  new base T (2/5);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CTCTGTCCCGCGCACACGAACCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCA  >  NC_002947/2939226‑2939414
                                         |                                                                                                                                                   
agttccgcagcgacgctgacattttgaaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGc                                          <  8:184646‑M1/117‑1 (MQ=255)
            acgctgacattttgaaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTcc                              >  7:184646‑M1/21‑149 (MQ=255)
                atctattttgaaaaagAAAAGGCGGGTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACa                          >  7:155340‑M1/17‑149 (MQ=255)
                     ttttgaaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTcctacgggc                       <  4:288278‑M1/136‑5 (MQ=255)
                     ttttgaaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACa                          <  8:155340‑M1/133‑1 (MQ=255)
                     ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCGACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTAc                     <  2:162785/149‑1 (MQ=17)
                     ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTAc                     >  1:102554/1‑149 (MQ=32)
                       tAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActg                   >  3:125595/1‑149 (MQ=33)
                       tAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActg                   >  4:209068/1‑149 (MQ=33)
                        aaCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATTTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgc                  <  1:175087/149‑1 (MQ=21)
                         aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActggt                 >  1:177754/1‑147 (MQ=18)
                         aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgct                 >  5:236535/1‑149 (MQ=35)
                         aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGCACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgct                 <  7:72265/149‑1 (MQ=25)
                         aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTAAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgct                 >  2:284803/1‑149 (MQ=25)
                          aaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTg                                        >  5:357777‑M1/7‑125 (MQ=255)
                          aaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTg                                        <  6:357777‑M1/119‑1 (MQ=255)
                           aaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgctg                <  6:113444‑M1/143‑1 (MQ=255)
                           aaaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGt               >  8:304477/1‑149 (MQ=37)
                            aaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTa              <  4:68346/149‑1 (MQ=38)
                                aaaaGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGaa          >  1:112921/1‑149 (MQ=40)
                                 aaaGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAAc         <  7:325382/149‑1 (MQ=255)
                                   aGGCGGCTGGCCCTACCGGGCAAACCGCTTTTTTTTTTTAAAAATGTCAGCGTCGCTGCGGAACTGAAGCCTTACACCTCATTTTTGCTGCGCTACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACgg       <  2:351773/149‑1 (MQ=14)
                                    ggcggcTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACggg      >  8:133047/1‑149 (MQ=255)
                                       ggcTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTc   <  7:65007/149‑1 (MQ=255)
                                        gcTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCa  >  4:89695/1‑149 (MQ=255)
                                         |                                                                                                                                                   
CTCTGTCCCGCGCACACGAACCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCA  >  NC_002947/2939226‑2939414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: