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 155 25.8 2285122 88.0% 2010907 104.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 1,166,676 +G intergenic (+150/‑109) h16_B1032 → / → h16_B1033 conserved hypothetical protein/probable extra‑cytoplasmic solute receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604801,166,6741.G90.9% 73.2 / 0.9 22intergenic (+148/‑111)h16_B1032/h16_B1033conserved hypothetical protein/probable extra‑cytoplasmic solute receptor
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (10/10);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

ACTGACGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGAA  >  AM260480/1166542‑1166797
                                                                                                                                     |                                                                                                                           
aCTGACGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAg                                                                                                                            <  2:112312/134‑1 (MQ=255)
    aCGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAgggtt                                                                                                                        >  2:744950/1‑131 (MQ=255)
                   ggACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTgcg                                                                                                          <  1:744950/134‑1 (MQ=255)
                     aCCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGcgcg                                                                                                        <  2:656933/134‑1 (MQ=255)
                          aCATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATc                                                                                                    >  2:457228/1‑133 (MQ=255)
                               tgATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGgtttt                                                                                                                      <  2:956909/109‑2 (MQ=255)
                               tgATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGgtttt                                                                                                                      >  1:956909/1‑108 (MQ=255)
                                       aaGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGc                                                                                      >  2:613660/1‑134 (MQ=255)
                                                         gtgtGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGtt                                                                                                              <  1:659079/92‑1 (MQ=255)
                                                         gtgtGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGtt                                                                                                              >  2:659079/1‑92 (MQ=255)
                                                              ggTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGAcg                                                               >  1:832704/1‑134 (MQ=255)
                                                              ggTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGAcg                                                               >  1:131830/1‑134 (MQ=255)
                                                                        tttCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATc                                                     >  1:925730/1‑134 (MQ=255)
                                                                          tCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATcc                                                    <  2:131830/133‑1 (MQ=255)
                                                                          tCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATcc                                                    <  2:832704/133‑1 (MQ=255)
                                                                                                     ccTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAaga                        <  2:472834/134‑1 (MQ=255)
                                                                                                                    ggCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCaa                              >  2:1054762/1‑113 (MQ=255)
                                                                                                                    ggCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCaa                              <  1:1054762/113‑1 (MQ=255)
                                                                                                                       cGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAg      >  1:68076/1‑134 (MQ=255)
                                                                                                                         aTACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGc                                                                        <  1:1015871/66‑1 (MQ=255)
                                                                                                                         aTACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGc                                                                        >  2:1015871/1‑66 (MQ=255)
                                                                                                                            ccGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGaa  <  2:68076/133‑1 (MQ=255)
                                                                                                                               cgAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAg                                     >  1:378915/1‑95 (MQ=255)
                                                                                                                               cgAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAg                                     <  2:378915/95‑1 (MQ=255)
                                                                                                                                     |                                                                                                                           
ACTGACGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGAA  >  AM260480/1166542‑1166797

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTGGAGGAACTGACGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGAACCAATAACCTATG  >  AM260480/1166534‑1166810
                                                                                                                                               |                                                                                                                                      
ccttgttctttttcttccttcttatttttttttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCC                                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:1109:21081:51630/41‑1 (MQ=60)
   ggagatgtgtataagagacagCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTC                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1211:13964:45910/129‑1 (MQ=60)
                       gggctcggagatgtgtataagagacAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAA                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2107:16042:42896/124‑1 (MQ=60)
                        CCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGctgtctcttat                                                                                                          >  GWNJ‑0478:712:GW2002102894th:2:2206:1735:25485/1‑139 (MQ=60)
                          CGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAAC                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:2110:7548:31337/150‑1 (MQ=60)
                                                       ACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGC                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1109:21111:4812/1‑150 (MQ=60)
                                                       ACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGC                                                                           >  GWNJ‑0478:712:GW2002102894th:2:2114:3268:84950/1‑150 (MQ=60)
                                                            ggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGC                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2211:7526:63223/96‑1 (MQ=60)
                                                                 GTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCC                                                                 >  GWNJ‑0478:712:GW2002102894th:2:2205:4942:47267/1‑150 (MQ=60)
                                                                              agattagcgagtgtctcgtgggctcggagatgtgtataagagacagGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAA                                                    <  GWNJ‑0478:712:GW2002102894th:2:2115:19402:29227/104‑1 (MQ=60)
                                                                                                      gagacagCCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGAC                            <  GWNJ‑0478:712:GW2002102894th:2:2212:5931:67599/143‑1 (MQ=60)
                                                                                                                CAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGC                  >  GWNJ‑0478:712:GW2002102894th:2:1104:13982:32026/1‑150 (MQ=60)
                                                                                                                        GTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGctgtctcttatacacatctccgagc          >  GWNJ‑0478:712:GW2002102894th:2:1207:17918:91570/1‑125 (MQ=60)
                                                                                                                                 ATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaactgataaaaaaa  >  GWNJ‑0478:712:GW2002102894th:2:2105:2945:34043/1‑53 (MQ=60)
                                                                                                                                               |                                                                                                                                      
TTGGAGGAACTGACGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGAACCAATAACCTATG  >  AM260480/1166534‑1166810

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