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 F1 I1 R1 66 46.7 5280542 98.1% 5180211 47.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_017635 2,322,732 T→C intergenic (‑40/+60) ECW_RS11650 ← / ← ECW_RS11655 head completion/stabilization protein/terminase endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176352,322,7320TC100.0% 94.4 / NA 36intergenic (‑40/+60)ECW_RS11650/ECW_RS11655head completion/stabilization protein/terminase endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/15);  total (21/15)

CATCATGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGTGTCAGGTGATTCACCGTCACTGACCGGCGTGCCGCCCTGGCGC  >  NC_017635/2322687‑2322775
                                             |                                           
catcatGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAgcgt                                           >  1:1099140/1‑48 (MQ=255)
catcatGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAgcgt                                           >  2:1696692/1‑48 (MQ=255)
  tcatGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCa                                         <  1:952202/48‑1 (MQ=255)
   catGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAg                                        >  1:1179824/1‑48 (MQ=255)
       aCTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGa                                    <  1:1839959/48‑1 (MQ=255)
        cTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGAt                                   >  1:2086976/1‑48 (MQ=255)
        cTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGAt                                   >  2:492344/1‑48 (MQ=255)
        cTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGAt                                   >  2:409839/1‑48 (MQ=255)
         tACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGAtt                                  >  1:18702/1‑48 (MQ=255)
         tACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGAtt                                  >  2:830022/1‑48 (MQ=255)
           ccTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTc                                 <  1:1567258/47‑1 (MQ=255)
           ccTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCa                                <  2:2265541/48‑1 (MQ=255)
            ctctGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCAc                               <  2:839517/48‑1 (MQ=255)
            ctctGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCAc                               >  2:1687792/1‑48 (MQ=255)
            ctctGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCAc                               >  1:800604/1‑48 (MQ=255)
            ctctGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCAc                               <  1:473538/48‑1 (MQ=255)
               tGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCAcc                              <  2:1514182/46‑1 (MQ=255)
               tGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCAcc                              <  2:1874585/46‑1 (MQ=255)
                gAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCAccct                            >  2:2508710/1‑45 (MQ=11)
                 aaaaGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTc                           <  1:910531/47‑1 (MQ=11)
                 aaaaGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCa                          <  1:2279796/48‑1 (MQ=11)
                  aaaGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCAc                         >  2:2049858/1‑48 (MQ=11)
                     gTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCACTGa                      <  1:701589/48‑1 (MQ=11)
                     gTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCACTGa                      >  2:230051/1‑48 (MQ=11)
                     gTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCACTGa                      <  2:210158/48‑1 (MQ=11)
                     gTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCACTGa                      <  1:456607/48‑1 (MQ=11)
                              tGGACGCCGGTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTgc             >  2:1978709/1‑48 (MQ=11)
                                 aCGCCGGTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTgccgc          >  1:2488885/1‑48 (MQ=11)
                                  cGCCGGTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTgccgcc         <  2:153596/48‑1 (MQ=11)
                                   gCCGGTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGccc        <  2:2614180/48‑1 (MQ=11)
                                        tACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGgcg   >  1:1294654/1‑48 (MQ=11)
                                        tACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGgcg   >  1:2553648/1‑48 (MQ=11)
                                        tACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGgcg   >  1:1972418/1‑48 (MQ=11)
                                         aCAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGgcgc  >  2:2190977/1‑48 (MQ=11)
                                         aCAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGgcgc  >  1:1383250/1‑48 (MQ=11)
                                         aCAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGgcgc  >  2:286283/1‑48 (MQ=11)
                                         aCAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGgcgc  >  1:2534793/1‑48 (MQ=11)
                                         aCAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGgcgc  >  1:509743/1‑48 (MQ=11)
                                             |                                           
CATCATGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGTGTCAGGTGATTCACCGTCACTGACCGGCGTGCCGCCCTGGCGC  >  NC_017635/2322687‑2322775

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTCATCATGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGTGTCAGGTGATTCACCGTCACTGACCGGCGTGCCGCCCTGGCGCGGGGCG  >  NC_017635/2322685‑2322781
                                               |                                                 
GTCATCATGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTC                                                >  M01420:84:000000000‑ABVEH:1:1110:3836:13179/1‑51 (MQ=10)
   ATCATGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGG                                             >  M01420:84:000000000‑ABVEH:1:1111:20246:5201/1‑51 (MQ=13)
   ATCATGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGG                                             <  M01420:84:000000000‑ABVEH:1:1109:19209:7380/51‑1 (MQ=13)
        GACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATT                                        >  M01420:84:000000000‑ABVEH:1:2104:9544:22457/1‑51 (MQ=13)
        GACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATT                                        <  M01420:84:000000000‑ABVEH:1:2102:20051:18830/51‑1 (MQ=13)
        GACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATT                                        >  M01420:84:000000000‑ABVEH:1:2113:9832:16084/1‑51 (MQ=13)
         ACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTC                                       >  M01420:84:000000000‑ABVEH:1:1101:17225:7432/1‑51 (MQ=13)
            ACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACC                                    >  M01420:84:000000000‑ABVEH:1:1107:7421:24059/1‑51 (MQ=13)
            ACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCAC                                     <  M01420:84:000000000‑ABVEH:1:1114:10083:11367/50‑1 (MQ=13)
             CCTCTGAGAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCC                                   <  M01420:84:000000000‑ABVEH:1:2110:8762:16375/51‑1 (MQ=60)
             CCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCC                                   <  M01420:84:000000000‑ABVEH:1:2110:9001:6612/51‑1 (MQ=13)
             CCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCC                                   >  M01420:84:000000000‑ABVEH:1:2110:25164:15614/1‑51 (MQ=13)
             CCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCC                                   <  M01420:84:000000000‑ABVEH:1:1105:17107:4914/51‑1 (MQ=13)
                CTGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTC                                 >  M01420:84:000000000‑ABVEH:1:2110:12896:17451/1‑50 (MQ=13)
                 TGAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCAC                               >  M01420:84:000000000‑ABVEH:1:2111:9059:26283/1‑51 (MQ=12)
                  GAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCACT                              <  M01420:84:000000000‑ABVEH:1:2106:2593:14909/51‑1 (MQ=12)
                  GAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCAC                               <  M01420:84:000000000‑ABVEH:1:1108:6694:22411/50‑1 (MQ=12)
                  GAAAAGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCA                                <  M01420:84:000000000‑ABVEH:1:2114:25381:19588/49‑1 (MQ=23)
                      AGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCACTGACC                          <  M01420:84:000000000‑ABVEH:1:1107:20377:3085/51‑1 (MQ=60)
                      AGTCGGGCGGTGGACGCCGGTACAGCGTCAGGTGATTCACCCTCACTGACC                          <  M01420:84:000000000‑ABVEH:1:1104:24273:25854/51‑1 (MQ=60)
                                 GGACGCCGGTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCC               >  M01420:84:000000000‑ABVEH:1:2108:15796:9841/1‑51 (MQ=12)
                                 GGACGCCGGTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCC               >  M01420:84:000000000‑ABVEH:1:2113:4714:9095/1‑51 (MQ=12)
                                        GGTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGGCGC        >  M01420:84:000000000‑ABVEH:1:2108:22479:22189/1‑51 (MQ=6)
                                        GGTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGGCGC        >  M01420:84:000000000‑ABVEH:1:2103:17881:22443/1‑51 (MQ=6)
                                        GGTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGGCGC        >  M01420:84:000000000‑ABVEH:1:1112:8518:3920/1‑51 (MQ=6)
                                         GTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGGCGCG       >  M01420:84:000000000‑ABVEH:1:2108:13723:18491/1‑51 (MQ=3)
                                         GTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGGCGCG       >  M01420:84:000000000‑ABVEH:1:1112:26402:21964/1‑51 (MQ=3)
                                         GTACAGCGTCAGGTGATTCACCCTCACTGACCGGCGTGCCGCCCTGGCGCG       >  M01420:84:000000000‑ABVEH:1:1105:26464:12798/1‑51 (MQ=3)
                                              GTGTCAGGTGATTCACCGTCACTGACCGGCGTGCCGCCCTGGCGCGGGGCG  >  M01420:84:000000000‑ABVEH:1:1109:18847:6391/1‑51 (MQ=0)
                                               |                                                 
GTCATCATGACTACCTCTGAAAAGTCGGGCGGTGGACGCCGGTACAGTGTCAGGTGATTCACCGTCACTGACCGGCGTGCCGCCCTGGCGCGGGGCG  >  NC_017635/2322685‑2322781

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