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 F4 I88 R1 217 30.5 1742422 95.9% 1670982 87.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 736,034 A→T I358F (ATC→TTC)  phrB → deoxyribodipyrimidine photo‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273736,0340AT100.0% 116.2 / NA 30I358F (ATC→TTC) phrBdeoxyribodipyrimidine photo‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/12);  total (18/12)

CCGCTATGCGTCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGATCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAACGGTGGCTGGCAGTGGG  >  NZ_CP009273/735948‑736118
                                                                                      |                                                                                    
ccGCTATGCGTCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGa                                                                                  >  1:695976/1‑91 (MQ=255)
   cTATGCGTCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTg                                                                               >  2:61945/1‑91 (MQ=255)
   cTATGCGGCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTg                                                                               >  1:809304/1‑91 (MQ=255)
              cTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGgcg                                                                    <  1:369128/91‑1 (MQ=255)
               ttAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGgcga                                                                   <  2:10361/91‑1 (MQ=255)
                 aaCAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGgcgagc                                                                 <  2:567256/91‑1 (MQ=255)
                 aaCAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGgcgagc                                                                 <  2:17548/91‑1 (MQ=255)
                     gCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGata                                                             >  1:487971/1‑91 (MQ=255)
                     gCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGata                                                             >  2:210551/1‑91 (MQ=255)
                               gATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGgcg                                                                    >  2:825821/1‑74 (MQ=255)
                               gATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGgcg                                                                    <  1:825821/74‑1 (MQ=255)
                                aTGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGc                                                  <  2:400377/91‑1 (MQ=255)
                                 tGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCa                                                 <  2:442868/91‑1 (MQ=255)
                                   cATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGc                                               <  1:726441/91‑1 (MQ=255)
                                   cATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGc                                               <  2:176061/91‑1 (MQ=255)
                                    aTAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCt                                              >  2:319148/1‑91 (MQ=255)
                                    aTAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCt                                              >  2:864845/1‑91 (MQ=255)
                                          ggCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCtgattga                                        >  2:118170/1‑91 (MQ=255)
                                            cTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATg                                      >  2:309735/1‑91 (MQ=255)
                                            cTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATg                                      >  1:535671/1‑91 (MQ=255)
                                             tACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTACTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATgg                                     >  1:329897/1‑91 (MQ=255)
                                              aCGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGt                                    >  1:447554/1‑91 (MQ=255)
                                              aCGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGt                                    >  1:349493/1‑91 (MQ=255)
                                               cGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTg                                   >  1:482756/1‑91 (MQ=255)
                                                      ttACAGTCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGc                            <  1:748413/91‑1 (MQ=255)
                                                        aCAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAg                          <  1:745229/91‑1 (MQ=255)
                                                                         gAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAACGGtggctgg         >  1:823999/1‑91 (MQ=255)
                                                                                ttATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGt                                    <  1:599640/57‑1 (MQ=255)
                                                                                ttATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGt                                    >  2:599640/1‑57 (MQ=255)
                                                                                ttATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAACGGTGGCTGGCAGTggg  >  1:501415/1‑91 (MQ=255)
                                                                                      |                                                                                    
CCGCTATGCGTCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGATCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAACGGTGGCTGGCAGTGGG  >  NZ_CP009273/735948‑736118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ 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
CCCGATTGTTGATGCCGCTATGCGTCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGATCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAACGGTGGCTGGCAGTGGG  >  NZ_CP009273/735934‑736118
                                                                                                    |                                                                                    
CCCGATTGTTGATGCCGCTATGCGTCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGT                                                                                      >  SRR3722118.455767/1‑101 (MQ=60)
CCCGATTGTTGATGCCGCTATGCGTCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGT                                                                                      >  SRR3722118.204347/1‑101 (MQ=60)
    ATTGTTGATGCCGCTATGCGTCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGA                                                                                  >  SRR3722118.711148/1‑101 (MQ=60)
       GTTGATGCCGCTATGCGGCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTG                                                                               >  SRR3722118.827095/1‑101 (MQ=60)
                         CAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATA                                                             >  SRR3722118.498294/1‑101 (MQ=60)
                            gagacagGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTT                                                          <  SRR3722118.844007/94‑1 (MQ=60)
                            CTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTT                                                          <  SRR3722118.376825/101‑1 (MQ=60)
                                          acgagataggagtccgtctcgtgggctcggagatgtgtataagagacagGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGA                                            <  SRR3722118.847507/52‑1 (MQ=60)
                                                GCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATG                                      >  SRR3722118.547036/1‑101 (MQ=60)
                                                 CATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGG                                     <  SRR3722118.742400/101‑1 (MQ=60)
                                                 CATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTACTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGG                                     >  SRR3722118.336674/1‑101 (MQ=60)
                                                  ATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGT                                    >  SRR3722118.456994/1‑101 (MQ=60)
                                                  ATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGT                                    >  SRR3722118.356741/1‑101 (MQ=60)
                                                   TAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTG                                   >  SRR3722118.492972/1‑101 (MQ=60)
                                                            ctcggagatgtgtataagagacagGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAG                          <  SRR3722118.612429/77‑1 (MQ=60)
                                                                    TTACAGTCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAAC                  <  SRR3722118.764843/101‑1 (MQ=60)
                                                                      ACAGCCAGTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAACGG                <  SRR3722118.761605/101‑1 (MQ=60)
                                                                             GTTTTCTGGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAACGGTGGCTGG         >  SRR3722118.842152/1‑101 (MQ=60)
                                                                                    GGTGAAAGATTTATTGTTCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAACGGTGGCTGGCAGTGGG  >  SRR3722118.512052/1‑101 (MQ=60)
                                                                                                    |                                                                                    
CCCGATTGTTGATGCCGCTATGCGTCAGCTTAACAGCACTGGCTGGATGCATAACAGGCTACGGATGATTACAGCCAGTTTTCTGGTGAAAGATTTATTGATCGACTGGCGCGAAGGCGAGCGATATTTCATGTCGCAGCTGATTGATGGTGATTTGGCAGCCAATAACGGTGGCTGGCAGTGGG  >  NZ_CP009273/735934‑736118

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