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 I227 R1 223 21.3 1167702 97.0% 1132670 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,064,259 T→G V93G (GTG→GGG)  BW25113_RS10440 → GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,064,2590TG95.8% 69.8 / ‑3.3 24V93G (GTG→GGG) BW25113_RS10440GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (11/12);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCAC  >  NZ_CP009273/2064172‑2064345
                                                                                       |                                                                                      
tACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGc                                                                                      <  1:223004/90‑1 (MQ=255)
    gtgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGtt                                                                                  >  1:244394/1‑90 (MQ=255)
    gtgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGtt                                                                                  >  1:208945/1‑90 (MQ=255)
                    gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                                  >  2:137920/1‑90 (MQ=255)
                     aTGACCGTGCCCTGTATGTGGATGAGAATTTCTGGCGACAAATCCTTCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcca                                                                 >  1:482150/1‑89 (MQ=255)
                              cccTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACgg                                                        <  2:459247/90‑1 (MQ=255)
                                 tgtctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGc                                                     <  2:208945/90‑1 (MQ=255)
                                        ggATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCa                                              <  2:244394/90‑1 (MQ=255)
                                            gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                          >  2:460913/1‑90 (MQ=255)
                                            gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                          >  1:47268/1‑90 (MQ=255)
                                              gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATg                                        >  1:262355/1‑90 (MQ=255)
                                                   tCTGGCGACACATCCTGCAGTGCGGACCTCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGata                                   <  2:137362/90‑1 (MQ=255)
                                                        cGACACATCCTGCANTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGcc                              <  2:361299/90‑1 (MQ=255)
                                                        cGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGcc                              <  2:403737/90‑1 (MQ=255)
                                                               tCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTc                       <  1:365872/90‑1 (MQ=255)
                                                                    cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt                  <  2:61575/90‑1 (MQ=255)
                                                                      gTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGAc                                                               >  2:258023/1‑43 (MQ=255)
                                                                      gTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGAc                                                               <  1:258023/43‑1 (MQ=255)
                                                                        gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                                  <  1:538691/38‑1 (MQ=255)
                                                                        gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                                  >  2:538691/1‑38 (MQ=255)
                                                                        gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTtctc              >  1:543675/1‑90 (MQ=255)
                                                                        gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTtctc              >  1:434236/1‑90 (MQ=255)
                                                                               tCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcaca       <  2:371802/90‑1 (MQ=255)
                                                                                    agGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcacaggcac  <  2:169290/90‑1 (MQ=255)
                                                                                       |                                                                                      
TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCAC  >  NZ_CP009273/2064172‑2064345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 29 ≤ 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
ACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATA  >  NZ_CP009273/2064166‑2064351
                                                                                             |                                                                                            
ACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTT                                                                                        >  SRR3722111.211507/1‑100 (MQ=60)
ACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTT                                                                                        >  SRR3722111.247289/1‑100 (MQ=60)
 tcgtcggcagcgtcagatgtgtataagagacaGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTT                                                                                       <  SRR3722111.13766/68‑1 (MQ=60)
      TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGT                                                                                  <  SRR3722111.225705/100‑1 (MQ=60)
                 ATTAAAGCCGATGACCGTGCCCTGTATGTGGATGAGAATTTCTGGCGACAAATCCTTCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCA                                                                       >  SRR3722111.489542/1‑100 (MQ=60)
                          cattgacattcgtcggcagcgtcagatgtgtataagagacagATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGG                                                              <  SRR3722111.547155/58‑1 (MQ=60)
                             acattcgtcggcagcgtcagatgtgtataagagacagACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGC                                                           <  SRR3722111.261060/63‑1 (MQ=60)
                                        GTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAA                                                >  SRR3722111.47846/1‑100 (MQ=60)
                                          CTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATG                                              >  SRR3722111.265455/1‑100 (MQ=60)
                                                                    ATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTC                    >  SRR3722111.440727/1‑100 (MQ=60)
                                                                    ATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTC                    >  SRR3722111.552226/1‑100 (MQ=60)
                                                                     TCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCC                   <  SRR3722111.371014/100‑1 (MQ=60)
                                                                                      CAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATA  >  SRR3722111.81929/1‑100 (MQ=60)
                                                                                             |                                                                                            
ACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATA  >  NZ_CP009273/2064166‑2064351

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