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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,273,721 A→G E49G (GAA→GGA)  agaV → PTS N‑acetylgalactosamine transporter subunit IIB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,273,7210AG100.0% 98.0 / NA 28E49G (GAA→GGA) agaVPTS N‑acetylgalactosamine transporter subunit IIB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/16);  total (12/16)

TTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGAC  >  NZ_CP009273/3273634‑3273804
                                                                                       |                                                                                   
ttCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGgaa                                                                                   <  2:258268/90‑1 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                            <  2:606652/90‑1 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                            >  2:581124/1‑90 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                            >  1:499683/1‑90 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                            <  1:723944/90‑1 (MQ=255)
          cGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGc                                                                         <  1:455557/90‑1 (MQ=255)
                 gCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAggg                                                                  >  2:42778/1‑90 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                        <  2:499683/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                        <  1:415405/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                        <  1:776392/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                        <  1:805087/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                        <  2:837229/90‑1 (MQ=255)
                                gTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                        <  2:452281/85‑1 (MQ=255)
                                gTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                        >  1:452281/1‑85 (MQ=255)
                                     ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAc                                              >  1:843493/1‑90 (MQ=255)
                                       tAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGc                                            <  1:184982/90‑1 (MQ=255)
                                       tAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGc                                            <  2:200664/90‑1 (MQ=255)
                                                       tGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGa                            >  1:137065/1‑90 (MQ=255)
                                                       tGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGa                            >  2:830134/1‑90 (MQ=255)
                                                         ccGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGaca                          <  2:843493/90‑1 (MQ=255)
                                                                 ccGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCAt                  <  1:206248/90‑1 (MQ=255)
                                                                  cGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATc                                                               >  2:305505/1‑44 (MQ=255)
                                                                  cGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATc                                                               <  1:305505/44‑1 (MQ=255)
                                                                       aaacaaAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGcgcc            <  1:411328/89‑1 (MQ=255)
                                                                             aaCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGc                                            >  2:908838/1‑52 (MQ=255)
                                                                             aaCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGc                                            <  1:908838/52‑1 (MQ=255)
                                                                              aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc     >  2:896941/1‑90 (MQ=255)
                                                                                 tGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGAc  >  2:131638/1‑90 (MQ=255)
                                                                                 tGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGAc  >  1:382445/1‑90 (MQ=255)
                                                                                       |                                                                                   
TTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGAC  >  NZ_CP009273/3273634‑3273804

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ 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
AGGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGT  >  NZ_CP009273/3273625‑3273820
                                                                                                |                                                                                                   
AGGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAAT                                                                                                  <  SRR3722072.77712/100‑1 (MQ=60)
   TCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGT                                                                                               <  SRR3722072.840587/100‑1 (MQ=60)
     GGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAC                                                                                             <  SRR3722072.198093/100‑1 (MQ=60)
     GGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAC                                                                                             <  SRR3722072.801041/100‑1 (MQ=60)
     GGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAC                                                                                             <  SRR3722072.596932/100‑1 (MQ=60)
     GGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAC                                                                                             <  SRR3722072.46586/100‑1 (MQ=60)
      GCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACT                                                                                            >  SRR3722072.508711/1‑100 (MQ=60)
                GGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGG                                                                                  <  SRR3722072.736807/100‑1 (MQ=60)
                   CGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATC                                                                               <  SRR3722072.463902/100‑1 (MQ=60)
                             acgtctcgtgggctcggagatgtgtataagagacagGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTT                                                                     <  SRR3722072.310798/64‑1 (MQ=60)
                               GGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTC                                                                   >  SRR3722072.460560/1‑100 (MQ=60)
                                    cggcatacgagataatcggacgtctcgtgggctcggagatgtgtataagagacagGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              <  SRR3722072.595172/45‑1 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              <  SRR3722072.422923/100‑1 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              <  SRR3722072.790030/100‑1 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              <  SRR3722072.819131/100‑1 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                              >  SRR3722072.858095/1‑100 (MQ=60)
                                                tgggctcggagatgtgtataagagacagGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGT                                                  <  SRR3722072.924477/72‑1 (MQ=60)
                                                TAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGT                                                  <  SRR3722072.188196/100‑1 (MQ=60)
                                                      ACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGA                                            >  SRR3722072.139511/1‑100 (MQ=60)
                                                                          CCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCG                        <  SRR3722072.209780/100‑1 (MQ=60)
                                                                                CAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGAC                  >  SRR3722072.389285/1‑100 (MQ=60)
                                                                                AAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGAC                  <  SRR3722072.418775/100‑1 (MQ=60)
                                                                                                GAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGT  >  SRR3722072.675263/1‑100 (MQ=60)
                                                                                                |                                                                                                   
AGGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGT  >  NZ_CP009273/3273625‑3273820

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