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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,840,928 T→C E97G (GAA→GGA)  topB ← DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,840,9280TC100.0% 126.3 / NA 37E97G (GAA→GGA) topBDNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/18);  total (19/18)

CGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCG  >  NZ_CP009273/1840849‑1841009
                                                                               |                                                                                 
cGGTGCCAGTTGCAGATAGTCCAGCCCTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGAGTGAACGATTCCGCTGGCTTc                                                                         <  1:1229346/90‑1 (MQ=255)
cGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTc                                                                         <  2:856754/90‑1 (MQ=255)
cGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTc                                                                         <  1:382612/90‑1 (MQ=255)
cGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTc                                                                         <  1:568775/90‑1 (MQ=255)
        gTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGa                                                                 <  2:724463/90‑1 (MQ=255)
                  gTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTtgat                                                       <  2:1016150/90‑1 (MQ=255)
                  gTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGGTtgat                                                       >  1:27275/1‑90 (MQ=255)
                          cTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTaa                                               >  1:203795/1‑90 (MQ=255)
                            tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                             >  2:480809/1‑90 (MQ=255)
                            tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                             >  1:910878/1‑90 (MQ=255)
                            tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                             >  1:70027/1‑90 (MQ=255)
                            tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                             >  1:576143/1‑90 (MQ=255)
                            tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                             >  1:380626/1‑90 (MQ=255)
                            tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                             >  1:30749/1‑90 (MQ=255)
                            tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATACCGCTGGCTTc                                                                         <  1:710081/62‑1 (MQ=255)
                            tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATACCGCTGGCTTc                                                                         >  2:710081/1‑62 (MQ=255)
                                ccaccaGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgttgtt                                         <  1:1026837/90‑1 (MQ=255)
                                 caccaGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGttt                                        <  1:491524/90‑1 (MQ=255)
                                 caccaGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGttt                                        <  2:110514/90‑1 (MQ=255)
                                 caccaGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGttt                                        <  2:684648/90‑1 (MQ=255)
                                      gCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCa                                   >  2:669178/1‑90 (MQ=255)
                                         aTTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACgg                                >  2:820698/1‑90 (MQ=255)
                                          ttGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGa                               <  1:480809/90‑1 (MQ=255)
                                               cTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTc                          <  2:1108496/90‑1 (MQ=255)
                                               cTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTc                          <  2:123887/90‑1 (MQ=255)
                                                ttCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCg                         >  2:966837/1‑90 (MQ=255)
                                                    cGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGc                     <  2:49558/90‑1 (MQ=255)
                                                       tCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGt                  >  2:86282/1‑90 (MQ=255)
                                                           ggTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAAtt              >  1:1071910/1‑90 (MQ=255)
                                                              ccccGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGcc           >  1:894916/1‑90 (MQ=255)
                                                                ccGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGa                               >  2:1106179/1‑68 (MQ=255)
                                                                ccGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGa                               <  1:1106179/68‑1 (MQ=255)
                                                                ccGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCAc         >  1:228431/1‑90 (MQ=255)
                                                                 cGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACt        <  1:478956/90‑1 (MQ=255)
                                                                 cGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACt        <  2:933819/90‑1 (MQ=255)
                                                                     gTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTc    >  2:247552/1‑90 (MQ=255)
                                                                       gAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCg  >  2:855095/1‑90 (MQ=255)
                                                                               |                                                                                 
CGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCG  >  NZ_CP009273/1840849‑1841009

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ 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
TGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCG  >  NZ_CP009273/1840829‑1841027
                                                                                                   |                                                                                                   
TGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTC                                                                                                     >  SRR3722089.299640/1‑100 (MQ=60)
                    gatgtgtataagagacagGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATACCGCTGGCTTCATGCAGGAAC                                                                                 <  SRR3722089.719327/82‑1 (MQ=60)
                    CGGTGCCAGTTGCAGATAGTCCAGCCCTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGAGTGAACGATTCCGCTGGCTTCATGCAGGAAC                                                                                 <  SRR3722089.1248029/100‑1 (MQ=60)
                    CGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAAC                                                                                 <  SRR3722089.575493/100‑1 (MQ=60)
                    CGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAAC                                                                                 <  SRR3722089.387274/100‑1 (MQ=60)
                            GTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGGTTGAT                                                                         >  SRR3722089.27605/1‑100 (MQ=60)
                             nnnnnnnnnnnnnnnnnnnnCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATG                                                                        <  SRR3722089.1064522/80‑1 (MQ=60)
                                    TAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAA                                                                 >  SRR3722089.206368/1‑100 (MQ=60)
                                      GTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGT                                                               >  SRR3722089.924129/1‑100 (MQ=60)
                                      GTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGT                                                               >  SRR3722089.70879/1‑100 (MQ=60)
                                      GTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGT                                                               >  SRR3722089.582980/1‑100 (MQ=60)
                                      GTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGT                                                               >  SRR3722089.385276/1‑100 (MQ=60)
                                      GTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGT                                                               >  SRR3722089.31120/1‑100 (MQ=60)
                                                    CCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGA                                                 <  SRR3722089.1042157/100‑1 (MQ=60)
                                                     CACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAG                                                <  SRR3722089.497346/100‑1 (MQ=60)
                                                              tataagagacaGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGC                                       <  SRR3722089.1122904/89‑1 (MQ=60)
                                                              TTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGC                                       <  SRR3722089.486520/100‑1 (MQ=60)
                                                                     TCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATT                                >  SRR3722089.1088006/1‑100 (MQ=60)
                                                                        CGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCC                             >  SRR3722089.907860/1‑100 (MQ=60)
                                                                          ATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCAC                           >  SRR3722089.231332/1‑100 (MQ=60)
                                                                                     CGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGAC                <  SRR3722089.484649/100‑1 (MQ=60)
                                                                                       GCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAA              >  SRR3722089.747342/1‑100 (MQ=60)
                                                                                           GAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGG          >  SRR3722089.652047/1‑100 (MQ=60)
                                                                                                   CCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCG  >  SRR3722089.416159/1‑100 (MQ=60)
                                                                                                   |                                                                                                   
TGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCG  >  NZ_CP009273/1840829‑1841027

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