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 4,004,326 T→G V256G (GTG→GGG)  yigL → sugar/pyridoxal phosphate phosphatase YigL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,004,3260TG100.0% 137.6 / NA 39V256G (GTG→GGG) yigLsugar/pyridoxal phosphate phosphatase YigL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (12/27);  total (12/27)

AGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATA  >  NZ_CP009273/4004240‑4004410
                                                                                      |                                                                                    
aGGCTGCATTATGGGCAGTGCGCATCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGcc                                                                                   <  2:110614/90‑1 (MQ=255)
     gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCAtt                                                                              >  1:852934/1‑90 (MQ=255)
       aTTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTAt                                                                            <  2:695011/90‑1 (MQ=255)
         tATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCt                                                                          <  2:769762/90‑1 (MQ=255)
         tATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCt                                                                          <  1:576708/90‑1 (MQ=255)
            gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCg                                                                       <  1:502395/90‑1 (MQ=255)
             ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGt                                                                      <  1:665475/90‑1 (MQ=255)
              gCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTa                                                                     <  1:646437/90‑1 (MQ=255)
               cAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTaa                                                                    >  2:559895/1‑90 (MQ=255)
                aGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTaaa                                                                   >  1:848748/1‑90 (MQ=255)
                 gTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAc                                                                  >  2:264603/1‑90 (MQ=255)
                 gTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAc                                                                  <  1:65923/90‑1 (MQ=255)
                 gTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAc                                                                  <  2:821260/90‑1 (MQ=255)
                  tGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAct                                                                 <  1:563569/90‑1 (MQ=255)
                  tGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAct                                                                 >  2:57835/1‑90 (MQ=255)
                   gcgcACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAActc                                                                >  2:214377/1‑90 (MQ=255)
                      cACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTAt                                                             >  2:310897/1‑90 (MQ=255)
                           gCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATc                                                        <  1:703711/90‑1 (MQ=255)
                           gCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATc                                                        >  1:455045/1‑90 (MQ=255)
                                tGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAAt                                                   <  1:316565/90‑1 (MQ=255)
                                 gAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  <  1:253844/90‑1 (MQ=255)
                                 gAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  <  2:40142/90‑1 (MQ=255)
                                 gAAAGACCTTCATCCCGAGCTGGAAGTAATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  <  2:326448/90‑1 (MQ=255)
                                     gACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTc                                              <  1:383958/90‑1 (MQ=255)
                                     gACCTTAATCCCGAGCTGGTTGGGATTGGCACTAATGCCGACGACGCGGGGCCGCTTAATCTTCGTAGAGTCTCTTTATCGTAATGGTTc                                              >  2:876883/1‑90 (MQ=255)
                                         ttCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  >  1:309899/1‑82 (MQ=255)
                                         ttCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  <  2:309899/82‑1 (MQ=255)
                                                 gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCa                                  <  1:739440/90‑1 (MQ=255)
                                                       gAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTc                            <  2:455045/90‑1 (MQ=255)
                                                                   aCTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACtt                <  2:920567/90‑1 (MQ=255)
                                                                    cTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTc               <  2:289949/90‑1 (MQ=255)
                                                                      aaTGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGc             >  2:790224/1‑90 (MQ=255)
                                                                        tGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTa           <  1:457551/90‑1 (MQ=255)
                                                                        tGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTa           <  1:358831/90‑1 (MQ=255)
                                                                         gCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTAc          <  2:154551/90‑1 (MQ=255)
                                                                           cgacgaCGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACaa        >  2:92741/1‑90 (MQ=255)
                                                                                acgCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGAt   <  2:401416/90‑1 (MQ=255)
                                                                                acgCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGAt   <  1:564556/90‑1 (MQ=255)
                                                                                 cgcgGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATAGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATa  <  1:518157/90‑1 (MQ=255)
                                                                                      |                                                                                    
AGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATA  >  NZ_CP009273/4004240‑4004410

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
GGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATC  >  NZ_CP009273/4004235‑4004420
                                                                                           |                                                                                              
GGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATT                                                                                        >  SRR3722072.867666/1‑100 (MQ=60)
  GAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTAT                                                                                      <  SRR3722072.841735/100‑1 (MQ=60)
   agaagacggcatacgagataatcggacgtctcgtgggctcggagatgtgtataagagacAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATC                                                                                     <  SRR3722072.510842/41‑1 (MQ=60)
           CATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAA                                                                             >  SRR3722072.863422/1‑100 (MQ=60)
              TATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTC                                                                          <  SRR3722072.587007/100‑1 (MQ=60)
                 GGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTAT                                                                       <  SRR3722072.511463/100‑1 (MQ=60)
                  GGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATT                                                                      <  SRR3722072.677373/100‑1 (MQ=60)
                   GCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTT                                                                     <  SRR3722072.657974/100‑1 (MQ=60)
                      GTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATC                                                                  <  SRR3722072.67140/100‑1 (MQ=60)
                      GTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATC                                                                  >  SRR3722072.463377/1‑100 (MQ=60)
                       TGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCG                                                                 <  SRR3722072.573672/100‑1 (MQ=60)
                                ggacgtctcgtgggctcggagatgtgtataagagacaGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTC                                                        <  SRR3722072.615996/63‑1 (MQ=60)
                                GCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTC                                                        <  SRR3722072.716269/100‑1 (MQ=60)
                                    CTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTA                                                    >  SRR3722072.315285/1‑100 (MQ=60)
                                     TGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTAT                                                   <  SRR3722072.322110/100‑1 (MQ=60)
                                      GAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATT                                                  <  SRR3722072.258149/100‑1 (MQ=60)
                                          GACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGT                                              <  SRR3722072.390821/100‑1 (MQ=60)
                                                      GAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACC                                  <  SRR3722072.752538/100‑1 (MQ=60)
                                                                             TGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATAC           <  SRR3722072.465929/100‑1 (MQ=60)
                                                                             TGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATAC           <  SRR3722072.365183/100‑1 (MQ=60)
                                                                               CCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACNC         <  SRR3722072.12001/100‑1 (MQ=60)
                                                                                     ACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAAT   <  SRR3722072.574674/100‑1 (MQ=60)
                                                                                      CGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATAGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATC  <  SRR3722072.527483/100‑1 (MQ=60)
                                                                                           |                                                                                              
GGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATC  >  NZ_CP009273/4004235‑4004420

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