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
A0 F0 I1 R1 61 39.0 4165548 97.4% 4057243 49.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_017635 2,331,461 (T)8→7 intergenic (‑905/+261) ECW_RS11680 ← / ← ECW_RS11685 AAA family ATPase/replication endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176352,331,4540T.100.0% 182.0 / NA 42intergenic (‑898/+268)ECW_RS11680/ECW_RS11685AAA family ATPase/replication endonuclease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/26);  total (16/26)

CATTATGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGAC  >  NC_017635/2331406‑2331497
                                                |                                           
cATTATGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtt                                            >  1:1493376/1‑50 (MQ=255)
 aTTATGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGttt                                           >  1:485816/1‑50 (MQ=255)
 aTTATGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGttt                                           >  1:364425/1‑50 (MQ=255)
 aTTATGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGttt                                           >  2:1302278/1‑50 (MQ=255)
 aTTATGCTTTAATATGCTTCAAACAAAAAAACCACCTGACAGGGTGGttt                                           <  1:1235061/50‑1 (MQ=255)
  ttATGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtttt                                          <  2:1938401/50‑1 (MQ=255)
  ttATGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtttt                                          <  2:197430/50‑1 (MQ=255)
   tATGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGttttt                                         >  1:1456554/1‑50 (MQ=255)
    aTGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtttttt                                        >  2:1614998/1‑50 (MQ=255)
     tGCTTTAATATTCTTCAAACAAAAAAACCACCTGCCAGGGTGGttttttt                                       <  2:485668/50‑1 (MQ=255)
     tGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGttttttt                                       <  1:560346/50‑1 (MQ=255)
     tGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGttttttt                                       <  1:1897434/50‑1 (MQ=255)
     tGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtttttt                                        <  1:140718/49‑1 (MQ=255)
      gCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtttttttg                                      >  1:2054785/1‑49 (MQ=255)
      gCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtttttttg                                      >  2:592947/1‑49 (MQ=255)
       cTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtttttttgc                                     >  2:126814/1‑48 (MQ=255)
       cTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtttttttgc                                     >  1:158885/1‑48 (MQ=255)
       cTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGtttttttgc                                     >  1:1277572/1‑48 (MQ=255)
           aaTATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGcc                                 <  1:1811290/49‑1 (MQ=38)
           aaTATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGcc                                 <  2:906521/49‑1 (MQ=38)
            atatGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCca                               >  2:940876/1‑50 (MQ=39)
             tatGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCcat                              >  2:1054335/1‑50 (MQ=39)
                gCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCcatcat                           >  1:1637546/1‑50 (MQ=39)
                 cTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCcatcat                           <  2:1109628/49‑1 (MQ=38)
                 cTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCatcatc                          <  1:491897/50‑1 (MQ=39)
                 cTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCatcatc                          <  1:1037585/50‑1 (MQ=39)
                    caaacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCAtcatca                         <  2:131714/48‑1 (MQ=38)
                    caaacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCAtcatca                         <  2:1061545/48‑1 (MQ=38)
                    caaacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCaa                        >  1:682425/1‑49 (MQ=38)
                    caaacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCACCaa                        >  1:1561148/1‑49 (MQ=38)
                     aaacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAAt                       <  1:775296/49‑1 (MQ=38)
                     aaacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAAt                       <  2:2016740/49‑1 (MQ=38)
                     aaacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAAt                       <  2:1907338/49‑1 (MQ=38)
                     aaacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATa                      >  2:1647838/1‑50 (MQ=39)
                     aaacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATa                      >  2:1791541/1‑50 (MQ=39)
                      aacaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATAc                     <  2:610271/50‑1 (MQ=39)
                       acaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATAc                     >  1:28938/1‑49 (MQ=38)
                       acaaaAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACg                    <  1:75009/50‑1 (MQ=39)
                         aaaaaaaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGAt                  >  1:689885/1‑50 (MQ=39)
                         aaaaaaaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGAt                  <  2:1706968/50‑1 (MQ=39)
                         aaaaaaaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGAt                  <  2:1824585/50‑1 (MQ=39)
                          aaaaaaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATg                 <  2:626042/50‑1 (MQ=39)
                           aaaaaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATg                 >  2:2081770/1‑49 (MQ=39)
                           aaaaaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGa                <  2:343077/50‑1 (MQ=255)
                            aaaaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGaa               <  2:1241771/50‑1 (MQ=255)
                            aaaaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGa                >  2:938076/1‑49 (MQ=39)
                             aaaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGa                >  1:1922961/1‑48 (MQ=39)
                              aaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGa                >  1:1186375/1‑47 (MQ=39)
                              aaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAg             <  1:2055819/50‑1 (MQ=255)
                              aaCCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAg             <  1:1140090/50‑1 (MQ=255)
                               accaccTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGa            <  2:1263916/50‑1 (MQ=255)
                                ccaccTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGAc           <  2:1415801/50‑1 (MQ=255)
                                     tGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGAc       >  1:885817/1‑49 (MQ=39)
                                      gCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGAc       >  2:575949/1‑48 (MQ=39)
                                       ccAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGAc       >  2:1557647/1‑47 (MQ=39)
                                       ccAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACaa     <  2:1652045/49‑1 (MQ=39)
                                       ccAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACa      <  2:1686783/48‑1 (MQ=39)
                                       ccAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACa      <  2:1654669/48‑1 (MQ=39)
                                       ccAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACa      <  2:297310/48‑1 (MQ=39)
                                         aGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGAc  <  1:338894/50‑1 (MQ=255)
                                                |                                           
CATTATGCTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGAC  >  NC_017635/2331406‑2331497

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGCTCAAGAC  >  NC_017635/2331413‑2331509
                                                |                                                
CTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCG                                               >  M01420:85:000000000‑AC936:1:1103:21896:6967/1‑51 (MQ=60)
CTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCG                                               >  M01420:85:000000000‑AC936:1:2104:10295:11314/1‑51 (MQ=60)
   TAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCC                                             <  M01420:85:000000000‑AC936:1:2111:5225:14817/50‑1 (MQ=60)
         GCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATC                                      >  M01420:85:000000000‑AC936:1:2109:15635:7804/1‑51 (MQ=60)
         GCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATC                                      <  M01420:85:000000000‑AC936:1:1107:21974:18033/51‑1 (MQ=60)
         GCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATC                                      <  M01420:85:000000000‑AC936:1:2101:9459:4277/51‑1 (MQ=60)
             CAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAAT                                   >  M01420:85:000000000‑AC936:1:1110:12031:9177/1‑50 (MQ=60)
             CAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAAT                                   <  M01420:85:000000000‑AC936:1:1111:16951:15241/50‑1 (MQ=60)
             CAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCACCAAT                                   >  M01420:85:000000000‑AC936:1:2108:18777:7081/1‑50 (MQ=54)
               AACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACG                                <  M01420:85:000000000‑AC936:1:1102:12703:2359/51‑1 (MQ=60)
                ACAAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACG                                >  M01420:85:000000000‑AC936:1:1101:4446:12919/1‑50 (MQ=60)
                  AAAAAAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATG                             >  M01420:85:000000000‑AC936:1:1110:5720:11512/1‑51 (MQ=60)
                      AAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAG                         <  M01420:85:000000000‑AC936:1:2114:4131:18751/51‑1 (MQ=60)
                      AAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAG                         <  M01420:85:000000000‑AC936:1:2102:2956:14556/51‑1 (MQ=60)
                      AAACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAA                           >  M01420:85:000000000‑AC936:1:2112:10326:26263/1‑49 (MQ=60)
                       AACCACCTGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAA                           >  M01420:85:000000000‑AC936:1:2103:20403:5587/1‑48 (MQ=60)
                              TGCCAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACA                  >  M01420:85:000000000‑AC936:1:1113:14004:1718/1‑50 (MQ=60)
                                 CAGGGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGAC              <  M01420:85:000000000‑AC936:1:1105:7479:16579/51‑1 (MQ=60)
                                    GGTGGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTC           <  M01420:85:000000000‑AC936:1:1102:26913:15048/51‑1 (MQ=60)
                                        GGTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGC         <  M01420:85:000000000‑AC936:1:2107:4583:8582/50‑1 (MQ=60)
                                         GTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGCTC       >  M01420:85:000000000‑AC936:1:1111:14389:11011/1‑51 (MQ=60)
                                             TTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGCTCAAG    <  M01420:85:000000000‑AC936:1:2101:24401:6099/50‑1 (MQ=60)
                                             TTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGCTCAAG    <  M01420:85:000000000‑AC936:1:1102:18018:4388/50‑1 (MQ=60)
                                             TTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGCTCA      >  M01420:85:000000000‑AC936:1:1106:27204:18673/1‑48 (MQ=60)
                                             TTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGCTCA      >  M01420:85:000000000‑AC936:1:1103:11953:13332/1‑48 (MQ=60)
                                              TTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGCTCAAGAC  >  M01420:85:000000000‑AC936:1:1112:7587:20398/1‑51 (MQ=60)
                                               TTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGCTCAAGAC  <  M01420:85:000000000‑AC936:1:2104:23516:5521/50‑1 (MQ=60)
                                                |                                                
CTTTAATATGCTTCAAACAAAAAAACCACCTGCCAGGGTGGTTTTTTTTGCGCCCATCATCAATACGATGAAAGACTGACAAGACTTCGCTCAAGAC  >  NC_017635/2331413‑2331509

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