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 1,011,464 T→C I155V (ATC→GTC)  fabA ← bifunctional 3‑hydroxydecanoyl‑ACP dehydratase/trans‑2‑decenoyl‑ACP isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,011,4640TC100.0% 173.7 / NA 46I155V (ATC→GTC) fabAbifunctional 3‑hydroxydecanoyl‑ACP dehydratase/trans‑2‑decenoyl‑ACP isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/24);  total (22/24)

TCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGA  >  NZ_CP009273/1011379‑1011540
                                                                                     |                                                                            
tCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAgac                                                                         <  1:1114068/91‑1 (MQ=255)
      tttGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCAt                                                                   <  1:123716/91‑1 (MQ=255)
      tttGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCAt                                                                   <  2:1095920/91‑1 (MQ=255)
      tttGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCAt                                                                   <  1:141396/91‑1 (MQ=255)
           tACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAAcc                                                              <  2:1075001/91‑1 (MQ=255)
             ctctGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAg                                                            >  2:310047/1‑91 (MQ=255)
                 gtctgATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACt                                                        >  2:229253/1‑91 (MQ=255)
                       ttATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCa                                                  >  2:187017/1‑91 (MQ=255)
                        tataATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCAt                                                 >  1:774507/1‑91 (MQ=255)
                        tataATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCAt                                                 <  1:388504/91‑1 (MQ=255)
                           aaTCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCg                                              >  1:797735/1‑91 (MQ=255)
                           aaTCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCg                                              >  1:719250/1‑91 (MQ=255)
                               aGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAg                                          <  1:353922/91‑1 (MQ=255)
                               aGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAg                                          >  2:74409/1‑91 (MQ=255)
                                gAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg                                         >  2:144436/1‑91 (MQ=255)
                                gAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg                                         >  1:69429/1‑91 (MQ=255)
                                gAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg                                         >  1:69428/1‑91 (MQ=255)
                                gAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg                                         >  1:210630/1‑91 (MQ=255)
                                   ggCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGccc                                      <  1:24142/91‑1 (MQ=255)
                                     cAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCAt                                    <  1:368357/91‑1 (MQ=255)
                                      aGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATa                                   <  2:956590/91‑1 (MQ=255)
                                            aTCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAga                             >  2:1088817/1‑91 (MQ=255)
                                             tCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAgac                            <  1:317390/91‑1 (MQ=255)
                                             tCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAgac                            <  2:760461/91‑1 (MQ=255)
                                             tCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAgac                            <  2:839886/91‑1 (MQ=255)
                                             tCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAgac                            >  2:856207/1‑91 (MQ=255)
                                             tCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAgac                            <  1:187017/91‑1 (MQ=255)
                                             tCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAgac                            <  1:123924/91‑1 (MQ=255)
                                              ccTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAgacg                           <  2:126492/91‑1 (MQ=255)
                                               cTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAgacga                          <  1:1098890/91‑1 (MQ=255)
                                                    aCAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGtt                     >  2:582941/1‑91 (MQ=255)
                                                       gACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGccc                                      >  2:594349/1‑71 (MQ=255)
                                                       gACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGccc                                      <  1:594349/71‑1 (MQ=255)
                                                       gACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAAc                  >  2:137126/1‑91 (MQ=255)
                                                       gACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAAc                  >  1:861529/1‑91 (MQ=255)
                                                       gACCGACTTTCAGGTCGCCGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAAc                  >  1:261166/1‑91 (MQ=255)
                                                          cGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAAt               <  2:378161/91‑1 (MQ=255)
                                                          cGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAAt               <  2:702913/91‑1 (MQ=255)
                                                               ttCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGtt          >  2:1068721/1‑91 (MQ=255)
                                                                 cAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTa        <  1:189998/91‑1 (MQ=255)
                                                                   ggTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTaaa      <  1:1068721/91‑1 (MQ=255)
                                                                    gTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAg     <  2:721439/91‑1 (MQ=255)
                                                                      cGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTg   <  2:44830/91‑1 (MQ=255)
                                                                       gCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGa  >  2:1108854/1‑91 (MQ=255)
                                                                       gCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGa  >  1:269949/1‑91 (MQ=255)
                                                                       gCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGa  >  2:872992/1‑91 (MQ=255)
                                                                                     |                                                                            
TCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGA  >  NZ_CP009273/1011379‑1011540

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ 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
CGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATACGGTAGGTCACTTT  >  NZ_CP009273/1011364‑1011557
                                                                                                    |                                                                                             
CGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGAC                                                                                               <  SRR3722125.742220/101‑1 (MQ=60)
      GCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACG                                                                                         <  SRR3722125.941399/101‑1 (MQ=60)
               TCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAAC                                                                                <  SRR3722125.1139212/101‑1 (MQ=60)
                     TTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCAC                                                                          <  SRR3722125.126048/101‑1 (MQ=60)
                     TTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCAC                                                                          <  SRR3722125.144090/101‑1 (MQ=60)
                             TCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCAT                                                                  >  SRR3722125.791338/1‑101 (MQ=60)
                                GTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCG                                                               >  SRR3722125.734766/1‑101 (MQ=60)
                                GTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCG                                                               >  SRR3722125.815154/1‑101 (MQ=60)
                                     ATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGG                                                          >  SRR3722125.214796/1‑101 (MQ=60)
                                     ATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGG                                                          >  SRR3722125.70715/1‑101 (MQ=60)
                                     AATATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGG                                                          >  SRR3722125.70716/1‑101 (MQ=60)
                                       TATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCC                                                        <  SRR3722125.396642/101‑1 (MQ=60)
                                              AGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATC                                                 <  SRR3722125.361321/101‑1 (MQ=60)
                                                  taagagacagTCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGAC                                             <  SRR3722125.607130/91‑1 (MQ=60)
                                                  GGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGAC                                             <  SRR3722125.24565/101‑1 (MQ=60)
                                                    CAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGA                                           <  SRR3722125.376078/101‑1 (MQ=60)
                                                            TCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAGC                                   <  SRR3722125.126259/101‑1 (MQ=60)
                                                            TCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAAC                                   <  SRR3722125.323917/101‑1 (MQ=60)
                                                            TCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAAC                                   >  SRR3722125.880533/1‑101 (MQ=60)
                                                            TCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAAC                                   <  SRR3722125.190647/101‑1 (MQ=60)
                                                            TCCTGGAACAGACCGACTTTCAGGTCGCCGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAAC                                   >  SRR3722125.266417/1‑101 (MQ=60)
                                                              CTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAA                                 <  SRR3722125.1123640/101‑1 (MQ=60)
                                                                         gagatgtgtataagagacaGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAG                      <  SRR3722125.850449/82‑1 (MQ=60)
                                                                            CTTTCAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGA                   >  SRR3722125.275384/1‑101 (MQ=60)
                                                                                CAGGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATAC               <  SRR3722125.193709/101‑1 (MQ=60)
                                                                                  GGTCGCTGGCGGTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATACGG             <  SRR3722125.1092811/101‑1 (MQ=60)
                                                                                             GTATAGACCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATACGGTAGGTCACTTT  >  SRR3722125.345897/1‑101 (MQ=60)
                                                                                                    |                                                                                             
CGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATACGGTAGGTCACTTT  >  NZ_CP009273/1011364‑1011557

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: