Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,464,055 A→T 22.4% intergenic (+62/‑180) pmbA → / → cybC predicted peptidase required for the maturation and secretion of the antibiotic peptide MccB17/ECK4231:JW4195:b4236; cytochrome b562, truncated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,464,0550AT22.4% 97.0 / 25.1 49intergenic (+62/‑180)pmbA/cybCpredicted peptidase required for the maturation and secretion of the antibiotic peptide MccB17/ECK4231:JW4195:b4236; cytochrome b562, truncated
Reads supporting (aligned to +/- strand):  ref base A (22/16);  new base T (6/5);  total (28/21)
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

GGACAGTAAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAG  >  W3110S.gb/4463985‑4464123
                                                                      |                                                                    
ggACAGTAAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACGCAAAATCATTCGCACTGCa                                                                      >  1:3482955/1‑71 (MQ=255)
ggACAGTAAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCa                                                                      <  1:1069378/71‑1 (MQ=255)
   cAGTAAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCTTc                                                                    >  1:1805779/1‑70 (MQ=255)
   cAGTAAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATc                                                                    >  1:2535018/1‑70 (MQ=255)
   cAGTAAAAAAACTGGCGCGACATGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATc                                                                    >  1:3030306/1‑70 (MQ=255)
     gTAAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATc                                                                    >  1:67818/1‑68 (MQ=255)
      tAAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCTTCgaga                                                                >  1:3029174/1‑71 (MQ=255)
      tAAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCgaga                                                                >  1:391410/1‑71 (MQ=255)
          aaaaCTGGCGCGACCTGCCACGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGa                                                            <  1:2706649/71‑1 (MQ=255)
            aaCTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGa                                                            >  1:799823/1‑69 (MQ=255)
             aCTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACg                                                              <  1:3460029/66‑1 (MQ=255)
                     gACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAg                                                 >  1:3364580/1‑71 (MQ=255)
                        cTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTg                                                       <  1:2214918/62‑1 (MQ=255)
                          gCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCg                                            >  1:4437136/1‑71 (MQ=255)
                                cTATAAACTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGc                                      >  1:1338068/1‑71 (MQ=255)
                                  ataAACTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCg                                     >  1:801282/1‑70 (MQ=255)
                                  ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGaa                                                     >  1:3471404/1‑54 (MQ=255)
                                  ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGc                                      >  1:2939152/1‑69 (MQ=255)
                                  ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCg                                     >  1:4079410/1‑70 (MQ=255)
                                   taAACTAAAAAATTACACAAAATCATTCGGACTGCTTCgaga                                                                >  1:336580/1‑42 (MQ=38)
                                   taAACTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCg                                     <  1:3147579/69‑1 (MQ=255)
                                   taAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGa                                                            <  1:4535935/46‑1 (MQ=255)
                                   taAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTg                                                       <  1:2695392/51‑1 (MQ=255)
                                    aaaCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTAc                                  >  1:2612859/1‑71 (MQ=255)
                                     aaCTAAAAAATTACACAAAATCATTCGCACTGCTTCga                                                                  >  1:54912/1‑38 (MQ=255)
                                     aaCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGaa                                                     >  1:787630/1‑51 (MQ=255)
                                     aaCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGaa                                                     >  1:3413300/1‑51 (MQ=255)
                                     aaCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGaa                                                     >  1:3089511/1‑51 (MQ=255)
                                        tAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCg                               >  1:3305593/1‑70 (MQ=255)
                                         aaaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGaa                                                     >  1:4459808/1‑47 (MQ=255)
                                          aaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGaa                                                           <  1:605858/40‑1 (MQ=255)
                                          aaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAAt                                                    >  1:557630/1‑47 (MQ=255)
                                          aaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGc                              <  1:3174645/69‑1 (MQ=255)
                                           aaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGt                           <  1:3328805/71‑1 (MQ=255)
                                           aaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGt                           <  1:3841443/71‑1 (MQ=255)
                                            aaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGa                                                      >  1:2440883/1‑43 (MQ=255)
                                               ttACACAAAATCATTCTCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGt                       <  1:1445216/71‑1 (MQ=255)
                                               ttACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGt                       <  1:492209/71‑1 (MQ=255)
                                                        aTCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTa                                   >  1:2241478/1‑50 (MQ=255)
                                                        aTCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTc              >  1:1465758/1‑71 (MQ=255)
                                                           aTTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGa            <  1:1718445/70‑1 (MQ=255)
                                                            ttCGCACTGCTTCGAGACGGTAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAAt          <  1:3368836/71‑1 (MQ=255)
                                                            ttCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAAt          <  1:4258942/71‑1 (MQ=255)
                                                            ttCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAAt          <  1:4293236/71‑1 (MQ=255)
                                                            ttCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAAt          <  1:4374478/71‑1 (MQ=255)
                                                               gCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGa            <  1:949597/66‑1 (MQ=255)
                                                                  cTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCg     >  1:3481702/1‑70 (MQ=255)
                                                                   tGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTa                                   <  1:1110574/39‑1 (MQ=255)
                                                                   tGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTa                          <  1:4436322/48‑1 (MQ=255)
                                                                    gCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAg  >  1:860525/1‑71 (MQ=255)
                                                                     cATCGAGACGGAAAGTGAATGAGGGCCGGGATCGTACTCGc                              <  1:932743/41‑1 (MQ=37)
                                                                      |                                                                    
GGACAGTAAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAG  >  W3110S.gb/4463985‑4464123

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