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 F5 I100 R1 10 30.5 2136404 85.7% 1830898 81.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 859,587 G→A A164T (GCA→ACA)  fsa → fructose‑6‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273859,5870GA100.0% 90.0 / NA 24A164T (GCA→ACA) fsafructose‑6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/11);  total (13/11)

TATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGC  >  NZ_CP009273/859499‑859669
                                                                                        |                                                                                  
tATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACa                                                                                  >  2:801225/1‑91 (MQ=255)
                 gTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGt                                                                            >  1:40721/1‑80 (MQ=255)
                 gTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGt                                                                            <  2:40721/80‑1 (MQ=255)
                 gTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAAcccc                                                                 >  2:933141/1‑91 (MQ=255)
                      ggCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTc                                                            >  2:748798/1‑91 (MQ=255)
                      ggCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTc                                                            >  1:866523/1‑91 (MQ=255)
                         aTTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAgg                                                         >  2:95551/1‑91 (MQ=255)
                         aTTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAgg                                                         >  1:467382/1‑91 (MQ=255)
                                 tgtgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGAc                                                 <  1:962110/91‑1 (MQ=255)
                                  gtgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACt                                                >  1:532321/1‑91 (MQ=255)
                                    gACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGc                                              <  1:95551/91‑1 (MQ=255)
                                    gACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGc                                              <  1:140858/91‑1 (MQ=255)
                                        gACTTACACCAATTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTAc                                          >  1:1021695/1‑91 (MQ=255)
                                             acacCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCa                                     <  2:270896/91‑1 (MQ=255)
                                              cacCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAg                                    <  1:748798/91‑1 (MQ=255)
                                                   gTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtgt                               >  1:472499/1‑91 (MQ=255)
                                                        tGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATc                          >  2:810392/1‑91 (MQ=255)
                                                         gAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCa                         <  1:801225/91‑1 (MQ=255)
                                                           aaaTGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAAccc                                                                  >  2:481164/1‑48 (MQ=255)
                                                           aaaTGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAAccc                                                                  <  1:481164/48‑1 (MQ=255)
                                                             atgcatgcGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTa                     >  2:1036481/1‑91 (MQ=255)
                                                                      ccgcAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCAc            <  1:665156/91‑1 (MQ=255)
                                                                           ggCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGAt       <  1:706655/91‑1 (MQ=255)
                                                                                aaGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGc  <  2:440956/91‑1 (MQ=255)
                                                                                        |                                                                                  
TATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGC  >  NZ_CP009273/859499‑859669

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
TACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAAC  >  NZ_CP009273/859488‑859674
                                                                                                   |                                                                                       
ttcgtcggcagcgtcagatgtgtataagagacaGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGAC                                                                                        <  SRR3722123.459421/68‑1 (MQ=60)
                  GCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCC                                                                      >  SRR3722123.42484/1‑101 (MQ=60)
                       GGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTC                                                                 >  SRR3722123.907931/1‑101 (MQ=60)
                          CGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGG                                                              >  SRR3722123.489528/1‑101 (MQ=60)
                           tcgtcggcagcgtcagatgtgtataagagacagCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGC                                                             <  SRR3722123.503958/68‑1 (MQ=60)
                                   CATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACT                                                     >  SRR3722123.557650/1‑101 (MQ=60)
                                         GACTGTGACCGACTTACACCAATTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTAC                                               >  SRR3722123.1070925/1‑101 (MQ=60)
                                            TGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGG                                            <  SRR3722123.1008464/101‑1 (MQ=60)
                                               GACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAG                                         <  SRR3722123.147180/101‑1 (MQ=60)
                                               GACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAG                                         <  SRR3722123.99714/101‑1 (MQ=60)
                                                    ACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGT                                    >  SRR3722123.494887/1‑101 (MQ=60)
                                                         CACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATC                               <  SRR3722123.784437/101‑1 (MQ=60)
                                                                    GAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGC                    <  SRR3722123.839368/101‑1 (MQ=60)
                                                                                 CCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGC       <  SRR3722123.696720/101‑1 (MQ=60)
                                                                                      GGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAAC  <  SRR3722123.740213/101‑1 (MQ=60)
                                                                                                   |                                                                                       
TACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAAC  >  NZ_CP009273/859488‑859674

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