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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,174,381 T→G C270G (TGC→GGC)  yegU → ADP‑ribosylglycohydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,174,3810TG93.3% 91.3 / ‑0.4 30C270G (TGC→GGC) yegUADP‑ribosylglycohydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (17/11);  total (18/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 = 6.80e-01

GTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATG  >  NZ_CP009273/2174297‑2174462
                                                                                    |                                                                                 
gTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGcc                                                                              >  2:211956/1‑90 (MQ=255)
gTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGcc                                                                              >  1:388571/1‑90 (MQ=255)
 tCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCATAAACCGACCCGAATCGCGGCGCCg                                                                             <  1:15033/90‑1 (MQ=255)
 tCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCg                                                                             <  2:90156/90‑1 (MQ=255)
   gTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTc                                                                           >  2:61313/1‑90 (MQ=255)
         gcAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTgc                                                                     >  2:458393/1‑90 (MQ=255)
                 cAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCt                                                             >  2:637631/1‑90 (MQ=255)
                           gAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGac                                                   <  2:388571/90‑1 (MQ=255)
                             gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                 >  1:485809/1‑90 (MQ=255)
                             gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                 >  1:658681/1‑90 (MQ=255)
                             gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                 >  2:60351/1‑90 (MQ=255)
                             gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                 >  2:373222/1‑90 (MQ=255)
                                cgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGa                                              >  2:625047/1‑90 (MQ=255)
                                cgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGa                                              >  2:289575/1‑90 (MQ=255)
                                     cgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCa                                         <  1:61313/90‑1 (MQ=255)
                                      gtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCAt                                        >  2:497607/1‑90 (MQ=255)
                                       tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                       <  1:637631/90‑1 (MQ=255)
                                       tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                       <  1:142458/90‑1 (MQ=255)
                                         cgcCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCgg                                     <  2:523756/90‑1 (MQ=255)
                                                gcgcTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGgcggcg                                                     >  2:510859/1‑67 (MQ=255)
                                                gcgcTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGgcggcg                                                     <  1:510859/67‑1 (MQ=255)
                                                               gCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                       <  1:16479/66‑1 (MQ=255)
                                                               gCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                       >  2:16479/1‑66 (MQ=255)
                                                               gCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgc               >  2:229001/1‑90 (MQ=255)
                                                               gCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgc               >  2:477848/1‑90 (MQ=255)
                                                                 acaAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgcgg             <  2:187037/90‑1 (MQ=255)
                                                                    aaCCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGgcgc          >  2:18726/1‑90 (MQ=255)
                                                                      ccGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGt        <  1:625047/90‑1 (MQ=255)
                                                                        gACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTg      >  1:546335/1‑90 (MQ=255)
                                                                            cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg  <  1:152335/90‑1 (MQ=255)
                                                                                    |                                                                                 
GTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATG  >  NZ_CP009273/2174297‑2174462

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGCTATCGA  >  NZ_CP009273/2174286‑2174478
                                                                                               |                                                                                                 
ccgagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacaGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGC                                                                                               <  SRR3722112.660922/47‑1 (MQ=60)
 GCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCC                                                                                              >  SRR3722112.393635/1‑100 (MQ=60)
            TCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCATAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGC                                                                                   <  SRR3722112.15225/100‑1 (MQ=60)
                              GCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACAC                                                                 >  SRR3722112.492802/1‑100 (MQ=60)
                              GCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACAC                                                                 >  SRR3722112.668961/1‑100 (MQ=60)
                                    gtataagagacagGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACAC                                                           <  SRR3722112.518322/87‑1 (MQ=60)
                                                CGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTAT                                               <  SRR3722112.62052/100‑1 (MQ=60)
                                                  tgtataagagacagGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG                                             <  SRR3722112.16686/86‑1 (MQ=60)
                                                  tcgtcggcagcgtcagatgtgtataagagacagGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG                                             <  SRR3722112.27060/67‑1 (MQ=60)
                                                  TGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG                                             <  SRR3722112.144242/100‑1 (MQ=60)
                                                  TGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG                                             <  SRR3722112.647538/100‑1 (MQ=60)
                                                                         GGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTG                      >  SRR3722112.554553/1‑100 (MQ=60)
                                                                                 CCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGT              <  SRR3722112.634697/100‑1 (MQ=60)
                                                                                       CGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGC        <  SRR3722112.154214/100‑1 (MQ=60)
                                                                                             GCGGCGCCGGCCTGTGCGCTAACCTTGGCGGGGACACCGGCACCATCGGTGCTATGGCGACGTCAATTTGCGGCGCGTTGCATGGCGTTAACGCTATCGA  <  SRR3722112.292928/100‑1 (MQ=60)
                                                                                               |                                                                                                 
AGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGCTATCGA  >  NZ_CP009273/2174286‑2174478

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