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,292,091 +G coding (36/39 nt) ECW_RS27980 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176352,292,0911.G100.0% 85.0 / NA 41V12V (GTA→GTCECW_RS27980hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (10/31);  total (10/31)

GCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACC  >  NC_017635/2292045‑2292130
                                               |                                       
gCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTtatg                                         <  1:2427676/48‑2 (MQ=255)
gCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTtatg                                         <  2:2386008/48‑2 (MQ=255)
 cGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTtatga                                        <  1:2261679/48‑3 (MQ=255)
           cccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGcac                              <  2:808384/48‑1 (MQ=17)
           cccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGcac                              <  1:1920437/48‑1 (MQ=17)
               cccTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACg                          <  1:2479118/48‑1 (MQ=17)
                 cTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgat                        <  2:2145027/48‑1 (MQ=17)
                    gTAGTCTTATTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  1:1824611/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatga                      <  2:2115733/47‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  1:1005078/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  2:1391168/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  2:2134369/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  2:2218872/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  2:2408773/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  2:275619/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  2:381574/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  2:574630/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  2:610048/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  2:943499/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  1:1020343/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  1:2197041/46‑1 (MQ=17)
                    gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatg                       <  1:874044/46‑1 (MQ=17)
                     tAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatgat                     >  2:2616223/1‑47 (MQ=17)
                     tAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatgat                     >  1:263201/1‑47 (MQ=17)
                     tAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACgatgat                     >  2:2035247/1‑47 (MQ=17)
                       gTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGAtgatga                   <  2:16121/47‑1 (MQ=17)
                       gTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAg                  <  2:1285544/48‑1 (MQ=17)
                        tCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAg                  >  2:1582430/1‑47 (MQ=17)
                            aGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGc              <  1:1536138/47‑1 (MQ=17)
                             gTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCtt            >  2:1156647/1‑48 (MQ=17)
                              tAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCtt            >  2:1007007/1‑47 (MQ=17)
                              tAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtc           <  1:2567127/48‑1 (MQ=17)
                                aGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtca          <  1:1233865/47‑1 (MQ=17)
                                aGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtca          <  2:2042113/47‑1 (MQ=17)
                                 gCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtcatc        <  1:2514836/48‑1 (MQ=14)
                                 gCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtcatc        <  2:333998/48‑1 (MQ=14)
                                 gCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTtcatc        <  1:433668/48‑1 (MQ=14)
                                     ggAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGa    >  2:214543/1‑48 (MQ=11)
                                     ggAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGa    >  2:1358007/1‑48 (MQ=11)
                                     ggAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGa    >  1:1185055/1‑48 (MQ=11)
                                      gAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGAc   <  2:2202672/48‑1 (MQ=11)
                                      gAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGAc   <  1:985085/48‑1 (MQ=11)
                                       aaGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGAcc  >  1:2132748/1‑48 (MQ=11)
                                       aaGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGAcc  <  2:667906/48‑1 (MQ=11)
                                               |                                       
GCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACC  >  NC_017635/2292045‑2292130

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ 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
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTT  >  NC_017635/2292042‑2292138
                                                 |                                                
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATG                                                <  M01420:84:000000000‑ABVEH:1:2114:10722:12629/51‑1 (MQ=52)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:2108:14927:25489/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:1106:25580:24047/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:1107:12968:7976/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:2112:20620:17734/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:1109:7881:18619/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:1110:10109:19425/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:2110:10311:13258/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:1112:12939:24073/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:2104:6333:7286/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:2101:7526:15311/50‑1 (MQ=53)
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:2109:9340:26353/50‑1 (MQ=53)
AGAGCGGGGTTCCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTAT                                                  <  M01420:84:000000000‑ABVEH:1:2109:22907:18601/50‑1 (MQ=49)
  AGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATgac                                               <  M01420:84:000000000‑ABVEH:1:2107:15997:20613/51‑4 (MQ=60)
   GCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACT                                              <  M01420:84:000000000‑ABVEH:1:2106:7077:11474/51‑1 (MQ=49)
    CGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTA                                             <  M01420:84:000000000‑ABVEH:1:1113:12737:20825/51‑1 (MQ=49)
     GGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTA                                             >  M01420:84:000000000‑ABVEH:1:1114:26368:22831/1‑50 (MQ=49)
             CCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCA                                    <  M01420:84:000000000‑ABVEH:1:2103:19972:12247/51‑1 (MQ=52)
                 GCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGAT                                <  M01420:84:000000000‑ABVEH:1:2108:18494:7827/51‑1 (MQ=60)
                    CTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGAT                             >  M01420:84:000000000‑ABVEH:1:2111:17573:7001/1‑51 (MQ=52)
                    CTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGAT                             >  M01420:84:000000000‑ABVEH:1:2114:3492:9432/1‑51 (MQ=52)
                      GGTAGTCTTATTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGAT                             <  M01420:84:000000000‑ABVEH:1:2102:21566:15878/49‑1 (MQ=60)
                      GGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATG                            >  M01420:84:000000000‑ABVEH:1:1103:19714:10118/1‑50 (MQ=53)
                      GGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGAT                             <  M01420:84:000000000‑ABVEH:1:2105:23203:21348/49‑1 (MQ=60)
                      GGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGAT                             <  M01420:84:000000000‑ABVEH:1:2112:21182:27310/49‑1 (MQ=60)
                      GGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGAT                             <  M01420:84:000000000‑ABVEH:1:2113:18802:13400/49‑1 (MQ=60)
                      GGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGAT                             <  M01420:84:000000000‑ABVEH:1:1109:7404:20704/49‑1 (MQ=60)
                      GGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGAT                             <  M01420:84:000000000‑ABVEH:1:1109:19210:17739/49‑1 (MQ=60)
                      GGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGAT                             <  M01420:84:000000000‑ABVEH:1:1108:16528:15244/49‑1 (MQ=60)
                              TAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTT                    <  M01420:84:000000000‑ABVEH:1:1114:9107:4599/50‑1 (MQ=40)
                                GTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCAT                 <  M01420:84:000000000‑ABVEH:1:2108:21810:25304/51‑1 (MQ=34)
                                GTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCAT                 >  M01420:84:000000000‑ABVEH:1:2110:17775:21428/1‑51 (MQ=34)
                                  AAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATC                <  M01420:84:000000000‑ABVEH:1:1111:27918:16175/50‑1 (MQ=28)
                                  AAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATC                <  M01420:84:000000000‑ABVEH:1:2114:26341:23692/50‑1 (MQ=28)
                                   AGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCAT              <  M01420:84:000000000‑ABVEH:1:2108:13606:14738/51‑1 (MQ=25)
                                   AGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCAT              <  M01420:84:000000000‑ABVEH:1:1104:11717:20598/51‑1 (MQ=25)
                                      GGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGAC           >  M01420:84:000000000‑ABVEH:1:1111:14208:6529/1‑51 (MQ=12)
                                       GGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACC          >  M01420:84:000000000‑ABVEH:1:2113:9472:12292/1‑51 (MQ=9)
                                        GGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCC         >  M01420:84:000000000‑ABVEH:1:2105:17333:8766/1‑51 (MQ=6)
                                        GGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCC         <  M01420:84:000000000‑ABVEH:1:1109:28027:13979/51‑1 (MQ=6)
                                            GCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTC     <  M01420:84:000000000‑ABVEH:1:2112:12688:17243/51‑1 (MQ=0)
                                               TATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTT  >  M01420:84:000000000‑ABVEH:1:2109:17507:8892/1‑51 (MQ=0)
                                                 |                                                
AGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTT  >  NC_017635/2292042‑2292138

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