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 F1 I188 R2 3 28.7 1681193 95.1% 1598814 84.5

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% 122.3 / NA 31A164T (GCA→ACA) fsafructose‑6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/21);  total (10/21)

ATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACA  >  NZ_CP009273/859504‑859672
                                                                                   |                                                                                     
aTGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGa                                                                                 <  1:63776/90‑1 (MQ=255)
    tCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGttt                                                                             >  1:202759/1‑90 (MQ=255)
    tCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGttt                                                                             >  1:717871/1‑90 (MQ=255)
            gTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAAccc                                                                     >  2:291772/1‑90 (MQ=255)
                 ggCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGt                                                                >  2:654070/1‑90 (MQ=255)
                   cATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCa                                                              <  1:19822/90‑1 (MQ=255)
                        aGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGgcgc                                                         <  1:686770/90‑1 (MQ=255)
                            tgtgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGt                                                                               <  2:180799/64‑1 (MQ=255)
                            tgtgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGt                                                                               >  1:180799/1‑64 (MQ=255)
                               gACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTg                                                  >  1:350367/1‑90 (MQ=255)
                                         cacCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCa                                        <  2:643584/90‑1 (MQ=255)
                                              gTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtg                                   >  1:78639/1‑90 (MQ=255)
                                              gTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtg                                   >  2:50279/1‑90 (MQ=255)
                                              gTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtg                                   <  1:291772/90‑1 (MQ=255)
                                                tATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtgtg                                 <  1:825940/90‑1 (MQ=255)
                                                tATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtgtg                                 <  2:56885/90‑1 (MQ=255)
                                                  ttGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGaa                               >  2:170640/1‑90 (MQ=255)
                                                  ttGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGaa                               <  1:50279/90‑1 (MQ=255)
                                                        atgcatgcGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAAtt                         >  1:523154/1‑90 (MQ=255)
                                                            atgcGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTActc                     <  2:626713/90‑1 (MQ=255)
                                                            atgcGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTActc                     <  2:548862/90‑1 (MQ=255)
                                                              gcgcCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTg                   <  1:368315/90‑1 (MQ=255)
                                                                    cAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCTACTg             <  1:326771/90‑1 (MQ=255)
                                                                    cAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTg             <  1:750219/90‑1 (MQ=255)
                                                                    cAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTg             <  1:655363/90‑1 (MQ=255)
                                                                       gCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGAt          <  2:406344/90‑1 (MQ=255)
                                                                       gCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGAt          <  2:202578/90‑1 (MQ=255)
                                                                         gAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGAtgt        <  2:350367/90‑1 (MQ=255)
                                                                            aGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGc     <  2:717872/90‑1 (MQ=255)
                                                                              tGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGcac   <  1:654069/90‑1 (MQ=255)
                                                                               gCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGcaca  <  2:672274/90‑1 (MQ=255)
                                                                                   |                                                                                     
ATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACA  >  NZ_CP009273/859504‑859672

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGC  >  NZ_CP009273/859498‑859685
                                                                                         |                                                                                                  
GTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTT                                                                                          >  SRR3722064.206364/1‑100 (MQ=60)
GTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTT                                                                                          >  SRR3722064.732418/1‑100 (MQ=60)
      ATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAAC                                                                                    <  SRR3722064.64890/100‑1 (MQ=60)
                        GCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCtgtctcttatacaca                                                                  >  SRR3722064.184033/1‑85 (MQ=60)
                         CATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGAC                                                                 <  SRR3722064.20172/100‑1 (MQ=60)
                           TTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTG                                                               >  SRR3722064.356185/1‑100 (MQ=60)
                              AGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTT                                                            <  SRR3722064.700670/100‑1 (MQ=60)
                                          ACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATG                                                >  SRR3722064.80038/1‑100 (MQ=60)
                                                    GTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATT                                      <  SRR3722064.296739/100‑1 (MQ=60)
                                                    GTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATT                                      >  SRR3722064.533052/1‑100 (MQ=60)
                                                      TATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTAC                                    <  SRR3722064.842850/100‑1 (MQ=60)
                                                        TTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTC                                  <  SRR3722064.51176/100‑1 (MQ=60)
                                                                    GCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATG                      <  SRR3722064.374603/100‑1 (MQ=60)
                                                                          CAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCTACTGGATGTGGCAC                <  SRR3722064.332228/100‑1 (MQ=60)
                                                                          CAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCAC                <  SRR3722064.668554/100‑1 (MQ=60)
                                                                          CAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCAC                <  SRR3722064.765519/100‑1 (MQ=60)
                                                                                    TGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGAT      <  SRR3722064.667231/100‑1 (MQ=60)
                                                                                        GACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGC  <  SRR3722064.605333/100‑1 (MQ=60)
                                                                                         |                                                                                                  
GTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGC  >  NZ_CP009273/859498‑859685

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