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 F2 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 568,752 A→G G71G (GGA→GGG ybcO → DUF1364 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273568,7520AG97.6% 129.8 / ‑3.8 42G71G (GGA→GGGybcODUF1364 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (15/26);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

TACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568665‑568838
                                                                                       |                                                                                      
tACGAACCCGCCAGCCCTGATTGCCACCATTGCATGTTCCGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGgta                                                                                      <  1:314517/90‑1 (MQ=255)
tACGAAACCGCCAGACCTGATTGCCCCCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGgta                                                                                      <  1:672364/90‑1 (MQ=255)
     aaCCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCaa                                                                                 >  2:317343/1‑90 (MQ=255)
         gCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGa                                                                             >  2:624829/1‑90 (MQ=255)
         gCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGa                                                                             >  2:603378/1‑90 (MQ=255)
         gCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGa                                                                             >  1:800475/1‑90 (MQ=255)
         gCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGa                                                                             >  2:44395/1‑90 (MQ=255)
             gACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgc                                                                         <  1:719161/90‑1 (MQ=255)
             gACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgc                                                                         <  1:501007/90‑1 (MQ=255)
               ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  2:194049/90‑1 (MQ=255)
                  gATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCt                                                                    >  2:390498/1‑90 (MQ=255)
                     tGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGa                                                                 <  2:273767/90‑1 (MQ=255)
                     tGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGa                                                                 <  1:197603/90‑1 (MQ=255)
                            aTTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATg                                                          <  1:596733/90‑1 (MQ=255)
                            aTTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATg                                                          >  1:800747/1‑90 (MQ=255)
                                 aTGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCgag                                                     >  1:433843/1‑90 (MQ=255)
                                 aTGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCgag                                                     <  2:450868/90‑1 (MQ=255)
                                      ctgcctgcCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAacac                                                <  2:109639/90‑1 (MQ=255)
                                      ctgcctgcCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAacac                                                <  1:704551/90‑1 (MQ=255)
                                        gcctgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAg                                              >  2:772100/1‑90 (MQ=255)
                                        gcctgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAg                                              >  2:51784/1‑90 (MQ=255)
                                              cACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATc                                        <  1:772100/90‑1 (MQ=255)
                                              cACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATc                                        <  1:603378/90‑1 (MQ=255)
                                              cACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATc                                        <  2:223884/90‑1 (MQ=255)
                                                        tCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAag                              <  1:624829/90‑1 (MQ=255)
                                                         cGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAaga                             <  1:404373/90‑1 (MQ=255)
                                                          gACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAagag                            <  2:53395/90‑1 (MQ=255)
                                                           aCCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAgg                           <  1:390498/90‑1 (MQ=255)
                                                           aCCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTGATCTGGCTGAAAGAgg                           >  1:695356/1‑90 (MQ=255)
                                                               ccgcACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       >  1:174435/1‑58 (MQ=255)
                                                               ccgcACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       <  2:174435/58‑1 (MQ=255)
                                                               ccgcACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGt                       <  1:455147/90‑1 (MQ=255)
                                                                     gcaTTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTaa                 <  2:710792/90‑1 (MQ=255)
                                                                     gcaTTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTaa                 <  2:714582/90‑1 (MQ=255)
                                                                        ttttGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGc              <  2:199451/90‑1 (MQ=255)
                                                                          ttgttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGt            <  1:51784/90‑1 (MQ=255)
                                                                          ttgttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGt            <  1:677917/90‑1 (MQ=255)
                                                                             ttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGaa         <  2:800475/90‑1 (MQ=255)
                                                                              tgACGCTGGGTATGCAAAAGAAAGGGCGCTGGAAGGTATGGCGAGAACCCCGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAAt        <  2:511519/90‑1 (MQ=255)
                                                                               gACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATa       >  2:732988/1‑90 (MQ=255)
                                                                                aCGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATAc      >  1:675328/1‑90 (MQ=255)
                                                                                aCGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAAAACAGGTTATCTGGCTGAAAGGGGGGGTTATTAAGGCGTGAATAc      >  2:372590/1‑90 (MQ=255)
                                                                                  gCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCt    <  2:466870/90‑1 (MQ=255)
                                                                                    tGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAc  <  1:295259/90‑1 (MQ=255)
                                                                                       |                                                                                      
TACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568665‑568838

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ 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
GTGCGGCACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTAC  >  NZ_CP009273/568653‑568851
                                                                                                   |                                                                                                   
GTGCGGCACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGG                                                                                                     >  SRR3722079.201175/1‑100 (MQ=60)
       ACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAA                                                                                              <  SRR3722079.343555/100‑1 (MQ=60)
           GTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGA                                                                                          >  SRR3722079.809051/1‑100 (MQ=60)
            TACGAACCCGCCAGCCCTGATTGCCACCATTGCATGTTCCGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAA                                                                                         <  SRR3722079.317059/100‑1 (MQ=60)
            TACGAAACCGCCAGACCTGATTGCCCCCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAA                                                                                         <  SRR3722079.679397/100‑1 (MQ=60)
                         GACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAG                                                                            <  SRR3722079.726841/100‑1 (MQ=60)
                         GACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAG                                                                            <  SRR3722079.505505/100‑1 (MQ=60)
                              GATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATG                                                                       >  SRR3722079.809326/1‑100 (MQ=60)
                                 TGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCG                                                                    <  SRR3722079.199306/100‑1 (MQ=60)
                                   CCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAG                                                                  >  SRR3722079.437600/1‑100 (MQ=60)
                                        ATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACAC                                                             <  SRR3722079.602714/100‑1 (MQ=60)
                                                  CTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTG                                                   <  SRR3722079.712027/100‑1 (MQ=60)
                                                          CACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAG                                           <  SRR3722079.780366/100‑1 (MQ=60)
                                                          CACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAG                                           <  SRR3722079.609466/100‑1 (MQ=60)
                                                             GACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTGATCTGGCTGAAAGAGG                                        >  SRR3722079.702675/1‑100 (MQ=60)
                                                                 AAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGctgtctcttatacacatctccg                                    >  SRR3722079.175951/1‑78 (MQ=60)
                                                                    TCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTAT                                 <  SRR3722079.631199/100‑1 (MQ=60)
                                                                     CGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATT                                <  SRR3722079.407797/100‑1 (MQ=60)
                                                                       ACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAA                              <  SRR3722079.393736/100‑1 (MQ=60)
                                                                           CCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCG                          <  SRR3722079.459129/100‑1 (MQ=60)
                                                                                  CATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATAC                   >  SRR3722079.682413/1‑100 (MQ=60)
                                                                                      TTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC               <  SRR3722079.685032/100‑1 (MQ=60)
                                                                                      TTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC               <  SRR3722079.52192/100‑1 (MQ=60)
                                                                                                TGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACAT     <  SRR3722079.297684/100‑1 (MQ=60)
                                                                                                 GGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATT    >  SRR3722079.641630/1‑100 (MQ=60)
                                                                                                   GTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTAC  >  SRR3722079.559223/1‑100 (MQ=60)
                                                                                                   |                                                                                                   
GTGCGGCACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTAC  >  NZ_CP009273/568653‑568851

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: