Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,939,275 A→G 29.6% intergenic (‑221/+219) 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,2750AG29.6% 26.2 / 15.0 27intergenic (‑221/+219)PP_2570/PP_2572transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/8);  new base G (3/5);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

CTCTGTCCCGCGCACACGAACCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCG  >  NC_002947/2939226‑2939420
                                                 |                                                                                                                                                 
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)
                                             cccTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGGCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCGGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCcgc   <  5:121593/149‑1 (MQ=255)
                                               cTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCcgcg  >  6:142090/1‑148 (MQ=255)
                                                 |                                                                                                                                                 
CTCTGTCCCGCGCACACGAACCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCG  >  NC_002947/2939226‑2939420

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: