Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,756,138 Δ1 bp intergenic (+183/‑21) ssrA → / → intA tmRNA/CP4‑57 prophage; integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,756,1370C.98.4% 286.9 / ‑1.8 63intergenic (+182/‑22)ssrA/intAtmRNA/CP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (29/33);  total (29/34)
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

AGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAA  >  NC_000913/2756103‑2756165
                                  |                            
aGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                              >  1:6464109/1‑35 (MQ=255)
aGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                              <  1:3028392/35‑1 (MQ=255)
 ggCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                             >  1:4072998/1‑34 (MQ=255)
   cccAACGAAAAGCTCTATTGTTTACGTTGGGCCTa                           <  1:21117/35‑1 (MQ=255)
     cAACGAAAAGCTCTATTGTTTACGTTGGGCTAAAc                        <  1:5554666/35‑1 (MQ=38)
     cAACGAAAAGCTCTATTGTTTACGTTGGGCTAAAc                        <  1:3994589/35‑1 (MQ=38)
      aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACg                       <  1:5702252/35‑1 (MQ=38)
      aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACg                       <  1:6897843/35‑1 (MQ=38)
       aCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGc                      >  1:3532443/1‑35 (MQ=38)
        cGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCa                     >  1:22900/1‑35 (MQ=38)
        cGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCa                     >  1:3440743/1‑35 (MQ=38)
         gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAg                    >  1:7163316/1‑35 (MQ=38)
         gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAg                    <  1:6708253/35‑1 (MQ=38)
          aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAgg                   >  1:3470847/1‑35 (MQ=38)
          aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAgg                   <  1:3458340/35‑1 (MQ=38)
          aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAgg                   >  1:6254666/1‑35 (MQ=38)
          aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAgg                   >  1:5671070/1‑35 (MQ=38)
          aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAgg                   <  1:6655517/35‑1 (MQ=38)
          aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAgg                   >  1:7191728/1‑35 (MQ=38)
           aaaGCTCTATTGTTTACGTTGGGCTAAACGCAggg                  >  1:6641884/1‑35 (MQ=38)
           aaaGCTCTATTGTTTACGTTGGGCTAAACGCAggg                  >  1:2237433/1‑35 (MQ=38)
           aaaGCTCTATTGTTTACGTTGGGCTAAACGCAggg                  >  1:7738162/1‑35 (MQ=38)
           aaaGCTCTATTGTTTACGTTGGGCTAAACGCAggg                  >  1:5594346/1‑35 (MQ=38)
           aaaGCTCTATTGTTTACGTTGGGCTAAACGCAggg                  >  1:820872/1‑35 (MQ=38)
           aaaGCTCCATTGTTTACGTTGGGCTAAACGCAggg                  >  1:6392150/1‑35 (MQ=25)
            aaGCTCTATTGTTTACGTTGGGCTAAACGCAGGga                 >  1:5821559/1‑35 (MQ=38)
            aaGCTCTATTGTTTACGTTGGGCTAAACGCAGGga                 <  1:7334261/35‑1 (MQ=38)
            aaGCTCTATTGTTTACGTTGGGCTAAACGCAGGga                 >  1:7437888/1‑35 (MQ=38)
             aGCTCTATTGTTTACGTTGGGCTAAACGCAGGgag                >  1:539053/1‑35 (MQ=38)
             aGCTCTATTGTTTACGTTGGGCTAAACGCAGGgag                >  1:5662128/1‑35 (MQ=38)
             aGCTCTATTGTTTACGTTGGGCTAAACGCAGGgag                >  1:4468827/1‑35 (MQ=38)
              gCTCTATTGTTTACGTTGGGCTAAACGCAGGgaga               <  1:2662085/35‑1 (MQ=38)
              gCTCTATTGTTTACGTTGGGCTAAACGCAGGgaga               <  1:4544512/35‑1 (MQ=38)
                 ctATTGTTTACGTTGGGCTAAACGCAGGGAGACTc            >  1:6622273/1‑35 (MQ=38)
                  tATTGTTTACGTTGGGCTAAACGCAGGGAGACTcc           >  1:3600057/1‑35 (MQ=38)
                   aTTGTTTACGTTGGGCTAAACGCAGGGAGACTccc          <  1:6248019/35‑1 (MQ=38)
                    ttGTTTACGTTGGGCTAAACGCAGGGAGACTcccc         <  1:7078055/35‑1 (MQ=38)
                    ttGTTTACGTTGGGCTAAACGCAGGGAGACTcccc         <  1:164904/35‑1 (MQ=38)
                    ttGTTTACGTTGGGCTAAACGCAGGGAGACTcccc         >  1:5472525/1‑35 (MQ=38)
                    ttGTTTACGTTGGGCTAAACGCAGGGAGACTcccc         <  1:3416270/35‑1 (MQ=38)
                    ttGTTTACGTTGGGCTAAACGCAGGGAGACTcccc         <  1:5153791/35‑1 (MQ=38)
                    ttGTTTACGTTGGGCTAAACGCAGGGAGACTcccc         <  1:731672/35‑1 (MQ=38)
                     tGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA        <  1:5038894‑M1/35‑1 (MQ=255)
                     tGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA        <  1:4517438‑M1/35‑1 (MQ=255)
                      gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAc       <  1:5188291‑M1/35‑2 (MQ=255)
                      gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAc       <  1:7889841‑M1/35‑2 (MQ=255)
                      gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAc       <  1:5524957‑M1/35‑2 (MQ=255)
                       tttACGTTGGGCTAAACGCAGTGAGACTCCCCAcc      >  1:2450369‑M1/1‑33 (MQ=255)
                       tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcc      <  1:5900939‑M1/35‑3 (MQ=255)
                       tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcc      <  1:3587960‑M1/35‑3 (MQ=255)
                       tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcc      <  1:3506016‑M1/35‑3 (MQ=255)
                       tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcc      >  1:6803093‑M1/1‑33 (MQ=255)
                        ttACGTTGGGCTAAACGCAGGGAGACTCCCCAcct     >  1:1519758‑M1/1‑32 (MQ=255)
                         tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt    <  1:2095215‑M1/35‑5 (MQ=255)
                         tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt    >  1:7361378‑M1/1‑31 (MQ=255)
                         tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt    >  1:5312253‑M1/1‑31 (MQ=255)
                          aCGTTGGGCTAAACGCAGGGAGACTCCCCAccttt   <  1:7061829‑M1/35‑6 (MQ=255)
                          aCGTTGGGCTAAACGCAGGGAGACTCCCCAccttt   <  1:3465330‑M1/35‑6 (MQ=255)
                          aCGTTGGGCTAAACGCAGGGAGACTCCCCAccttt   <  1:4042105‑M1/35‑6 (MQ=255)
                          aCGTTGGGCTAAACGCAGGGAGACTCCCCAccttt   >  1:5313097‑M1/1‑30 (MQ=255)
                           cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg  >  1:3314158‑M1/1‑29 (MQ=255)
                           cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg  <  1:2397150‑M1/35‑7 (MQ=255)
                           cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg  <  1:6471714‑M1/35‑7 (MQ=255)
                           cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg  <  1:1952363‑M1/35‑7 (MQ=255)
                           cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg  <  1:1814632‑M1/35‑7 (MQ=255)
                           cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg  <  1:5641922‑M1/35‑7 (MQ=255)
                                  |                            
AGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAA  >  NC_000913/2756103‑2756165

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: