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 4,153,816 T→C S84P (TCA→CCA)  btuB → TonB‑dependent vitamin B12 receptor BtuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,153,8160TC100.0% 65.6 / NA 20S84P (TCA→CCA) btuBTonB‑dependent vitamin B12 receptor BtuB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/11);  total (9/11)

GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153729‑4153901
                                                                                       |                                                                                     
gACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCa                                                                                     <  1:521715/90‑1 (MQ=255)
                    ggTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTAc                                                                 <  1:106925/90‑1 (MQ=255)
                      tCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACaa                                                               <  1:433463/90‑1 (MQ=255)
                          tgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGc                                                           <  2:279509/90‑1 (MQ=255)
                           gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                          <  2:782004/90‑1 (MQ=255)
                           gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                          <  1:769288/90‑1 (MQ=255)
                             tgtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAg                                                        <  1:407309/90‑1 (MQ=255)
                                  tgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCCTCTTTTTTTATTCGCGGTACAAAAGCCCGTCAtg                                                   >  1:274282/1‑90 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAAAGCCAGTCAttt                                                  >  2:832273/1‑88 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGGACAAGTGCCAGTCAtgt                                                  >  1:913024/1‑90 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCCGCTCCCCTCTTTTTTTTTTCGCGGGACAAATGCCCGTCAtgt                                                  >  2:801656/1‑90 (MQ=255)
                                       cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                              <  2:2363/90‑1 (MQ=255)
                                                   gTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTCGCGGGACAAATGCCAGCCATGTTTTGGGGTTAATTGAt                                  >  1:52342/1‑90 (MQ=255)
                                                   gTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCCGTCATGTGTTGGGGTTAATTGAt                                  >  1:557101/1‑90 (MQ=255)
                                                         atCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTa                            <  2:199685/90‑1 (MQ=255)
                                                                 aaaCGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGGGTTAATTGATGGCGTACGCCTGaa                    >  1:401586/1‑90 (MQ=255)
                                                                       cggTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGc              <  2:605078/90‑1 (MQ=255)
                                                                                   gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGt                                                       <  2:20636/37‑1 (MQ=38)
                                                                                   gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGt                                                       >  1:20636/1‑37 (MQ=38)
                                                                                   gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTgg  >  1:184459/1‑90 (MQ=255)
                                                                                       |                                                                                     
GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153729‑4153901

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGAC  >  NZ_CP009273/4153721‑4153911
                                                                                               |                                                                                               
AGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATC                                                                                             <  SRR3722072.911428/100‑1 (MQ=60)
     ATCGACCGCTAGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTT                                                                                        <  SRR3722072.509493/100‑1 (MQ=60)
      TCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCATCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTT                                                                                       <  SRR3722072.844024/100‑1 (MQ=60)
        GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTA                                                                                     <  SRR3722072.531103/100‑1 (MQ=60)
                            GGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGT                                                                 <  SRR3722072.108897/100‑1 (MQ=60)
                              TCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCA                                                               <  SRR3722072.441352/100‑1 (MQ=60)
                                AATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCCTCTTTTTTTATTCGCGGTACAAAAGCCCGTCATG                                                             >  SRR3722072.278972/1‑100 (MQ=60)
                                 ATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGGACAAGTGCCAGTCATGT                                                            >  SRR3722072.928714/1‑100 (MQ=60)
                                   GATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGT                                                          <  SRR3722072.782825/100‑1 (MQ=60)
                                     TGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTG                                                        <  SRR3722072.414678/100‑1 (MQ=60)
                                         cggagatgtgtataagagacaGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGT                                                    <  SRR3722072.562664/79‑1 (MQ=60)
                                                 TCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTCGCGGGACAAATGCCAGCCATGTTTTGGGGTTAATTGAT                                            >  SRR3722072.53291/1‑100 (MQ=60)
                                                 TCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCCGTCATGTGTTGGGGTTAATTGAT                                            >  SRR3722072.567113/1‑100 (MQ=60)
                                                               ATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTctgtctcttatacacatctccgagcccacgagacgtccgattttctcgtatgccgtct                              >  SRR3722072.626216/1‑42 (MQ=60)
                                                               ATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGGGTTAATTGATGGCGTACGCCTGAA                              >  SRR3722072.408837/1‑100 (MQ=60)
                                                                                 GTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGctgtctcttatacacatctccgagcccacgagacgtccgatta            >  SRR3722072.21011/1‑57 (MQ=60)
                                                                                 GTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG            >  SRR3722072.187658/1‑100 (MQ=60)
                                                                                      GGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTG       >  SRR3722072.35035/1‑100 (MQ=60)
                                                                                           GCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGAC  >  SRR3722072.816808/1‑100 (MQ=60)
                                                                                               |                                                                                               
AGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGAC  >  NZ_CP009273/4153721‑4153911

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