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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,717,498 G→C G374A (GGC→GCC)  ydhK → FUSC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,717,4980GC96.3% 89.1 / ‑3.0 27G374A (GGC→GCC) ydhKFUSC family protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (14/12);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACT  >  NZ_CP009273/1717414‑1717575
                                                                                    |                                                                             
ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                          <  1:148546/90‑1 (MQ=255)
ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                          >  2:81664/1‑90 (MQ=255)
ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                          <  1:81664/90‑1 (MQ=255)
ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                          <  1:56807/90‑1 (MQ=255)
ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                          <  1:470037/90‑1 (MQ=255)
ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                          <  1:438681/90‑1 (MQ=255)
  aCATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATgg                                                                        <  2:335823/90‑1 (MQ=255)
         gCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTg                                                                 >  1:181462/1‑90 (MQ=255)
                tCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCa                                                          <  2:211081/90‑1 (MQ=255)
                      tACCGATAAGGCCGAAGCTATGTGGGGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGa                                                    <  1:425105/90‑1 (MQ=255)
                           aTAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCg                                               >  2:46296/1‑90 (MQ=255)
                            tAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCgg                                              >  2:636019/1‑90 (MQ=255)
                                  cGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCaa                                        <  1:127147/90‑1 (MQ=255)
                                   gAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAAt                                       >  2:12889/1‑90 (MQ=255)
                                         aTGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTa                                 <  1:587921/90‑1 (MQ=255)
                                                gCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTgg                          >  2:70134/1‑90 (MQ=255)
                                                gCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTgg                          >  1:704936/1‑90 (MQ=255)
                                                       gCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCt                            <  1:507483/81‑1 (MQ=255)
                                                       gCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCt                            >  2:507483/1‑81 (MQ=255)
                                                       gCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAt                   >  1:356387/1‑90 (MQ=255)
                                                          tACATTTTGTACGTTGATGATGATCGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAg                <  1:763111/90‑1 (MQ=255)
                                                           aCATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGc               <  2:618174/90‑1 (MQ=255)
                                                                      gTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactcta    >  2:258420/1‑90 (MQ=255)
                                                                      gTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactcta    >  2:144048/1‑90 (MQ=255)
                                                                      gTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactcta    >  1:634645/1‑90 (MQ=255)
                                                                      gTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactcta    >  1:313438/1‑90 (MQ=255)
                                                                        tgatgatgaTTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactctact  >  2:691163/1‑90 (MQ=255)
                                                                                    |                                                                             
GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACT  >  NZ_CP009273/1717414‑1717575

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ 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
CGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTT  >  NZ_CP009273/1717413‑1717595
                                                                                     |                                                                                                 
CGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTG                                                                                     >  SRR3722073.184967/1‑100 (MQ=60)
 GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGC                                                                                    <  SRR3722073.151598/100‑1 (MQ=60)
 GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGC                                                                                    <  SRR3722073.448060/100‑1 (MQ=60)
 GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGC                                                                                    <  SRR3722073.480048/100‑1 (MQ=60)
 GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGC                                                                                    <  SRR3722073.58149/100‑1 (MQ=60)
 GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGC                                                                                    <  SRR3722073.83509/100‑1 (MQ=60)
                       TACCGATAAGGCCGAAGCTATGTGGGGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCC                                                              <  SRR3722073.434201/100‑1 (MQ=60)
                                   CGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACG                                                  <  SRR3722073.129823/100‑1 (MQ=60)
                                       GCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGG                                              >  SRR3722073.719208/1‑100 (MQ=60)
                                          ATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAG                                           <  SRR3722073.600204/100‑1 (MQ=60)
                                              GGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAT                                       >  SRR3722073.363925/1‑100 (MQ=60)
                                               GAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATT                                      <  SRR3722073.518189/100‑1 (MQ=60)
                                                           TACATTTTGTACGTTGATGATGATCGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTC                          <  SRR3722073.778519/100‑1 (MQ=60)
                                                             CATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTA                        >  SRR3722073.320084/1‑100 (MQ=60)
                                                             CATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTA                        >  SRR3722073.647761/1‑100 (MQ=60)
                                                                                   TGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTT  <  SRR3722073.655853/100‑1 (MQ=60)
                                                                                     |                                                                                                 
CGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTT  >  NZ_CP009273/1717413‑1717595

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