Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,012,376 G→A 38.3% A500T (GCC→ACC)  rbbA → fused ribosome‑associated ATPases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,012,3760GA38.3% 39.4 / 49.6 47A500T (GCC→ACC) rbbAfused ribosome‑associated ATPases
Reads supporting (aligned to +/- strand):  ref base G (11/18);  new base A (5/13);  total (16/31)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CATCTCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAACGAAGCCGAAGC  >  W3110S.gb/4012306‑4012446
                                                                      |                                                                      
cATCTCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAg                                                                        <  1:5024293/71‑1 (MQ=255)
    tcACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAgccgc                                                                    >  1:255265/1‑71 (MQ=255)
    tcACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAaccgc                                                                    >  1:334467/1‑71 (MQ=255)
        tGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAg                                                                 <  1:2032104/70‑1 (MQ=255)
          aTGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGt                                                                <  1:3066322/69‑1 (MQ=255)
           tGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTg                                                             >  1:1274045/1‑71 (MQ=255)
             cACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAaccgc                                                                    >  1:3497927/1‑62 (MQ=255)
              aCGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGaa                                                          <  1:1135382/71‑1 (MQ=255)
               cGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTgg                                                            >  1:3114172/1‑68 (MQ=255)
                 ccGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAagag                                                       >  1:2094731/1‑71 (MQ=255)
                 ccGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAagag                                                       >  1:2100638/1‑71 (MQ=255)
                  cGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTgg                                                            <  1:1037934/65‑1 (MQ=255)
                    gAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCa                                                    >  1:2614174/1‑71 (MQ=255)
                       aaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAttt                                                 <  1:4728771/71‑1 (MQ=255)
                        aGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTa                                                <  1:4297851/71‑1 (MQ=255)
                        aGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTa                                                <  1:3160320/71‑1 (MQ=255)
                        aGTGCTTGCCAGCGGTACACAGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCa                                                    <  1:4967222/67‑1 (MQ=255)
                        aGTGCTTCCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTa                                                <  1:4608984/71‑1 (MQ=255)
                            cTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAtt                                                  <  1:1054447/65‑1 (MQ=255)
                            cTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGc                                            >  1:3880827/1‑71 (MQ=255)
                               gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCt                                          <  1:1681853/70‑1 (MQ=255)
                                ccAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTa                                         <  1:3867258/70‑1 (MQ=255)
                                  aGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTAttt                                      <  1:1603358/71‑1 (MQ=255)
                                   gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAaga                                                        <  1:2617626/52‑1 (MQ=255)
                                   gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAaga                                                        <  1:676506/52‑1 (MQ=255)
                                   gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAaga                                                        <  1:1406316/52‑1 (MQ=255)
                                    cGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGc                                    <  1:1546306/71‑1 (MQ=255)
                                    cGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGc                                    <  1:3857376/71‑1 (MQ=255)
                                    cGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGc                                    <  1:4763457/71‑1 (MQ=255)
                                     ggTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAttt                                                 <  1:4068386/57‑1 (MQ=255)
                                     ggTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCa                                   <  1:1707215/71‑1 (MQ=255)
                                      gTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAg                                  <  1:3919409/71‑1 (MQ=255)
                                       tACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCAtt                                                  <  1:2561131/54‑1 (MQ=255)
                                         cacCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAgg                                 >  1:3050060/1‑69 (MQ=255)
                                            cGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAGAGGCATTTATCGCCTATTTGCAGGAAGCgg                           >  1:3995486/1‑71 (MQ=255)
                                             gCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTAtt                                       >  1:3105268/1‑59 (MQ=255)
                                               aGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCt                                          <  1:3276122/54‑1 (MQ=255)
                                                  aaCTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                      >  1:1289206/1‑71 (MQ=255)
                                                   aCTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAaga                                                        >  1:3608028/1‑36 (MQ=255)
                                                    cTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                      >  1:3310428/1‑69 (MQ=255)
                                                     tGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTAttt                                      <  1:2139214/52‑1 (MQ=255)
                                                       gTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCa                                                    <  1:640100/36‑1 (MQ=255)
                                                           agaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAc             <  1:4635623/71‑1 (MQ=255)
                                                            gaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAc             <  1:855202/70‑1 (MQ=255)
                                                            gaAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcg            <  1:124894/71‑1 (MQ=255)
                                                             aaaCGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcg            >  1:3127693/1‑70 (MQ=255)
                                                                cgcgGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggca                         >  1:181098/1‑54 (MQ=255)
                                                                cgcgGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggca                         >  1:1522819/1‑54 (MQ=255)
                                                                cgcgGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggca                         >  1:638170/1‑54 (MQ=255)
                                                                cgcgGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAc             <  1:3024142/66‑1 (MQ=255)
                                                                   gtaGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGc                             <  1:3095655/45‑1 (MQ=255)
                                                                      gccgccAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccgaagc  >  1:1852442/1‑71 (MQ=255)
                                                                      |                                                                      
CATCTCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAACGAAGCCGAAGC  >  W3110S.gb/4012306‑4012446

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