Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,530,025 T→A 10.3% intergenic (+31/‑137) 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,0250TA10.3% 115.9 / 2.0 45intergenic (+31/‑137)PP_4018/topBacetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (27/13);  new base A (3/2);  total (30/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATG  >  NC_002947/4529889‑4530164
                                                                                                                                          |                                                                                                                                           
gCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGGCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGGCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                 <  7:259089/149‑1 (MQ=255)
 cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                 <  4:119521/148‑1 (MQ=255)
      cacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                           >  8:3565/1‑149 (MQ=255)
             ttCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCCGGTTTACCAGCGAAGAGGCCGGCTCAGGc                                                                                                                    >  7:237711/1‑149 (MQ=255)
              gctaggggCGGCTGGGCTTTCGCGTTGGCCTGGTGGAATTCCTCGATGGTGGCATCCGCCATGTCTCCATGGTGCGCGACGCGCGCGCCTGATCCATCACCGGCAGGGGCCTCTTCGCGGTTTTCCCCGCTATTAGGCCGGCTCAGGc                                                                                                                    <  2:303560/141‑1 (MQ=255)
                           gggCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCg                                                                                                      >  2:25862/1‑149 (MQ=255)
                               tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGGTTTACCGCGAAGAGGGCGGCTCAGGCGAATTACATCCCCGcgg                                                                                                   >  1:169367/1‑146 (MQ=255)
                                    cgTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGGGAAGAGGCCGGCTCAGGCGATATACATACCCGCTGACaaa                                                                                              >  8:79803/1‑148 (MQ=255)
                                    cgTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGGTTAACCGCGAAGAGGCCGGCTCAGGCGCTATACATCCCCGCTGACaaa                                                                                              >  1:172869/1‑148 (MQ=255)
                                     gTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGGTTACCCGCGAAGAGGCCGGCTCAGGCGATATAAATCCACGCTGACaaaa                                                                                             >  5:27829/1‑148 (MQ=255)
                                      tAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGAAGACGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACaaaa                                                                                             >  1:293231/1‑148 (MQ=255)
                                        agcttcGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAAATACATCCCCGCTGACAAAAAAcc                                                                                         >  7:139314/6‑149 (MQ=255)
                                         gccgCGAGTAATTCCTGGAAGCTGGCAGCCCGCATGGCGACTTGGGGCGCGACTGGCGCGCCTTATCCATCCCCTCCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCt                                                                                        <  6:163245/146‑1 (MQ=255)
                                           aGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTg                                                                                       >  1:90896/1‑148 (MQ=255)
                                               aGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcat                                                                                  >  5:221431/1‑149 (MQ=255)
                                                                 gaaTCCCGAATGTCGAATTGGTGGGGGCCTCGGGCGCCTTATCCATCACCTGCACGGGCATCTGGGGGTTTTTACCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACtt                                                                 <  4:317757/146‑1 (MQ=255)
                                                                       cGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCGTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCctct                                                           >  1:268700/1‑149 (MQ=255)
                                                                        gCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTgc                                                         >  5:314152/1‑149 (MQ=255)
                                                                                aCATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGc                                                 <  7:3565/149‑1 (MQ=255)
                                                                                aCATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCATCGCGGGTTTACCCGCGAAGACGCCGGCACAGGCGATATAAAACCCCCCTTACAAAAAACCAGTTCATCCATCCAGATAAAACAAGCCTATGCGCCCCCGc                                                 >  4:126747/1‑149 (MQ=255)
                                                                                           gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCCCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgcgatgcag                                      >  3:239468/1‑143 (MQ=255)
                                                                                            cacggccgcgcgTCATCCATCCCCTGCACTGCCCTCTGCGGGGGTTTCCCGGCTAATTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCt                                      <  6:27829/143‑1 (MQ=255)
                                                                                                cgcgcgCCTTACCCAACACCTGCCCGGGCCTGTTCGCGGTTTTCCCCGCTACGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATgg                                  <  8:139314/149‑1 (MQ=255)
                                                                                                cgcgcgCCTGATCCATCACCTGCACTTGCCTCTTCGGGGTTTAACCGGGGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATgg                                  <  1:25862/149‑1 (MQ=255)
                                                                                                     gcCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACTTCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCccac                             >  4:37589/1‑149 (MQ=255)
                                                                                                       cTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCccacca                           >  3:342335/1‑149 (MQ=255)
                                                                                                          ttcaaTCACGTGCACTGGCCTCTTCGCGGGTTTACCCGCTATGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                        <  2:293231/145‑1 (MQ=255)
                                                                                                          aTCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                        >  3:308687/1‑149 (MQ=255)
                                                                                                              aTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGAt                    >  6:143828/1‑149 (MQ=255)
                                                                                                              aTCACCTGCACTGGCCTCTTAGCGGGTTTAACCGCGAAGATGCCCGCAAAGGCGATATACATACCCGCGGACAAAAAAACGGTACAAACATCCAGATAAAACATGGCTATGCGCCCCCCCTTGCGCAACCTATGGGCCACcccccccag                    >  4:207561/1‑146 (MQ=255)
                                                                                                                 aCCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                 <  4:308687/149‑1 (MQ=255)
                                                                                                                 aCCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                 <  2:90896/149‑1 (MQ=255)
                                                                                                                   cTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                        >  5:121668/1‑140 (MQ=255)
                                                                                                                   cTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                        <  6:121668/140‑1 (MQ=255)
                                                                                                                   cTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGc               >  6:305570/1‑149 (MQ=255)
                                                                                                                            gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCCTCCTATGGCCCACCCCCCCGATGAAGCGTCCCCgg       >  7:103191/1‑149 (MQ=255)
                                                                                                                              ctctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGa                     >  1:312905/1‑133 (MQ=255)
                                                                                                                              ctctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGa                     <  2:312905/133‑1 (MQ=255)
                                                                                                                               tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa    >  5:213824/1‑149 (MQ=255)
                                                                                                                               tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa    >  8:76912/1‑149 (MQ=255)
                                                                                                                               tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt   >  4:370843/1‑149 (MQ=255)
                                                                                                                               tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTTACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa    >  1:275871/1‑149 (MQ=255)
                                                                                                                                 ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATg  >  1:251403/1‑149 (MQ=255)
                                                                                                                                     cgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGcc     >  6:267852/1‑142 (MQ=255)
                                                                                                                                     cgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGcc     <  5:267852/142‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATG  >  NC_002947/4529889‑4530164

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