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 F1 I1 R1 154 26.6 2069536 90.2% 1866721 108.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 951,073 +G intergenic (‑29/+230) h16_B0838 ← / ← h16_B0839 conserved hypothetical protein/transcriptional regulator, CRP‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260480951,0731.G100.0% 74.3 / NA 20intergenic (‑29/+230)h16_B0838/h16_B0839conserved hypothetical protein/transcriptional regulator, CRP‑family
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (10/10);  total (10/10)

CGACGGCATGGCTTGCACGGGCTGGGCGGACTGCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAGGGAACAGTCACTGTAGCGCCG  >  AM260480/950964‑951197
                                                                                                              |                                                                                                                            
cGACGGCATGGCTTGCACGGGCTGGGCGGACTGCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGAcg                                                                                                       >  2:383064/1‑134 (MQ=255)
                  gggCTGGGCGGACTGCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCgccagcca                                                                                     >  2:913120/1‑134 (MQ=255)
                        ggCGGACTGCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGAc                                                                                 >  1:7640/1‑132 (MQ=255)
                                 cACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGAttt                                                                        <  1:978142/132‑1 (MQ=255)
                                 cACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGAttt                                                                        <  1:73905/132‑1 (MQ=255)
                                  aCCGGCTGGCTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCgcca                                                                                         <  2:233264/114‑1 (MQ=255)
                                  aCCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCgcca                                                                                         >  1:233264/1‑114 (MQ=255)
                                                   catcggctgcatcggctgCGTCATCGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTc                                                    <  2:7640/134‑1 (MQ=255)
                                                                       tCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGcc                                 >  1:461992/1‑133 (MQ=255)
                                                                           tGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGAttt                                                                        <  1:410280/90‑1 (MQ=255)
                                                                           tGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGAttt                                                                        >  2:410280/1‑90 (MQ=255)
                                                                             gCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAgg                           <  1:383064/133‑1 (MQ=255)
                                                                                   gctggctCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAgg                      <  2:361587/132‑1 (MQ=255)
                                                                                   gctggctCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAgg                      >  1:361587/1‑132 (MQ=255)
                                                                                                 aGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAGGGAACAGTCACTGTAgc      <  1:913120/134‑1 (MQ=255)
                                                                                                  ggggCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTa                                        >  1:204394/1‑99 (MQ=255)
                                                                                                  ggggCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTa                                        <  2:204394/99‑1 (MQ=255)
                                                                                                     gCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCCCCTAGAACGCCTACAGGCAAGGGAACAGTCACTGTAGcgccg  >  2:326048/1‑134 (MQ=255)
                                                                                                       tGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCt                                                          <  1:830020/76‑1 (MQ=255)
                                                                                                       tGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCt                                                          >  2:830020/1‑76 (MQ=255)
                                                                                                              |                                                                                                                            
CGACGGCATGGCTTGCACGGGCTGGGCGGACTGCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAGGGAACAGTCACTGTAGCGCCG  >  AM260480/950964‑951197

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ 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
GGCTTGCACGGGCTGGGCGGACTGCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAGGGAACAGTCACTGTAGCGCCGCAGCACCTGAC  >  AM260480/950973‑951208
                                                                                                    |                                                                                                                                        
GGCTTGCACGGGCTGGGCGGACTGCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACC                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1101:9441:57362/1‑148 (MQ=60)
          GGCTGGGCGGACTGCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATctgtct                                                                               >  GWNJ‑0478:712:GW2002102894th:2:1114:12863:66960/1‑144 (MQ=60)
                     cgtgggctcggagatgtgtataagagacaGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTT                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1211:17583:99635/121‑1 (MQ=60)
                       GCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTT                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1107:4767:9959/148‑1 (MQ=60)
                       GCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTT                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2213:13105:89846/148‑1 (MQ=60)
                                   atgtagctccgtctcgtgggctcggagatgtgtataagagacaGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTC                                                      <  GWNJ‑0478:712:GW2002102894th:2:2204:3115:90679/107‑1 (MQ=60)
                                               GCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCT                                           >  GWNJ‑0478:712:GW2002102894th:2:1212:4513:92408/1‑149 (MQ=60)
                                                           GCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAGGGA                               >  GWNJ‑0478:712:GW2002102894th:2:1208:5877:11665/1‑149 (MQ=60)
                                                                   GGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAGGGAACAGTCAC                       <  GWNJ‑0478:712:GW2002102894th:2:1209:11414:44231/149‑1 (MQ=60)
                                                                          GCTGGCTCGGCTGCAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCActgtctcttatacacatctcc               >  GWNJ‑0478:712:GW2002102894th:2:1115:3281:4253/1‑129 (MQ=60)
                                                                                       CAGGGGCTGAGCCTGCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAGGGAACAGTCACTGTAGCGCCGCAGCACCTGAC  <  GWNJ‑0478:712:GW2002102894th:2:2211:5791:15531/150‑1 (MQ=60)
                                                                                                    |                                                                                                                                        
GGCTTGCACGGGCTGGGCGGACTGCACCGGCTGACTGGATTGCATCGGCTGCATCGGCTGCGTCATTGGCATCGGCTGGCTCGGCTGCAGGGGCTGAGCCTCTGCGCAAACAGCGGGTAGGACGCGATCATCGCGCCAGCCACGACCAGGGATTTCGCCAGATGTTGCTTCATTCGTGTACCTCCTAGAACGCCTACAGGCAAGGGAACAGTCACTGTAGCGCCGCAGCACCTGAC  >  AM260480/950973‑951208

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