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 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 728,467 T→C I1143I (ATT→ATC rhsC → rhs element protein RhsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273728,4670TC100.0% 106.7 / NA 32I1143I (ATT→ATCrhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/17);  total (15/17)

AAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAAT  >  NZ_CP009273/728379‑728546
                                                                                        |                                                                               
aaaGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCa                                                                               <  2:401780/91‑1 (MQ=255)
    ccgccgCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTg                                                                           >  1:916711/1‑91 (MQ=37)
    ccgccgCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGGATACACGCCGGCGCGGAAAATCCACCTg                                                                           >  2:1055548/1‑91 (MQ=255)
     cgccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGt                                                                          <  1:199342/91‑1 (MQ=255)
           tggctggcATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACt                                                                    >  1:449577/1‑91 (MQ=37)
                      gTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCAcc                                                                             <  1:1088270/71‑1 (MQ=255)
                      gTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCAcc                                                                             >  2:1088270/1‑71 (MQ=255)
                               gACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCg                                                >  1:403223/1‑91 (MQ=255)
                                 cGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCCCTGCGACCATCGCGGCCTACCGCt                                              <  2:928623/91‑1 (MQ=255)
                                 cGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCt                                              <  2:649178/91‑1 (MQ=37)
                                         cAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTa                                      >  1:200392/1‑91 (MQ=18)
                                          aGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTAt                                     <  1:704713/91‑1 (MQ=18)
                                                aaaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCAc                               <  2:895634/91‑1 (MQ=21)
                                                  aGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACgg                             <  2:591272/91‑1 (MQ=25)
                                                            aCCCGGTATACACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGacaac                   <  2:952254/91‑1 (MQ=16)
                                                            aCCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGacaac                   <  2:840975/91‑1 (MQ=21)
                                                                 gTATACACGCCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGt              <  2:121942/91‑1 (MQ=255)
                                                                 gTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGt              >  1:108443/1‑91 (MQ=37)
                                                                      cacGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTAc         >  2:707153/1‑91 (MQ=255)
                                                                       acGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACg        >  1:1048948/1‑91 (MQ=255)
                                                                        cGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGc       <  2:200392/91‑1 (MQ=255)
                                                                        cGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGc       <  2:449577/91‑1 (MQ=255)
                                                                        cGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGc       <  1:734156/91‑1 (MQ=255)
                                                                        cGCCGGCGCGAAAAATCCACCTGTACCACGGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGc       <  2:875984/91‑1 (MQ=255)
                                                                         gCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCa      <  2:876826/91‑1 (MQ=255)
                                                                         gCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCa      >  2:1091552/1‑91 (MQ=255)
                                                                         gCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCa      >  2:556640/1‑91 (MQ=255)
                                                                          ccGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAg     >  2:95682/1‑91 (MQ=255)
                                                                           cGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGa    >  2:885671/1‑91 (MQ=255)
                                                                           cGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGa    >  1:54512/1‑91 (MQ=255)
                                                                             gcgcGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAAt  >  1:833419/1‑91 (MQ=255)
                                                                             gcgcGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAAt  >  1:805702/1‑91 (MQ=255)
                                                                             gcgcGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAAt  <  1:566436/91‑1 (MQ=255)
                                                                                        |                                                                               
AAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAAT  >  NZ_CP009273/728379‑728546

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 26 ≤ 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
GTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATG  >  NZ_CP009273/728373‑728556
                                                                                              |                                                                                         
GTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTG                                                                                     >  SRR3722125.937049/1‑101 (MQ=24)
    GGAAGGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACC                                                                                 <  SRR3722125.190341/101‑1 (MQ=23)
     GAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCA                                                                                <  SRR3722125.84361/101‑1 (MQ=24)
       AAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACT                                                                              >  SRR3722125.459062/1‑101 (MQ=24)
        taagagacagGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACCCCGGCGCGAAAAATCCACCTGTACCACTG                                                                             <  SRR3722125.1112766/91‑1 (MQ=24)
        AGCCGCCGCTGGCTGGCATCGTGCCGCCTGACGGTGGCGCAGAAGCAAAGCCAGATGGACCCGGTATACACCCCGGCGCGAAAAATCCACCTGTACCACTG                                                                             <  SRR3722125.529977/101‑1 (MQ=22)
           CGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGA                                                                          <  SRR3722125.203266/101‑1 (MQ=24)
                           CGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCG                                                          >  SRR3722125.411699/1‑101 (MQ=24)
                                     GACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTA                                                >  SRR3722125.204341/1‑101 (MQ=34)
                                                AGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAA                                     <  SRR3722125.719915/101‑1 (MQ=24)
                                                          gtgtataagagacagATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAG                           <  SRR3722125.605687/86‑1 (MQ=0)
                                                             GATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGT                        >  SRR3722125.110489/1‑101 (MQ=20)
                                                                   CCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACG                  >  SRR3722125.1072511/1‑101 (MQ=34)
                                                                    agGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGC                 <  SRR3722125.283003/99‑1 (MQ=34)
                                                                       GTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGA              >  SRR3722125.55515/1‑101 (MQ=39)
                                                                         ATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAAT            >  SRR3722125.823349/1‑101 (MQ=39)
                                                                         ATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAAT            >  SRR3722125.851732/1‑101 (MQ=39)
                                                                              CGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGAT       <  SRR3722125.749982/101‑1 (MQ=37)
                                                                                   GCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATG  <  SRR3722125.578558/101‑1 (MQ=37)
                                                                                              |                                                                                         
GTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATG  >  NZ_CP009273/728373‑728556

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