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,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,2590TG100.0% 75.0 / NA 23V93G (GTG→GGG) BW25113_RS10440GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (9/14);  total (9/14)

TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCAC  >  NZ_CP009273/2064172‑2064339
                                                                                       |                                                                                
tACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGc                                                                                <  1:143768/90‑1 (MQ=255)
    gtgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGtt                                                                            >  1:609270/1‑90 (MQ=255)
         tGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTgtgg                                                                       >  1:162188/1‑90 (MQ=255)
         tGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTgtgg                                                                       >  2:392788/1‑90 (MQ=255)
                    gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                            <  1:451209/90‑1 (MQ=255)
                    gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                            <  1:546357/90‑1 (MQ=255)
                           gTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACaaaa                                                     >  2:212005/1‑90 (MQ=255)
                                   tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                             <  2:607351/90‑1 (MQ=255)
                                   tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                             <  1:383278/90‑1 (MQ=255)
                                   tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                             <  1:618938/90‑1 (MQ=255)
                                   tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                             <  2:441696/90‑1 (MQ=255)
                                   tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                             <  2:339431/90‑1 (MQ=255)
                                       tgGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGcc                                         >  1:84602/1‑90 (MQ=255)
                                            gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                    >  2:508575/1‑90 (MQ=255)
                                                aTTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCAt                                       <  2:285541/83‑1 (MQ=255)
                                                aTTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCAt                                       >  1:285541/1‑83 (MQ=255)
                                                aTTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATggg                                >  2:674536/1‑90 (MQ=255)
                                                               tCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTc                 <  2:359315/90‑1 (MQ=255)
                                                               tCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTc                 <  2:626010/90‑1 (MQ=255)
                                                                    cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt            <  2:290735/90‑1 (MQ=255)
                                                                        gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAACGGGATATGGCCGGTATTCAGCCTtctc        >  1:106147/1‑90 (MQ=255)
                                                                         cGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTcc       <  2:288514/90‑1 (MQ=255)
                                                                              aTCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcac  <  2:609270/90‑1 (MQ=255)
                                                                                       |                                                                                
TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCAC  >  NZ_CP009273/2064172‑2064339

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
CTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTC  >  NZ_CP009273/2064162‑2064333
                                                                                                 |                                                                          
CTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGC                                                                          <  SRR3722112.298543/100‑1 (MQ=60)
  TGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTG                                                                        <  SRR3722112.449387/100‑1 (MQ=60)
  TGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTG                                                                        <  SRR3722112.733382/100‑1 (MQ=60)
    ACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTT                                                                      >  SRR3722112.618625/1‑100 (MQ=60)
         GTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGG                                                                 >  SRR3722112.164185/1‑100 (MQ=60)
          TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGT                                                                <  SRR3722112.145559/100‑1 (MQ=60)
                              GATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGG                                            <  SRR3722112.457520/100‑1 (MQ=60)
                              GATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGG                                            <  SRR3722112.554575/100‑1 (MQ=60)
                                       GCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCC                                   >  SRR3722112.85647/1‑100 (MQ=60)
                                             TCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAT                             <  SRR3722112.388258/100‑1 (MQ=60)
                                             TCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAT                             <  SRR3722112.628493/100‑1 (MQ=60)
                                                GTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGG                          >  SRR3722112.288854/1‑100 (MQ=60)
                                                                        ATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAACGGGATATGGCCGGTATTCAGCCTTCTC  >  SRR3722112.107466/1‑100 (MQ=60)
                                                                                                 |                                                                          
CTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTC  >  NZ_CP009273/2064162‑2064333

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