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 3,655,775 G→A S431F (TCT→TTT)  secY ← preprotein translocase subunit SecY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176353,655,7750GA100.0% 226.2 / NA 68S431F (TCT→TTT) secYpreprotein translocase subunit SecY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/56);  total (12/56)

CGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAA  >  NC_017635/3655728‑3655820
                                               |                                             
cGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAAc                                             <  1:1952464/50‑1 (MQ=255)
cGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAAc                                             <  2:2000178/50‑1 (MQ=255)
cGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAAc                                             <  2:1860164/50‑1 (MQ=255)
 gACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACt                                            <  1:1010827/50‑1 (MQ=255)
 gACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACt                                            >  2:144920/1‑50 (MQ=255)
  aCCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTc                                           <  1:273713/50‑1 (MQ=255)
   ccAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCa                                          <  1:1242117/50‑1 (MQ=255)
   ccAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCa                                          <  1:1253501/50‑1 (MQ=255)
   ccAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCa                                          <  1:942855/50‑1 (MQ=255)
    cAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCa                                          >  2:740377/1‑49 (MQ=255)
    cAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCAt                                         <  2:1511250/50‑1 (MQ=255)
      aTTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATac                                       <  1:927148/50‑1 (MQ=255)
       ttATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATact                                      <  2:727383/50‑1 (MQ=255)
       ttATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATact                                      <  2:1906565/50‑1 (MQ=255)
       ttATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATact                                      <  1:291388/50‑1 (MQ=255)
        tATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATact                                      >  2:1091852/1‑49 (MQ=255)
         aTCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATactg                                     <  1:836253/49‑1 (MQ=255)
          tCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATactgac                                   <  2:566232/50‑1 (MQ=255)
           cGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATactgact                                  >  2:703640/1‑50 (MQ=255)
           cGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATactgact                                  >  1:48328/1‑50 (MQ=255)
            ggCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATactgactg                                 <  2:1250162/50‑1 (MQ=255)
            ggCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATactgactg                                 >  1:362761/1‑50 (MQ=255)
            ggCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATactgactg                                 <  1:7694/50‑1 (MQ=255)
            ggCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATactgactg                                 <  1:1062486/50‑1 (MQ=255)
            ggCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATactgactg                                 <  1:1436063/50‑1 (MQ=255)
              ccGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGa                               <  1:28792/50‑1 (MQ=255)
              ccGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGa                               >  1:1561333/1‑50 (MQ=255)
               cGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGa                               >  1:1965886/1‑49 (MQ=255)
               cGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACCGACTGGAc                              <  1:1129199/50‑1 (MQ=39)
                gTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGAca                             <  2:357557/50‑1 (MQ=255)
                gTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGAca                             <  1:1332168/50‑1 (MQ=255)
                  aGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGAcat                            >  2:1704626/1‑49 (MQ=255)
                   gCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGAcatca                          <  2:932628/50‑1 (MQ=255)
                        tCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACAtcatca                       <  1:1605316/48‑1 (MQ=255)
                        tCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACAtcatca                       <  2:305215/48‑1 (MQ=255)
                        tCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACAtcatca                       <  1:1130845/48‑1 (MQ=255)
                         cAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCaga                     <  2:141551/49‑1 (MQ=255)
                          aGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGt                   <  1:336790/50‑1 (MQ=255)
                           ggTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGtt                  <  2:253145/50‑1 (MQ=255)
                            gTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGttt                 <  2:350461/50‑1 (MQ=255)
                             ttCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTg                <  1:1020845/50‑1 (MQ=255)
                             ttCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTg                <  2:1612839/50‑1 (MQ=255)
                               cGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGc               <  1:754398/49‑1 (MQ=255)
                               cGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGc               <  2:1734426/49‑1 (MQ=255)
                               cGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCa              >  1:764267/1‑50 (MQ=255)
                                gCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCAc             <  1:877906/50‑1 (MQ=255)
                                 ccTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACt            >  2:584188/1‑50 (MQ=255)
                                 ccTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACt            <  2:1401464/50‑1 (MQ=255)
                                  cttcttCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACtt           <  2:496992/50‑1 (MQ=255)
                                  cttcttCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACtt           <  1:1866373/50‑1 (MQ=255)
                                  cttcttCAACGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACtt           <  1:1489469/50‑1 (MQ=255)
                                     cttcAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTg          <  1:715820/48‑1 (MQ=255)
                                     cttcAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTg          <  2:1985388/48‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAg        <  1:595645/47‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAg        <  1:1296488/47‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGcc      <  2:2025231/49‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGcc      <  1:628981/49‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGcc      <  1:1838153/49‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGcc      <  2:1643361/49‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGcc      <  1:1469246/49‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGcc      <  2:1219185/49‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGcc      <  1:1386416/49‑1 (MQ=255)
                                        cAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGcc      <  2:1405021/49‑1 (MQ=255)
                                         aaTGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCa     <  1:1936429/49‑1 (MQ=255)
                                          aTGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCAt    >  2:1976742/1‑49 (MQ=255)
                                           tGCAAACTCATACTGACTGGACATCATCAGAGTTTGCat             <  1:817043/39‑2 (MQ=39)
                                           tGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATaa  <  1:1426669/50‑1 (MQ=255)
                                           tGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATaa  <  1:1341011/50‑1 (MQ=255)
                                               |                                             
CGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAA  >  NC_017635/3655728‑3655820

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
GGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCC  >  NC_017635/3655726‑3655825
                                                 |                                                  
GGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAA                                                    >  M01420:85:000000000‑AC936:1:2110:20627:22204/1‑50 (MQ=60)
 GCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAAC                                                  <  M01420:85:000000000‑AC936:1:2113:13180:11458/51‑1 (MQ=60)
  CGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACT                                                 <  M01420:85:000000000‑AC936:1:1114:24642:19571/51‑1 (MQ=60)
   GACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTC                                                <  M01420:85:000000000‑AC936:1:1104:3555:19333/51‑1 (MQ=60)
    ACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCA                                               <  M01420:85:000000000‑AC936:1:2103:23833:23262/51‑1 (MQ=60)
    ACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCA                                               <  M01420:85:000000000‑AC936:1:2103:14392:28965/51‑1 (MQ=60)
    ACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCA                                               <  M01420:85:000000000‑AC936:1:1113:26420:21445/51‑1 (MQ=60)
       AATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATAC                                            <  M01420:85:000000000‑AC936:1:1113:25281:16603/51‑1 (MQ=60)
        ATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACT                                           <  M01420:85:000000000‑AC936:1:1104:7403:25491/51‑1 (MQ=60)
          TATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTG                                          <  M01420:85:000000000‑AC936:1:1112:28088:11473/50‑1 (MQ=60)
             CGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTG                                      >  M01420:85:000000000‑AC936:1:1101:15204:18655/1‑51 (MQ=60)
             CGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTG                                      <  M01420:85:000000000‑AC936:1:2106:15115:14299/51‑1 (MQ=60)
             CGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTG                                      <  M01420:85:000000000‑AC936:1:2101:4827:13098/51‑1 (MQ=60)
             CGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTG                                      <  M01420:85:000000000‑AC936:1:1101:21461:5922/51‑1 (MQ=60)
              GGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGG                                     >  M01420:85:000000000‑AC936:1:1105:16420:24262/1‑51 (MQ=60)
               GCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGA                                    <  M01420:85:000000000‑AC936:1:1101:10521:12879/51‑1 (MQ=60)
                CCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGAC                                   >  M01420:85:000000000‑AC936:1:2108:25868:7159/1‑51 (MQ=60)
                CCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACCGACTGGAC                                   <  M01420:85:000000000‑AC936:1:2102:28536:11213/51‑1 (MQ=60)
                 CGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACA                                  <  M01420:85:000000000‑AC936:1:2105:12744:4391/51‑1 (MQ=60)
                 CGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGAC                                   >  M01420:85:000000000‑AC936:1:2113:19730:15348/1‑50 (MQ=60)
                         TTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCA                            <  M01420:85:000000000‑AC936:1:2102:10887:11738/49‑1 (MQ=60)
                         TTCAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCA                            <  M01420:85:000000000‑AC936:1:2108:12438:20663/49‑1 (MQ=60)
                           CAGGTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGT                        <  M01420:85:000000000‑AC936:1:1105:22936:15977/51‑1 (MQ=60)
                              GTTCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTG                     <  M01420:85:000000000‑AC936:1:1114:6799:22836/51‑1 (MQ=60)
                                TCGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGC                    <  M01420:85:000000000‑AC936:1:1111:28056:8746/50‑1 (MQ=60)
                                 CGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCAC                  <  M01420:85:000000000‑AC936:1:1112:23345:24679/51‑1 (MQ=60)
                                 CGCCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCAC                  >  M01420:85:000000000‑AC936:1:1111:23021:11899/1‑51 (MQ=60)
                                   CCTTCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTT                <  M01420:85:000000000‑AC936:1:2112:8150:8602/51‑1 (MQ=60)
                                   CCTTCTTCAACGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTT                <  M01420:85:000000000‑AC936:1:2107:5983:7849/51‑1 (MQ=60)
                                      TCTTCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTG               <  M01420:85:000000000‑AC936:1:1110:10900:19274/49‑1 (MQ=60)
                                         TCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCC           <  M01420:85:000000000‑AC936:1:1109:5051:15176/50‑1 (MQ=60)
                                         TCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCC           <  M01420:85:000000000‑AC936:1:2111:17176:23009/50‑1 (MQ=60)
                                         TCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCC           <  M01420:85:000000000‑AC936:1:2105:16455:21787/50‑1 (MQ=60)
                                         TCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCC           <  M01420:85:000000000‑AC936:1:2106:15502:24916/50‑1 (MQ=60)
                                         TCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAG             <  M01420:85:000000000‑AC936:1:1109:11028:3803/48‑1 (MQ=60)
                                         TCAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAG             <  M01420:85:000000000‑AC936:1:2104:14291:16975/48‑1 (MQ=60)
                                          CAATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCA          <  M01420:85:000000000‑AC936:1:2113:10605:6207/50‑1 (MQ=60)
                                            ATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCAT                  <  M01420:85:000000000‑AC936:1:1112:21182:4643/40‑1 (MQ=60)
                                            ATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAA       <  M01420:85:000000000‑AC936:1:2106:3515:11483/51‑1 (MQ=60)
                                            ATGCAAACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAA       <  M01420:85:000000000‑AC936:1:2105:13370:7863/51‑1 (MQ=60)
                                                 AACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCC  >  M01420:85:000000000‑AC936:1:1101:16153:23386/1‑51 (MQ=60)
                                                 AACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCC  >  M01420:85:000000000‑AC936:1:2101:22370:26513/1‑51 (MQ=60)
                                                 |                                                  
GGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCC  >  NC_017635/3655726‑3655825

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