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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,901,121 A→G D202D (GAT→GAC pstB ← phosphate ABC transporter ATP‑binding protein PstB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,901,1210AG100.0% 68.7 / NA 22D202D (GAT→GACpstBphosphate ABC transporter ATP‑binding protein PstB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/10);  total (12/10)

TCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCAC  >  NZ_CP009273/3901036‑3901207
                                                                                     |                                                                                      
tcaatTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCt                                                                                    <  2:332165/90‑1 (MQ=255)
  aatTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGc                                                                                  <  1:253860/90‑1 (MQ=255)
           aGGTACATAAACGCCGTGTGGTCGGGACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCg                                                                         >  2:228057/1‑90 (MQ=255)
               aCATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGa                                                                     <  2:270163/90‑1 (MQ=255)
                 aTAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAATTCGGTGATc                                                                   <  2:483558/90‑1 (MQ=255)
                      cGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGctc                                                              >  1:571519/1‑90 (MQ=255)
                       gCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGctct                                                             >  2:475319/1‑90 (MQ=255)
                       gCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGctct                                                             <  2:416967/90‑1 (MQ=255)
                                tCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGc                                                    >  2:298938/1‑90 (MQ=255)
                                      caacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTa                                              <  1:142065/90‑1 (MQ=255)
                                        acGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTcc                                                                                     <  2:484034/49‑1 (MQ=255)
                                        acGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTcc                                                                                     >  1:484034/1‑49 (MQ=255)
                                                    ctgcATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAg                                <  2:258013/90‑1 (MQ=255)
                                                           ttgtGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGaa                         >  2:496257/1‑90 (MQ=255)
                                                              tgGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAacac                      >  1:515458/1‑90 (MQ=255)
                                                                   gACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtc                 <  1:529276/90‑1 (MQ=255)
                                                                   gACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtc                 >  1:408565/1‑90 (MQ=255)
                                                                   gACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtc                 >  2:78703/1‑90 (MQ=255)
                                                                      gATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtcgtc              >  1:12217/1‑90 (MQ=255)
                                                                         caccacGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCg                                  >  1:584876/1‑67 (MQ=255)
                                                                         caccacGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCg                                  <  2:584876/67‑1 (MQ=255)
                                                                                  gtAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCAc  >  1:393160/1‑90 (MQ=255)
                                                                                     |                                                                                      
TCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCAC  >  NZ_CP009273/3901036‑3901207

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATG  >  NZ_CP009273/3901023‑3901220
                                                                                                  |                                                                                                   
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGT                                                                                                    <  SRR3722109.359051/100‑1 (MQ=60)
  GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                                  <  SRR3722109.130810/100‑1 (MQ=60)
  GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                                  <  SRR3722109.424421/100‑1 (MQ=60)
  GCCGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                                  <  SRR3722109.278712/100‑1 (MQ=60)
               AATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGG                                                                                     <  SRR3722109.256881/100‑1 (MQ=60)
                         GGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTC                                                                           >  SRR3722109.580187/1‑100 (MQ=60)
                                           GGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCtgtctcttatacacatctgacgctgccgac                                                         >  SRR3722109.491324/1‑70 (MQ=60)
                                                   CAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGT                                                 <  SRR3722109.143827/100‑1 (MQ=60)
                                                          attcgtcggcagcgtcagatgtgtataagagacaGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGC                                          <  SRR3722109.388311/66‑1 (MQ=60)
                                                                 CTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACAC                                   >  SRR3722109.523291/1‑100 (MQ=60)
                                                                      TGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTC                              >  SRR3722109.414517/1‑100 (MQ=60)
                                                                         TGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTC                           >  SRR3722109.12390/1‑100 (MQ=60)
                                                                            GGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACctgtctcttatac                        >  SRR3722109.593731/1‑87 (MQ=60)
                                                                                GACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGC                    <  SRR3722109.537331/100‑1 (MQ=60)
                                                                                     TCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCAC               >  SRR3722109.398798/1‑100 (MQ=60)
                                                                                          ACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCG          >  SRR3722109.227623/1‑100 (MQ=60)
                                                                                             GTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC       >  SRR3722109.70296/1‑100 (MQ=60)
                                                                                             GTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC       >  SRR3722109.79415/1‑100 (MQ=60)
                                                                                                  GTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATG  >  SRR3722109.533318/1‑100 (MQ=60)
                                                                                                  |                                                                                                   
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATG  >  NZ_CP009273/3901023‑3901220

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