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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

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% 85.7 / NA 25E97G (GAA→GGA) topBDNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/16);  total (9/16)

CTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCC  >  NZ_CP009273/1840843‑1841008
                                                                                     |                                                                                
cTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCt                                                                              <  2:232888/90‑1 (MQ=255)
cTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCt                                                                              <  1:756836/90‑1 (MQ=255)
    tCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCt                                                                          <  1:680717/90‑1 (MQ=255)
                tGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAAc                                                              <  1:150288/90‑1 (MQ=255)
                tGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAAc                                                              <  2:532320/90‑1 (MQ=255)
                tGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAAc                                                              <  2:420431/90‑1 (MQ=255)
                tGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAAc                                                              <  1:710965/90‑1 (MQ=255)
                       aatCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTtga                                                       <  1:469709/88‑1 (MQ=255)
                        gTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTtgat                                                      <  2:2217/90‑1 (MQ=255)
                             gCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGt                                                 <  1:173375/90‑1 (MQ=255)
                             gCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGt                                                 <  1:225183/90‑1 (MQ=255)
                                  tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                            >  2:710664/1‑90 (MQ=255)
                                  tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                            >  2:706182/1‑90 (MQ=255)
                                  tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                            >  2:613571/1‑90 (MQ=255)
                                  tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                            >  1:759858/1‑90 (MQ=255)
                                  tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                            >  1:684750/1‑90 (MQ=255)
                                  tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                            >  1:510932/1‑90 (MQ=255)
                                  tCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                            >  1:142625/1‑90 (MQ=255)
                                          caGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGt                                    >  1:694751/1‑90 (MQ=255)
                                                     cTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTc                         <  2:165662/90‑1 (MQ=255)
                                                          cGATCAGGGTCCCCGGCGTGAAGGATTCCACTGGCTACAGGCAGGAACCGTTTGATGACGTTAAGTCGTTTGGTCACGCAGGGTCGGGGc                    <  1:115687/90‑1 (MQ=255)
                                                                gggTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAAt              >  1:599087/1‑90 (MQ=255)
                                                                       cGGCGTGGACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGCCACGGAGGGTCGGGGCTGTAATTGCCACt       <  1:363351/90‑1 (MQ=255)
                                                                       cGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACt       <  2:599087/90‑1 (MQ=255)
                                                                            tGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTcc  <  2:182742/90‑1 (MQ=255)
                                                                                     |                                                                                
CTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCC  >  NZ_CP009273/1840843‑1841008

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ 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
CTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGAC  >  NZ_CP009273/1840837‑1841013
                                                                                           |                                                                                     
CTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCT                                                                               <  SRR3722073.779070/100‑1 (MQ=60)
      CTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGC                                                                         <  SRR3722073.772138/100‑1 (MQ=60)
          TCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGA                                                                     <  SRR3722073.694618/100‑1 (MQ=60)
                      TGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGA                                                         <  SRR3722073.153365/100‑1 (MQ=60)
                      TGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGA                                                         <  SRR3722073.725332/100‑1 (MQ=60)
                             AATCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAG                                                  <  SRR3722073.479716/100‑1 (MQ=60)
                              GTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGT                                                 >  SRR3722073.145582/1‑100 (MQ=60)
                              GTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGT                                                 >  SRR3722073.521700/1‑100 (MQ=60)
                              GTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGT                                                 >  SRR3722073.698709/1‑100 (MQ=60)
                              GTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGT                                                 >  SRR3722073.775208/1‑100 (MQ=60)
                                   GCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTT                                            <  SRR3722073.176761/100‑1 (MQ=60)
                                   GCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTT                                            <  SRR3722073.229878/100‑1 (MQ=60)
                                      CTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGT                                         >  SRR3722073.708873/1‑100 (MQ=60)
                                                            TTCACGATCCGGGTCCCCGGCGTGAACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAAT                   >  SRR3722073.611575/1‑100 (MQ=60)
                                                                CGATCAGGGTCCCCGGCGTGAAGGATTCCACTGGCTACAGGCAGGAACCGTTTGATGACGTTAAGTCGTTTGGTCACGCAGGGTCGGGGCTGTAATTGCG               <  SRR3722073.118165/100‑1 (MQ=60)
                                                                             CGGCGTGGACGATTCCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGCCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGAC  <  SRR3722073.371061/100‑1 (MQ=60)
                                                                                           |                                                                                     
CTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGAC  >  NZ_CP009273/1840837‑1841013

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