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 F5 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 715,111 A→G I269T (ATC→ACC)  speF ← ornithine decarboxylase SpeF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273715,1110AG100.0% 147.8 / NA 40I269T (ATC→ACC) speFornithine decarboxylase SpeF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/25);  total (15/25)

ATAAATCGTACCGTCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGATCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCCGTTTC  >  NZ_CP009273/715023‑715190
                                                                                        |                                                                               
ataAATCGTACCGTCGTAGGTGCCTAACTGAATAACGGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTc                                                                               <  2:657969/91‑1 (MQ=255)
       gTACCGTCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCAc                                                                        >  2:1019894/1‑91 (MQ=255)
        tACCGTCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACg                                                                       >  1:1067807/1‑91 (MQ=255)
                  ggTGCCTAACTGAATCACAGGGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACtt                                                             <  1:945505/91‑1 (MQ=255)
                  ggTGCCTAACTGAATCACAGCGAGGGGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACtt                                                             <  2:195798/91‑1 (MQ=255)
                  ggTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACtt                                                             <  1:1025933/91‑1 (MQ=255)
                  ggTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACtt                                                             <  1:509676/91‑1 (MQ=255)
                  ggTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACtt                                                             <  1:275839/91‑1 (MQ=255)
                    tGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTttc                                                           <  1:67616/91‑1 (MQ=255)
                      ccTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCa                                                                     >  2:956524/1‑79 (MQ=255)
                      ccTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCa                                                                     <  1:956524/79‑1 (MQ=255)
                      ccTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTttctt                                                         >  1:15136/1‑91 (MQ=255)
                       cTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTttcttc                                                        >  1:2479/1‑91 (MQ=255)
                         aaCTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCaa                                                      >  1:463318/1‑91 (MQ=255)
                          aCTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCaaa                                                     <  2:977474/91‑1 (MQ=255)
                            tGAATCACAGCGAGGCGGAAAGGGCGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAc                                                   <  2:464568/91‑1 (MQ=255)
                            tGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTATGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAc                                                   <  1:697554/91‑1 (MQ=255)
                             gAATCACAGCGAGGCGGAAAGGGCGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAc                                                   <  2:475414/90‑1 (MQ=255)
                             gAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAc                                                   >  1:475414/1‑90 (MQ=255)
                               aTCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGt                                                >  2:946695/1‑91 (MQ=255)
                                 cacaGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTgc                                              <  2:947576/91‑1 (MQ=255)
                                 cacaGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTgc                                              <  2:27784/91‑1 (MQ=255)
                                   caGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACg                                                                       >  2:546458/1‑64 (MQ=255)
                                   caGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACg                                                                       <  1:546458/64‑1 (MQ=255)
                                             gAAAGGGCGAGCGTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCa                                  <  2:512926/91‑1 (MQ=255)
                                             gAAAGGAGGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCa                                  <  2:863662/91‑1 (MQ=255)
                                              aaaGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCAc                                 <  2:1067807/91‑1 (MQ=255)
                                                  gACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGAt                             >  2:1095964/1‑91 (MQ=255)
                                                    ggAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATaa                           <  2:723586/90‑1 (MQ=255)
                                                     gAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATaaa                          <  2:463318/91‑1 (MQ=255)
                                                      aGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAg                         <  1:802043/91‑1 (MQ=255)
                                                       gCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAc                                                   >  1:725687/1‑64 (MQ=255)
                                                       gCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAc                                                   <  2:725687/64‑1 (MQ=255)
                                                       gCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGc                        <  2:215659/91‑1 (MQ=255)
                                                           ctTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATa                    >  1:108767/1‑91 (MQ=255)
                                                                    gCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACg           >  2:453512/1‑91 (MQ=255)
                                                                     cTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGt          >  1:1099920/1‑91 (MQ=255)
                                                                     cTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGt          >  1:321752/1‑91 (MQ=255)
                                                                             cGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCCGTTTc  <  1:1019894/91‑1 (MQ=255)
                                                                             cGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCCGTTTc  <  1:569645/91‑1 (MQ=255)
                                                                                        |                                                                               
ATAAATCGTACCGTCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGATCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCCGTTTC  >  NZ_CP009273/715023‑715190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ 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
TTATAAATCGTACCGTCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGATCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCCGTTTCCAGATAAACCGGTGTTGCACC  >  NZ_CP009273/715021‑715211
                                                                                          |                                                                                                    
TTATAAATCGTACCGTCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACG                                                                                            >  SRR3722125.1091879/1‑101 (MQ=60)
          atgtgtataagagacagAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTT                                                                                  <  SRR3722125.558140/84‑1 (MQ=60)
            agGTCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTC                                                                                <  SRR3722125.977821/99‑1 (MQ=60)
              GTCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTT                                                                              >  SRR3722125.15418/1‑101 (MQ=60)
               TCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTC                                                                             >  SRR3722125.2541/1‑101 (MQ=60)
                 GTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAA                                                                           >  SRR3722125.473088/1‑101 (MQ=60)
                    GGTGCCTAACTGAATCACAGGGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAC                                                                        <  SRR3722125.966528/101‑1 (MQ=60)
                    GGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAG                                                                        <  SRR3722125.1048928/101‑1 (MQ=60)
                    GGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAC                                                                        <  SRR3722125.281413/101‑1 (MQ=60)
                    GGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAAC                                                                        <  SRR3722125.520486/101‑1 (MQ=60)
                     GTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGCGCGACTTCCGCGGTCAGCTCACGCAGGCAACTTTCTTCAAAACA                                                                       >  SRR3722125.485439/1‑101 (MQ=60)
                      TGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAG                                                                      <  SRR3722125.68855/101‑1 (MQ=60)
                              TGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTATGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATC                                                              <  SRR3722125.712588/101‑1 (MQ=60)
                                               GAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCctgtctcttatacacat                                             >  SRR3722125.741328/1‑84 (MQ=60)
                                                   GGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATA                                         >  SRR3722125.110816/1‑101 (MQ=60)
                                                        AGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGT                                    <  SRR3722125.819572/101‑1 (MQ=60)
                                                             CTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGT                               >  SRR3722125.1124695/1‑101 (MQ=60)
                                                             CTTTTGCCCGCTGCGGTGCGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGT                               >  SRR3722125.328361/1‑101 (MQ=60)
                                                                               CGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCCGTTTCCAGATAAACC             <  SRR3722125.1042743/101‑1 (MQ=60)
                                                                               CGACTTCCGCGGTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCCGTTTCCAGATAAACC             <  SRR3722125.581824/101‑1 (MQ=60)
                                                                                          GTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCtgtctcttatacacatctccgagccca  >  SRR3722125.442370/1‑74 (MQ=60)
                                                                                          GTCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCCGTTTCCAGATAAACCGGTGTTGCACC  <  SRR3722125.370645/101‑1 (MQ=60)
                                                                                          |                                                                                                    
TTATAAATCGTACCGTCGTAGGTGCCTAACTGAATCACAGCGAGGCGGAAAGGACGAGCCTCTTTTGCCCGCTGCGGTGCGACTTCCGCGATCAGCTCACGCAGGTAACTTTCTTCAAAACAGTGCGCATCAATGCCACCGATAAAGCCATACGGGTTACGTGCCGTTTCCAGATAAACCGGTGTTGCACC  >  NZ_CP009273/715021‑715211

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