Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,530,026 T→A 26.1% 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,0260TA26.1% 52.2 / 18.4 46intergenic (+32/‑136)PP_4018/topBacetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (18/16);  new base A (3/9);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCT  >  NC_002947/4529887‑4530169
                                                                                                                                           |                                                                                                                                               
gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                        <  3:286481/149‑1 (MQ=255)
gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                        <  3:23285/149‑1 (MQ=255)
   cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGGGGGTTTACCCGCGAagag                                                                                                                                     <  7:252314/149‑1 (MQ=255)
   cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                      <  5:68790/148‑1 (MQ=255)
    acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAgg                                                                                                                                    >  3:276489/1‑149 (MQ=255)
       ccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                                 <  6:144262/149‑1 (MQ=255)
        cacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                                >  7:279744/1‑149 (MQ=255)
         acgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGc                                                                                                                               >  3:173042/1‑149 (MQ=255)
                 cTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGa                                                                                                                       >  8:29122/1‑149 (MQ=255)
                            tcgccTTTCGCGTACTAAGCGAGGAATTACTCGAAGCTGGCATCCCGCATGTCGACATTGTGCGCTACTCGACCGCCTGATACATCACCTGCACTGGCCTCTTCGCGGGTTTAACCGCGAAGATGCCGGCTAAGGGCAAATACATcccc                                                                                                            >  4:362361/5‑149 (MQ=255)
                            tGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGa                                                                                                                                         <  2:93754/120‑1 (MQ=255)
                            tGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGa                                                                                                                                         >  1:93754/1‑120 (MQ=255)
                            tGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCATCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCCATATACATaccc                                                                                                            >  4:295717/1‑149 (MQ=255)
                               gCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCt                                                                                                         <  7:23421/149‑1 (MQ=255)
                                 ttgcgcgTTGCCCGCGAGCAATCCCTCTACGCTGCCACCCCGCTTGTCCACATGGGGCGCGACGCGCGCGCCTGATCCATCACCTCCCCGGGCCTCTTCGCGGGTTACCCGGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGa                                                                                                       <  6:143025/146‑1 (MQ=255)
                                   tCGCGTCGGCAGCGAGGACTTCCCCGACGCGGGCATCCCGCATCTCGACATGGGGCGCGACTCGCGCGCCTGATCCATCACCTGCCCTGCCCGCTTCGCGGGTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGa                                                                                                       <  8:204324/147‑1 (MQ=255)
                                           tCAGCGAGGAATTCCTCGAAGCTGGCATCCCGAATGTCGACATGGGGCGCTACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGTTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCt                                                                                             <  7:18352/149‑1 (MQ=255)
                                            cttcGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCATAGCGGGTTTACCCGCGAAGAGGCCGGCGAAGGCGATATACATCCCCGCTGACAAAAAACCTg                                                                                            >  8:132230/4‑149 (MQ=255)
                                             cgCGAGGAATTCCTCTACGGTGGCACCCCGCATGTCTACATGGTGCGCGCCTCGCGCGCCTGACCCATCACCTCCACTGGCCTCTTCGCGGTTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGt                                                                                           <  4:344188/148‑1 (MQ=255)
                                             aGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTg                                                                                            >  2:11323/1‑148 (MQ=255)
                                              gctagtaaTTCATCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGAACTGGCCTCTTCGCGGTTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                          <  3:295717/143‑1 (MQ=255)
                                              gCGAGTAATTCCTCCAAGCTGCCATCCCCCATGTCCACATGGGCCGCGACTCGCGCGCCTGAGCCATCCCCTGCACTGCCCTCTTGGCGGGTTTACCCGCTATGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                          <  2:81041/149‑1 (MQ=255)
                                               cGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                          >  1:145002/1‑148 (MQ=255)
                                                         cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTCGCCTCTTCGCGGGATTACCCGGGAAGAGGCCGGCTCACGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                  >  7:120278/1‑146 (MQ=255)
                                                          tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                  >  3:155310/1‑145 (MQ=255)
                                                          tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTCCGCGGGTAAACCCGCGAagag                                                                                                                                     <  7:27114/94‑1 (MQ=255)
                                                          tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTCCGCGGGTAAACCCGCGAagag                                                                                                                                     >  8:27114/1‑94 (MQ=255)
                                                                gtgtaaTCCCGATTGTCTAATTGGTGCGCGTCTGGGGCGCTTTATCAATCACTTGCACTGGCTTTTTGGGGGTTTTACGGGTTATTTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaaaa                                                                         <  4:173042/143‑1 (MQ=255)
                                                                  ggCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                      >  5:332030/1‑85 (MQ=255)
                                                                  ggCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                      <  6:332030/85‑1 (MQ=255)
                                                                         cGCTTGTGGACATGGTGCGCGACGCGCGCTCCTTATCCATCACCTGCACTGGCCTCTTCGCGGGTTAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTg                                                               <  7:66158/149‑1 (MQ=255)
                                                                          ccaggtggACTTGGGGCGCGACTCGCGCGCCTCATCCATCCCCTGCACGGGCCTCTGCGCGGGTTAACCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTgc                                                              <  8:154694/142‑1 (MQ=255)
                                                                                     tGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTg                                                   <  7:29122/149‑1 (MQ=255)
                                                                                            cgACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATc                                             >  8:93296/1‑148 (MQ=255)
                                                                                              aCTCGCGCGCCTTATCCATCCCCTCCACGGGCCTCTTGGCGGTTTTACCGGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTa                                          <  4:155310/149‑1 (MQ=255)
                                                                                               cTCGCGCGCCTTATCCACCACCTCCACGGCCCTCTTCGCGGCTTACCCCGCTATGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTAt                                         <  3:6300/149‑1 (MQ=255)
                                                                                                    gcgcCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                  <  2:66256/103‑1 (MQ=255)
                                                                                                    gcgcCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                  >  1:66256/1‑103 (MQ=255)
                                                                                                           aTCCATCACCTCCACTGGCCTCTTCGCGGGTTTACCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCcaccac                               <  2:145002/147‑1 (MQ=255)
                                                                                                             ccATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAcccc                            >  2:341006/1‑148 (MQ=255)
                                                                                                                            gcatctTCGCGGGTATACACGCGAAGAGGCAGGCTCAGGCGATATACATCCCCGCAGACAAAAAACCAGAACATCGAACCAGATAAAACTAGACTCTGCGCCCACGATTGCGAATCCTAAGGCACACCAACCCGATGAAGCGTCCCCgg            >  8:291164/4‑149 (MQ=255)
                                                                                                                            gCCTGTTCGCGGGTTTACCGGCTAAGAGGCCGGCTCAGGCGATATACATCCACGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCgg            <  6:76951/149‑1 (MQ=255)
                                                                                                                               tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                  <  6:199962/108‑1 (MQ=255)
                                                                                                                               tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                  >  5:199962/1‑108 (MQ=255)
                                                                                                                                     cgGGGTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATATAACTTGCCTCTGCGGCCCCGCTTGCGCATCCTATGGCCCCCCACCCCGATGAAGCGTCCCCGGCCATGCGcc   >  6:237085/1‑149 (MQ=255)
                                                                                                                                      gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCt  >  1:188554/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                               
GTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCT  >  NC_002947/4529887‑4530169

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