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 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% 92.0 / NA 27E49G (GAA→GGA) agaVPTS N‑acetylgalactosamine transporter subunit IIB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/14);  total (13/14)

TCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATC  >  NZ_CP009273/3273635‑3273801
                                                                                      |                                                                                
tCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAAt                                                                               <  1:488112/90‑1 (MQ=255)
tCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAAt                                                                               <  2:1216529/90‑1 (MQ=255)
 cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                              >  2:1012499/1‑90 (MQ=255)
 cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                              >  1:925954/1‑90 (MQ=255)
      ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                         <  2:872101/90‑1 (MQ=255)
                gCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAggg                                                               >  1:311381/1‑90 (MQ=255)
                 cGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGa                                                              <  2:204407/90‑1 (MQ=255)
                          aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                     <  1:803426/90‑1 (MQ=255)
                          aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                     <  2:1214698/90‑1 (MQ=255)
                          aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                     <  2:445739/90‑1 (MQ=255)
                          aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                     <  2:343170/90‑1 (MQ=255)
                                    ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                     <  1:643756/80‑1 (MQ=255)
                                    ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                     >  2:643756/1‑80 (MQ=255)
                                    ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAc                                           >  2:1116264/1‑90 (MQ=255)
                                     gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                          <  2:925954/90‑1 (MQ=255)
                                     gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                          <  2:876088/90‑1 (MQ=255)
                                     gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                          <  1:150617/90‑1 (MQ=255)
                                         ccAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCGGc                                      >  1:532925/1‑90 (MQ=255)
                                                  aGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTa                             <  2:40180/90‑1 (MQ=255)
                                                    gTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATc                           <  2:1111094/90‑1 (MQ=255)
                                                                       aacaaAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGcgccg        >  2:290680/1‑90 (MQ=255)
                                                                       aacaaAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGcgccg        >  2:180949/1‑90 (MQ=255)
                                                                       aacaaAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGcgccg        >  2:1105155/1‑90 (MQ=255)
                                                                           aaaCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGa    >  2:263527/1‑90 (MQ=255)
                                                                             aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc  >  2:36503/1‑90 (MQ=255)
                                                                             aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc  >  2:502010/1‑90 (MQ=255)
                                                                             aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc  >  2:1167426/1‑90 (MQ=255)
                                                                                      |                                                                                
TCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATC  >  NZ_CP009273/3273635‑3273801

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
AGGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTG  >  NZ_CP009273/3273625‑3273815
                                                                                                |                                                                                              
AGGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAAT                                                                                             <  SRR3722089.189763/100‑1 (MQ=60)
 GGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATG                                                                                            >  SRR3722089.939462/1‑100 (MQ=60)
     GGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGTTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAC                                                                                        <  SRR3722089.844478/100‑1 (MQ=60)
     GGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAC                                                                                        <  SRR3722089.668439/100‑1 (MQ=60)
          TCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCA                                                                                   <  SRR3722089.493906/100‑1 (MQ=60)
                GGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGG                                                                             >  SRR3722089.315326/1‑100 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                         <  SRR3722089.651979/100‑1 (MQ=60)
                                    ATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAC                                                         <  SRR3722089.814581/100‑1 (MQ=60)
                                         GTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCGGC                                                    >  SRR3722089.539136/1‑100 (MQ=60)
                                               GTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAG                                              <  SRR3722089.152487/100‑1 (MQ=60)
                                                                                        CCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATC     >  SRR3722089.957391/1‑100 (MQ=60)
                                                                                           GATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTG  >  SRR3722089.1192676/1‑100 (MQ=60)
                                                                                                |                                                                                              
AGGTCGGCGTTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTG  >  NZ_CP009273/3273625‑3273815

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