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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,941,604 A→C G10G (GGT→GGG ruvC ← crossover junction endodeoxyribonuclease RuvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,941,6040AC100.0% 87.6 / NA 27G10G (GGT→GGGruvCcrossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (13/14);  total (13/14)

TCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC  >  NZ_CP009273/1941518‑1941686
                                                                                      |                                                                                  
tcCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAcccc                                                                                 >  1:313960/1‑90 (MQ=255)
tcCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAcccc                                                                                 <  1:535714/90‑1 (MQ=255)
     tttGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATc                                                                            <  2:250539/90‑1 (MQ=255)
      ttGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCa                                                                           >  1:622184/1‑90 (MQ=255)
         gTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATg                                                                        >  2:586291/1‑90 (MQ=255)
         gTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATg                                                                        >  1:366084/1‑90 (MQ=255)
            cGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg                                                                     <  1:131151/90‑1 (MQ=255)
            cGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg                                                                     <  2:322692/90‑1 (MQ=255)
            cGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg                                                                     <  1:639909/90‑1 (MQ=255)
             ggATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCga                                                                    >  1:160564/1‑90 (MQ=255)
                      cGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAg                                                           <  2:160107/90‑1 (MQ=255)
                              aGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCAcg                                                   >  2:336741/1‑90 (MQ=255)
                               ggTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCAcgc                                                  <  1:5401/90‑1 (MQ=255)
                                        gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGc                                                                       >  1:214141/1‑60 (MQ=255)
                                        gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGc                                                                       <  2:214141/60‑1 (MQ=255)
                                        gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  2:335969/1‑90 (MQ=255)
                                        gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  2:152612/1‑90 (MQ=255)
                                        gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  2:234730/1‑90 (MQ=255)
                                                cctacctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGtt                                 >  2:476330/1‑90 (MQ=255)
                                                     ctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCa                            >  1:399201/1‑90 (MQ=255)
                                                     ctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCa                            >  2:324141/1‑90 (MQ=255)
                                                      tGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAg                           <  2:244558/90‑1 (MQ=255)
                                                      tGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAg                           <  1:707417/90‑1 (MQ=255)
                                                        gCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGc                         <  2:149764/90‑1 (MQ=255)
                                                           gATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCtc                      <  2:499096/90‑1 (MQ=255)
                                                              gACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCtcatc                   <  2:101567/90‑1 (MQ=255)
                                                                cGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGGATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCag                 <  2:33045/90‑1 (MQ=255)
                                                                          ggTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTg       >  1:342322/1‑90 (MQ=255)
                                                                               cgcgcgACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAcc  <  1:476330/90‑1 (MQ=255)
                                                                                      |                                                                                  
TCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC  >  NZ_CP009273/1941518‑1941686

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
AGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACC  >  NZ_CP009273/1941505‑1941703
                                                                                                   |                                                                                                   
AGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAC                                                                                                     >  SRR3722112.64290/1‑100 (MQ=60)
   CGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCC                                                                                                  >  SRR3722112.317564/1‑100 (MQ=60)
        AATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATC                                                                                             <  SRR3722112.159062/100‑1 (MQ=60)
        AATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATC                                                                                             <  SRR3722112.226843/100‑1 (MQ=60)
         ATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCA                                                                                            >  SRR3722112.631789/1‑100 (MQ=60)
            ATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATG                                                                                         >  SRR3722112.370785/1‑100 (MQ=60)
             TCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGC                                                                                        <  SRR3722112.543696/100‑1 (MQ=60)
                ACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGA                                                                                     >  SRR3722112.162538/1‑100 (MQ=60)
                         CGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAG                                                                            <  SRR3722112.132779/100‑1 (MQ=60)
                         CGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAG                                                                            <  SRR3722112.649851/100‑1 (MQ=60)
                                           AGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATctgtctcttatacacatctg                                                          >  SRR3722112.216721/1‑80 (MQ=60)
                                            GGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTT                                                         <  SRR3722112.5464/100‑1 (MQ=60)
                                                        GCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCA                                             >  SRR3722112.404455/1‑100 (MQ=60)
                                                                   TGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAG                                  <  SRR3722112.718471/100‑1 (MQ=60)
                                                                             CGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTG                        >  SRR3722112.346562/1‑100 (MQ=60)
                                                                                       GGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATC              >  SRR3722112.613577/1‑100 (MQ=60)
                                                                                       GGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATC              >  SRR3722112.117633/1‑100 (MQ=60)
                                                                                            CGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGA         <  SRR3722112.483144/100‑1 (MQ=60)
                                                                                                   CCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACC  >  SRR3722112.572033/1‑100 (MQ=60)
                                                                                                   |                                                                                                   
AGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACC  >  NZ_CP009273/1941505‑1941703

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