Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F2 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,623,127 (T)5→4 intergenic (+287/‑24) yfgG → / → yfgH YfgG family protein/glycine zipper 2TM domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,623,1230T.92.3% 84.7 / 3.0 26intergenic (+283/‑28)yfgG/yfgHYfgG family protein/glycine zipper 2TM domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (13/11);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

TGCGTGAATTTTAATAAATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTC  >  NZ_CP009273/2623042‑2623205
                                                                                 |                                                                                  
tGCGTGAATTTTAATAAATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           <  1:142411/90‑1 (MQ=255)
                          aTCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGaa                                                 >  1:294495/1‑90 (MQ=255)
                               tttATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           >  1:622800/1‑59 (MQ=255)
                               tttATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           <  2:622800/59‑1 (MQ=255)
                                 tataCGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTaaa                                          <  1:167386/90‑1 (MQ=255)
                                  ataCGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTaaaa                                         >  2:883975/1‑90 (MQ=255)
                                       cAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATg                                     >  2:615270/1‑90 (MQ=255)
                                       cAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGt                                    >  1:396746/1‑90 (MQ=255)
                                        gaTACATTTAGTGTCCTCTTTTGGTGATGTTAAATGTCTTATTTTCGTAATGGGTATAACAAGGAATATGGAGAAAATTTAAAAAGTGtc                                   >  1:733103/2‑90 (MQ=255)
                                         aTACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGtct                                  >  1:415470/1‑90 (MQ=255)
                                           aCATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGtcttc                                <  1:130747/90‑1 (MQ=255)
                                            cATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGtcttc                                <  1:615270/90‑1 (MQ=255)
                                                 aCTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATaa                                                                  >  1:759454/1‑50 (MQ=255)
                                                 aCTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATaa                                                                  <  2:759454/50‑1 (MQ=255)
                                                 aCTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCtg                          >  1:670708/1‑90 (MQ=255)
                                                   tttCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCtgtg                        >  2:229897/1‑90 (MQ=255)
                                                    ttCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           >  1:915653/1‑38 (MQ=38)
                                                    ttCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           <  2:915653/38‑1 (MQ=38)
                                                     tCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAAt                                                                          >  1:696878/1‑38 (MQ=38)
                                                     tCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAAt                                                                          <  2:696878/38‑1 (MQ=38)
                                                      ccTCTTTTGATGATCTTAAATGTCTTATTTTCGTAAtg                                                                         >  2:380195/1‑38 (MQ=38)
                                                      ccTCTTTTGATGATCTTAAATGTCTTATTTTCGTAAtg                                                                         <  1:380195/38‑1 (MQ=38)
                                                                   tCTTAAATGTCTTATTTTCGTAATGTGTATAACaa                                                               >  2:43134/1‑35 (MQ=255)
                                                                   tCTTAAATGTCTTATTTTCGTAATGTGTATAACaa                                                               <  1:43134/35‑1 (MQ=255)
                                                                     ttAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTc      <  1:883975/90‑1 (MQ=255)
                                                                         aTGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCActc  <  2:586528/90‑1 (MQ=255)
                                                                                 |                                                                                  
TGCGTGAATTTTAATAAATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTC  >  NZ_CP009273/2623042‑2623205

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGCGTGAATTTTAATAAATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTCTGGCAGGATGCCAGTCAAAT  >  NZ_CP009273/2623042‑2623225
                                                                                     |                                                                                                  
TGCGTGAATTTTAATAAATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAAC                                                                                     <  SRR3722072.144958/100‑1 (MQ=60)
  aatcggacgtctcgtgggctcggagatgtgtataagagacagCATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAA                                                                                   <  SRR3722072.386988/58‑1 (MQ=60)
            gataatcggacgtctcgtgggctcggagatgtgtataagagacagCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGA                                                                         <  SRR3722072.43925/55‑1 (MQ=60)
                AATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAA                                                                     >  SRR3722072.299589/1‑100 (MQ=60)
                     AATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACctgtctcttatacacatctcc                                                                >  SRR3722072.633914/1‑79 (MQ=60)
                             CCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGT                                                        >  SRR3722072.403905/1‑100 (MQ=60)
                               TTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCT                                                      >  SRR3722072.422991/1‑100 (MQ=60)
                                 TATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTC                                                    <  SRR3722072.170280/100‑1 (MQ=60)
                                       CAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGctgtctcttatacacatctccgagcccacg                                              >  SRR3722072.772843/1‑70 (MQ=60)
                                       CAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTG                                              >  SRR3722072.682715/1‑100 (MQ=60)
                                          TACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACctgcctcttatacacatctccgagcccacgagacgtccgatt                                           >  SRR3722072.931391/1‑58 (MQ=60)
                                           ACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACActgtctcttatacacatctccgagcccacgagacgtccgatt                                          >  SRR3722072.709346/1‑58 (MQ=60)
                                           ACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGC                                          <  SRR3722072.133097/100‑1 (MQ=60)
                                            CATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGC                                          <  SRR3722072.626237/100‑1 (MQ=60)
                                                                   TCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTCTGGC                  >  SRR3722072.155112/1‑100 (MQ=60)
                                                                   TCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTCTGGC                  >  SRR3722072.903693/1‑100 (MQ=60)
                                                                    CTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTCTGGCA                 >  SRR3722072.2599/1‑100 (MQ=60)
                                                                     TTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTCTGGCAG                <  SRR3722072.899192/100‑1 (MQ=60)
                                                                           GTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGActgtctcttatacacatctccgagcccacgagacgtccgattatctcgtatgccgtcttctgc          >  SRR3722072.561160/1‑37 (MQ=60)
                                                                                    TTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTCTGGCAGGATGCCAGTCAAAT  >  SRR3722072.386542/1‑100 (MQ=60)
                                                                                     |                                                                                                  
TGCGTGAATTTTAATAAATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTCTGGCAGGATGCCAGTCAAAT  >  NZ_CP009273/2623042‑2623225

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: