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 I229 R1 214 17.4 943020 97.5% 919444 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,812,807 (C)7→8 intergenic (‑50/+55) celF ← / ← chbR 6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,812,8001.C94.1% 44.1 / ‑2.5 17intergenic (‑43/+62)celF/chbR6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (9/7);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAA  >  NZ_CP009273/1812718‑1812882
                                                                                   |                                                                                  
tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                             <  1:419027/90‑1 (MQ=255)
tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                             >  2:466014/1‑90 (MQ=255)
tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                             >  2:340054/1‑90 (MQ=255)
 cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                            <  1:26574/90‑2 (MQ=255)
 cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                            <  1:363059/90‑2 (MQ=255)
              gTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAaataat                                                                <  1:111790/90‑1 (MQ=255)
                 acgacTTTTAATTTCTGGCTCATAATTTCCCCCTTCAGTACGGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                             <  2:397423/90‑1 (MQ=255)
                   gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                           >  1:95756/1‑90 (MQ=255)
                     cTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                         >  1:391479/1‑90 (MQ=255)
                      ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                        <  1:58008/90‑1 (MQ=255)
                           aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                   <  1:102495/90‑1 (MQ=255)
                                       tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt                                       <  1:318463/90‑1 (MQ=255)
                                         aTTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGt                                    <  2:320807/90‑1 (MQ=255)
                                             ctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                 <  1:137240/90‑1 (MQ=255)
                                                           gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                   >  1:80468/1‑90 (MQ=255)
                                                            aaTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACt                  <  2:148771/90‑1 (MQ=255)
                                                              tACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGa                >  2:229000/1‑90 (MQ=255)
                                                                cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt              >  1:297105/1‑90 (MQ=255)
                                                                cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt              >  1:129948/1‑90 (MQ=255)
                                                                cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt              >  1:200049/1‑90 (MQ=255)
                                                                  gATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaa            >  2:311153/1‑90 (MQ=255)
                                                                            atatCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCaa  >  1:165901/1‑90 (MQ=255)
                                                                                   |                                                                                  
TCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAA  >  NZ_CP009273/1812718‑1812882

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ 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
TCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAG  >  NZ_CP009273/1812718‑1812902
                                                                                          |                                                                                               
TCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAA                                                                                        <  SRR3722113.423455/100‑1 (MQ=60)
 CCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAAT                                                                                       <  SRR3722113.366794/100‑1 (MQ=60)
 CCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAAT                                                                                       <  SRR3722113.26798/100‑1 (MQ=60)
         CAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGC                                                                               >  SRR3722113.96618/1‑100 (MQ=60)
           ATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAA                                                                             >  SRR3722113.395582/1‑100 (MQ=60)
              GTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAG                                                                          <  SRR3722113.112764/100‑1 (MQ=60)
                      TTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATA                                                                  <  SRR3722113.58512/100‑1 (MQ=60)
                           ATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGA                                                             <  SRR3722113.103413/100‑1 (MQ=60)
                                       TAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTAT                                                 <  SRR3722113.321594/100‑1 (MQ=60)
                                             CTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGAC                                           <  SRR3722113.138387/100‑1 (MQ=60)
                                                 CTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAC                                       >  SRR3722113.81175/1‑100 (MQ=60)
                                                      GTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATT                                  >  SRR3722113.201595/1‑100 (MQ=60)
                                                      GTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATT                                  >  SRR3722113.299895/1‑100 (MQ=60)
                                                      GTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATT                                  >  SRR3722113.131046/1‑100 (MQ=60)
                                                                  GATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAA                      >  SRR3722113.167251/1‑100 (MQ=60)
                                                                              ATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATA          >  SRR3722113.176354/1‑100 (MQ=60)
                                                                                      CCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAG  >  SRR3722113.250858/1‑100 (MQ=60)
                                                                                          |                                                                                               
TCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAG  >  NZ_CP009273/1812718‑1812902

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