Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,530,026 T→A 29.5% intergenic (+32/‑136) PP_4018 → / → topB acetyltransferase/DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,530,0260TA29.5% 56.1 / 27.2 44intergenic (+32/‑136)PP_4018/topBacetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (19/12);  new base A (5/8);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4529886‑4530171
                                                                                                                                            |                                                                                                                                                   
ggTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCcgcg                                                                                                                                              >  8:162306/1‑148 (MQ=255)
 gggCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                            <  6:213470/147‑1 (MQ=255)
 gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                            <  5:1147/149‑1 (MQ=39)
 gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                            <  5:158298/149‑1 (MQ=255)
 gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                            <  6:201168/149‑1 (MQ=255)
 gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                            <  7:162306/149‑1 (MQ=255)
    cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGGGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                          <  6:245299/148‑1 (MQ=255)
    cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAagag                                                                                                                                         <  2:186821/149‑1 (MQ=255)
    cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAagag                                                                                                                                         <  1:158927/149‑1 (MQ=255)
    cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                          <  1:311347/148‑1 (MQ=255)
    cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                          <  1:180515/148‑1 (MQ=255)
      cgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGGGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGc                                                                                                                                       <  5:194584/149‑1 (MQ=255)
       gccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                          <  6:185702/145‑1 (MQ=255)
       gccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                          >  5:185702/1‑145 (MQ=255)
       gccacgcaGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCACGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGcc                                                                                                                                      <  6:413585/149‑1 (MQ=255)
            gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                          <  3:393293/140‑1 (MQ=255)
            gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                          >  4:393293/1‑140 (MQ=255)
                                            gcttcGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCt                                                                                                >  7:273853/5‑149 (MQ=255)
                                               gCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGTGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                             >  1:106328/1‑149 (MQ=255)
                                               gCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCTAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                             <  2:289548/149‑1 (MQ=255)
                                               gCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                             >  3:222618/1‑149 (MQ=255)
                                                  aGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTAATCGCGGGGTTAACCCCGTATAGGCCGG‑CTCAGGCGATATACAGCCCCGCGGACAAAAAACCTGTAca                                                                                           >  8:152587/1‑148 (MQ=255)
                                                          cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                     >  2:100480/1‑146 (MQ=255)
                                                          cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGGTTACCCGCGCAGAGGCCGG‑CTCAGGCGATATAAATCCCCGCTGACAAAAAACCTGGAcatccctc                                                                                     >  4:355149/1‑146 (MQ=255)
                                                                  tGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACt                                                                         >  1:50946/1‑149 (MQ=255)
                                                                     cATCCGGCATGTTTAACTGGTGCTCGACTCGCGCGCCTGTTCCATAACCTGACGTGGACTCTTCGCGGTTTTACCGGCGATGAGGCCTG‑CTCAGGCGATATACATCCCAGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGc                                                                      <  8:273853/149‑1 (MQ=255)
                                                                                           cgcgACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCAt                                                >  4:169849/1‑149 (MQ=255)
                                                                                             cgACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATcc                                              >  7:64089/1‑149 (MQ=255)
                                                                                                         cTGATCCATCACCTGCACTGGCCTCTTCGCCGGTTTACCCGCGAAGACGCCGG‑CTCAGGCGCTATACATACCCGCAGACAAAAAACGTGGACAACCCACCACAAAA‑AAAATGCCTCTGCGCCCCCCCATTCCCATCCAATGGCccccca                                  >  8:170893/1‑149 (MQ=255)
                                                                                                          tGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCccacccc                                 >  6:295806/1‑147 (MQ=255)
                                                                                                           gATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCaccacc                                >  6:369505/1‑149 (MQ=255)
                                                                                                            aTCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                               >  7:300575/1‑149 (MQ=255)
                                                                                                                 tCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGAt                           >  1:76607/1‑148 (MQ=255)
                                                                                                                   aCCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAACAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATg                          <  1:100480/148‑1 (MQ=255)
                                                                                                                   aCCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGa                         <  2:76607/148‑1 (MQ=255)
                                                                                                                             gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCgg              >  6:101981/1‑149 (MQ=255)
                                                                                                                             gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCgg              >  3:340597/1‑149 (MQ=255)
                                                                                                                             gCCTCTTCGCGGGTAAACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCCTCCTATGGCCCACCACCCCGATGCAGCGTCCCCgg              >  6:411457/1‑149 (MQ=255)
                                                                                                                               ctctTCGCGGGTAAACCCGCGAAGAGGCCGGTCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGc             >  1:247862/1‑149 (MQ=255)
                                                                                                                                  ttCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATg         >  2:373554/1‑149 (MQ=255)
                                                                                                                                  ttCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATg         >  3:167100/1‑149 (MQ=255)
                                                                                                                                      cgGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     <  3:169849/149‑1 (MQ=255)
                                                                                                                                         ggttACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑ACCTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCCCCACCCCGAGGAAGCGTCCCCGGCCAACCGccagc  >  5:314375/3‑146 (MQ=255)
                                                                                                                                         gTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAg                       <  7:290498/128‑1 (MQ=255)
                                                                                                                                         gTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAg                       >  8:290498/1‑128 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4529886‑4530171

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: ‑