Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,464,055 A→T 25.0% 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,0550AT25.0% 105.8 / 40.0 60intergenic (+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 (26/19);  new base T (9/6);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

AAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAG  >  W3110S.gb/4463992‑4464123
                                                               |                                                                    
aaaaaaaCTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGAc                                                               >  1:2213020/1‑71 (MQ=255)
 aaaaaaCTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACg                                                              >  1:369641/1‑71 (MQ=255)
          gcgcGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTg                                                       <  1:3174909/69‑1 (MQ=255)
          gcgcGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTg                                                       <  1:157788/69‑1 (MQ=255)
          gcgcGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTg                                                       <  1:3297742/69‑1 (MQ=255)
          gcgcGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTg                                                       <  1:33129/69‑1 (MQ=255)
             cgACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAAt                                                    >  1:1821577/1‑69 (MQ=255)
             cgACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAAt                                                    >  1:4051346/1‑69 (MQ=255)
                ccTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCga                                                                  <  1:2255143/52‑1 (MQ=255)
                 cTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGa                                                      <  1:2687178/63‑1 (MQ=255)
                   gCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGa                                                      >  1:3485289/1‑61 (MQ=255)
                      gcgcTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCgggg                                         <  1:3588189/71‑1 (MQ=255)
                        gcTATAAACTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGAAAg                                                         >  1:2415497/1‑53 (MQ=255)
                        gcTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAg                                       <  1:2437494/71‑1 (MQ=255)
                         cTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTg                                                       <  1:606900/54‑1 (MQ=255)
                         cTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTg                                                       <  1:897289/54‑1 (MQ=255)
                         cTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTg                                                       <  1:3440253/54‑1 (MQ=255)
                           ataaaataaaAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAggg                                               >  1:3049661/1‑60 (MQ=255)
                           ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCgaga                                                                >  1:2720589/1‑43 (MQ=255)
                           ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAg                                                 >  1:282117/1‑58 (MQ=255)
                           ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAg                                                 >  1:3046939/1‑58 (MQ=255)
                           ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAg                                                 >  1:3658847/1‑58 (MQ=255)
                           ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAg                                                 >  1:701354/1‑58 (MQ=255)
                           ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCg                                     >  1:2567370/1‑70 (MQ=255)
                           ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCg                                     >  1:2909777/1‑70 (MQ=255)
                            taAACTAAAAAATTACACAAAATCATTCGCACTGCTTc                                                                    <  1:1504320/38‑1 (MQ=255)
                            taAACTAAAAAATTACACAAAATCATTCGCACTGCTTc                                                                    <  1:1489694/38‑1 (MQ=255)
                            taAACTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTg                                                       <  1:1893739/51‑1 (MQ=255)
                            taAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCg                                     <  1:2261040/69‑1 (MQ=255)
                              aaCTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATg                                                   >  1:1138494/1‑53 (MQ=255)
                              aaCTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACt                                 >  1:4039332/1‑71 (MQ=255)
                              aaCTAAAAAATTACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACt                                 >  1:748230/1‑71 (MQ=255)
                              aaCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGaaa                                                          >  1:2893918/1‑46 (MQ=255)
                              aaCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGaaa                                                          >  1:2560824/1‑46 (MQ=255)
                              aaCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCg                                     <  1:938310/67‑1 (MQ=255)
                              aaCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCg                                     <  1:2239805/67‑1 (MQ=255)
                                cTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAg                                                         >  1:968424/1‑45 (MQ=255)
                                cTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAg                                                         >  1:2128206/1‑45 (MQ=255)
                                cTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCg                               >  1:1517579/1‑71 (MQ=255)
                                cTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCg                               >  1:537974/1‑71 (MQ=255)
                                cTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGCAAAGTGAATGAGGGCCGGGGAGCGTACTCg                               >  1:3090922/1‑71 (MQ=255)
                                  aaaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGaa                                                     >  1:3244599/1‑47 (MQ=255)
                                  aaaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGaa                                                     >  1:1564449/1‑47 (MQ=255)
                                  aaaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGa                                                      >  1:1970842/1‑46 (MQ=255)
                                   aaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGaa                                                           <  1:1117179/40‑1 (MQ=255)
                                   aaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGaa                                                           <  1:1845397/40‑1 (MQ=255)
                                   aaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGaa                                                           <  1:3464146/40‑1 (MQ=255)
                                   aaaaaTTACACAAAATCATTAGCACTGCATCGAGACGGaa                                                           <  1:3354295/40‑1 (MQ=255)
                                      aaTTACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTAc                         >  1:965252/1‑71 (MQ=255)
                                       aTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACg                        >  1:2084106/1‑71 (MQ=255)
                                        ttACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGt                       <  1:3911050/71‑1 (MQ=255)
                                        ttACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGt                       >  1:3723892/1‑71 (MQ=255)
                                        ttACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGt                       >  1:2673864/1‑71 (MQ=255)
                                        ttACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGt                       <  1:1530298/71‑1 (MQ=255)
                                        ttACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGt                       >  1:2834764/1‑71 (MQ=255)
                                         tACACAAAATCATTCGCACTGCTTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGt                       >  1:2122799/1‑70 (MQ=255)
                                          acacAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCg                               <  1:2434753/61‑1 (MQ=255)
                                                            tGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGAccc                  >  1:442106/1‑56 (MQ=255)
                                                             gCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGc                              <  1:632652/43‑1 (MQ=255)
                                                              cATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAg  >  1:2312653/1‑70 (MQ=255)
                                                               |                                                                    
AAAAAAACTGGCGCGACCTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAG  >  W3110S.gb/4463992‑4464123

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: