Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,012,376 G→A 12.7% 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,3760GA12.7% 155.1 / 12.4 55A500T (GCC→ACC) rbbAfused ribosome‑associated ATPases
Reads supporting (aligned to +/- strand):  ref base G (21/27);  new base A (3/4);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CATCTCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAACGAAGCCGAAG  >  W3110S.gb/4012306‑4012445
                                                                      |                                                                     
cATCTCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAg                                                                       <  1:2530825/71‑1 (MQ=255)
     cACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAgccgc                                                                   >  1:3596028/1‑70 (MQ=255)
      aCTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAgcc                                                                     <  1:1971975/67‑1 (MQ=255)
      aCTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCa                                                                 <  1:1713591/71‑1 (MQ=255)
          aTGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGGACTGGTTGAGAAACGCGGAGCCGCCAGTCt                                                             <  1:1462034/71‑1 (MQ=255)
                gCCGGAAAAGTGCTTGCCAGAGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGaa                                                         >  1:422027/1‑69 (MQ=255)
                 ccGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGaa                                                         >  1:3515421/1‑68 (MQ=255)
                   ggAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAgg                                                     >  1:2321812/1‑70 (MQ=255)
                   ggAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGc                                                    <  1:3296945/71‑1 (MQ=255)
                       aaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAttt                                                <  1:2420878/71‑1 (MQ=255)
                       aaGTGCTTGCCAGCGGTACACCGCAGGAACTGGGTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAttt                                                <  1:77182/71‑1 (MQ=255)
                         gTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTAt                                              <  1:3573356/71‑1 (MQ=255)
                          tGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATc                                             <  1:2309158/71‑1 (MQ=255)
                            cTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTg                                                            <  1:1471800/54‑1 (MQ=255)
                               gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGaa                                                         <  1:4063298/54‑1 (MQ=255)
                               gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGcc                                          <  1:2778695/69‑1 (MQ=255)
                               gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCt                                         <  1:2901557/70‑1 (MQ=255)
                               gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTa                                        >  1:3500686/1‑71 (MQ=255)
                               gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTa                                        >  1:2144497/1‑71 (MQ=255)
                               gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTa                                        >  1:2260360/1‑71 (MQ=255)
                                ccAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAaga                                                       >  1:1997309/1‑55 (MQ=255)
                                  aGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAt                                                  <  1:1858668/58‑1 (MQ=255)
                                   gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                                    <  1:569528/71‑1 (MQ=255)
                                   gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                                    >  1:3880029/1‑71 (MQ=255)
                                   gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                                    <  1:1462063/71‑1 (MQ=255)
                                    cGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                                    >  1:878674/1‑70 (MQ=255)
                                       tACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTa                                               <  1:4004429/56‑1 (MQ=255)
                                        acacCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAgg                                <  1:82684/70‑1 (MQ=255)
                                        acacCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAgg                                >  1:3906360/1‑70 (MQ=255)
                                         cacCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAttt                                                <  1:91740/53‑1 (MQ=255)
                                               aGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggca                        <  1:3252179/71‑1 (MQ=255)
                                               aGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggca                        <  1:4107620/71‑1 (MQ=255)
                                                ggAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcag                       >  1:4146149/1‑71 (MQ=255)
                                                ggAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcag                       >  1:4141077/1‑71 (MQ=255)
                                                ggAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcag                       <  1:2940368/71‑1 (MQ=255)
                                                 gAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCgg                          <  1:1450121/67‑1 (MQ=255)
                                                 gAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCgg                          <  1:3169262/67‑1 (MQ=255)
                                                   aCTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                     >  1:2301031/1‑70 (MQ=255)
                                                   aCTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                     >  1:613314/1‑70 (MQ=255)
                                                     tGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                     >  1:1510548/1‑68 (MQ=255)
                                                     tGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                     >  1:2715643/1‑68 (MQ=255)
                                                       gTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGc                                   >  1:1330063/1‑52 (MQ=255)
                                                        ttGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGc               <  1:861158/71‑1 (MQ=255)
                                                         tGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGc               >  1:767856/1‑70 (MQ=255)
                                                         tGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGc               >  1:1013118/1‑70 (MQ=255)
                                                          gagaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGc                                           >  1:2770061/1‑41 (MQ=255)
                                                          gagaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTAttt                                     >  1:3507845/1‑47 (MQ=255)
                                                           agaAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGGAAc            <  1:2718823/71‑1 (MQ=255)
                                                           agaAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAc            <  1:1099475/71‑1 (MQ=255)
                                                            gaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcg           <  1:2108340/71‑1 (MQ=255)
                                                              aaCGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcga          <  1:1601151/70‑1 (MQ=255)
                                                                 gcgGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCt                                         >  1:1062140/1‑36 (MQ=255)
                                                                 gcgGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGa                               >  1:2840885/1‑46 (MQ=255)
                                                                  cgGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccg     <  1:1439531/71‑1 (MQ=255)
                                                                  cgGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccg     <  1:1093912/71‑1 (MQ=255)
                                                                   ggAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccg     <  1:691705/70‑1 (MQ=255)
                                                                   ggAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccga    <  1:3783706/71‑1 (MQ=255)
                                                                     aGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccgaag  >  1:1966307/1‑71 (MQ=255)
                                                                      |                                                                     
CATCTCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAACGAAGCCGAAG  >  W3110S.gb/4012306‑4012445

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