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

CAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCG  >  NC_002947/4529884‑4530166
                                                                                                                                              |                                                                                                                                                
cAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                           >  2:112865/1‑149 (MQ=255)
 aGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCcgc                                                                                                                                           >  4:304889/1‑148 (MQ=255)
    tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAag                                                                                                                                       <  1:351494/149‑1 (MQ=255)
      cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAagag                                                                                                                                    <  2:95691/149‑1 (MQ=255)
      cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAagag                                                                                                                                    <  7:359051/149‑1 (MQ=255)
      cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAagag                                                                                                                                    <  3:304889/149‑1 (MQ=255)
      cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACCGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                      <  6:358109/148‑1 (MQ=255)
           cacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                               >  8:212971/1‑149 (MQ=255)
               ccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCGGCt                                                                                                                             >  5:388736/1‑147 (MQ=255)
                   tCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGc                                                                                                                        >  8:335240/1‑148 (MQ=255)
                     tATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcg                                                                                                                                            >  3:139811/1‑127 (MQ=255)
                     tATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcg                                                                                                                                            <  4:139811/127‑1 (MQ=255)
                                    tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                        >  2:54124/1‑116 (MQ=255)
                                    tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                        <  1:54124/116‑1 (MQ=255)
                                              gcttcGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCACCGCTGACAAAAAACCt                                                                                           >  3:28904/5‑149 (MQ=255)
                                                aGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCAGACTCGCAGGGGTTACCCCGAAGACGCCCGGCGCAGGCAAATAAATACCCGCGGAAAAAAAAccctg                                                                                         >  5:233164/1‑146 (MQ=255)
                                                  cGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAc                                                                                      <  5:158690/149‑1 (MQ=255)
                                                   gAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAca                                                                                     >  6:181856/1‑149 (MQ=255)
                                                     ggAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGt                                                                                         >  5:67592/1‑144 (MQ=255)
                                                     ggAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGt                                                                                         <  6:67592/144‑1 (MQ=255)
                                                     ggAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAACTGTAcatc                                                                                   >  4:124848/1‑149 (MQ=255)
                                                        aTTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATCTACATCCCCGCTGACAAAAAACCGGTAcatccat                                                                                >  7:285967/1‑149 (MQ=255)
                                                             tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGTGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                               >  2:74951/1‑145 (MQ=255)
                                                             tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                               >  2:47439/1‑145 (MQ=255)
                                                             tCGAAGCTGGCATCCCGCATGTCGACATGGAGCGCGACTCGCGCGCCTCATCCATCACCTGCACTGGCCTATTCGCGGGTTTACCCGCGAAGAGGAAGGCTCAGGCGATATACATCCCCCCTGACAAAAAACCTGTACAACCAACCAg                                                                            >  7:71092/1‑148 (MQ=255)
                                                              cGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGGTTAACCGCGAAGCGGCCGGCTCAGGCGCATTACATGCCCGCGGACAAAAAAACAGTAGAtacaaccaaa                                                                           >  1:385225/1‑146 (MQ=255)
                                                                 aGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaa                                                                        >  3:193955/1‑148 (MQ=255)
                                                                  gCTGGCATCCCGCATGTCTACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGAACTGGCCTCTTCGCGGTTTAACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaaaa                                                                      <  7:335240/149‑1 (MQ=255)
                                                                  gCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTTATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGTCAAAAAACCTGTACATCCATCCAGATaaaa                                                                      <  1:33236/149‑1 (MQ=255)
                                                                     ggCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGTACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                               >  3:113513/1‑137 (MQ=255)
                                                                        aTCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAACCTGTACATCCATCCAGATAAAACTTGc                                                                 <  1:47439/149‑1 (MQ=255)
                                                                           ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCctct                                                             >  5:118731/1‑149 (MQ=255)
                                                                           ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCctct                                                             >  2:314993/1‑149 (MQ=255)
                                                                           ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATACCCGCGGCAAAAAAACCTGTAAATCCATCCAGATAAAACTTGCctct                                                             >  3:314/1‑149 (MQ=255)
                                                                           ccGCATGTCCACAGGGTCCGCGACTCCCCCGCCTGATCCACCACCTGCACTGCCCTCTTCGCGGGTTCACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCctct                                                             <  8:94322/149‑1 (MQ=255)
                                                                                           tGCGTGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTCACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgcg                                              <  1:74951/149‑1 (MQ=255)
                                                                                            gcgcGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCa                                            >  7:133240/1‑149 (MQ=255)
                                                                                                  cTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTAt                                      <  5:181856/149‑1 (MQ=255)
                                                                                                                      ccTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAg                  <  1:112865/149‑1 (MQ=255)
                                                                                                                             tGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTcccc           >  6:87143/1‑149 (MQ=255)
                                                                                                                               gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCCCCACCCCGATGAAGCGTCCCCgg         >  6:202092/1‑149 (MQ=255)
                                                                                                                               gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCgg         >  4:287878/1‑149 (MQ=255)
                                                                                                                                  tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa      >  2:27236/1‑149 (MQ=255)
                                                                                                                                    ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGccc                                                       <  4:116724/98‑1 (MQ=255)
                                                                                                                                    ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGccc                                                       >  3:116724/1‑98 (MQ=255)
                                                                                                                                    ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATg    >  3:203964/1‑149 (MQ=255)
                                                                                                                                      cgcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGACCACCACCCCGATGAAGCGTCCCCGGCCATgcg  >  2:276028/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCG  >  NC_002947/4529884‑4530166

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: