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 150 26.0 2053212 90.0% 1847890 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 628,642 Δ1 bp coding (756/1206 nt) h16_A0595 ← Uncharacterized Fe‑S protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479628,6410A.100.0% 80.5 / NA 19coding (757/1206 nt)h16_A0595Uncharacterized Fe‑S protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (10/9);  total (10/9)

CATGCGCGAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAAGCGGCGCGCGTCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGG  >  AM260479/628515‑628765
                                                                                                                              |                                                                                                                             
caTGCGCGAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGgcgc                                                                                                                       <  1:501059/134‑1 (MQ=255)
caTGCGCGAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGgcgc                                                                                                                       >  2:1004004/1‑134 (MQ=255)
                       cAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTc                                                                                                                 <  2:66181/116‑1 (MQ=255)
                       cAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTc                                                                                                                 >  1:66181/1‑116 (MQ=255)
                            cAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATgg                                                                                           >  1:581816/1‑134 (MQ=255)
                                        gtcgCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTgg                                                                                  >  1:193393/1‑131 (MQ=255)
                                        gtcgCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTgg                                                                                  <  2:193393/131‑1 (MQ=255)
                                                                  tGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGt                                                                                                          <  1:143323/81‑1 (MQ=255)
                                                                  tGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGt                                                                                                          >  2:143323/1‑81 (MQ=255)
                                                                                         aTCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAg                                                                     >  1:447282/1‑95 (MQ=255)
                                                                                         aTCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAg                                                                     <  2:447282/95‑1 (MQ=255)
                                                                                              gCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTcgg                         >  1:733110/1‑134 (MQ=255)
                                                                                              gCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTcgg                         >  1:547768/1‑134 (MQ=255)
                                                                                              gCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTcgg                         >  1:33110/1‑134 (MQ=255)
                                                                                                ccTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTcggcc                       >  2:378426/1‑134 (MQ=255)
                                                                                                              gTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGt          <  2:733110/133‑1 (MQ=255)
                                                                                                                   gTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCa    <  2:33110/134‑1 (MQ=255)
                                                                                                                   gTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCa    <  2:547768/134‑1 (MQ=255)
                                                                                                                      ggAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAgg  <  1:1004004/133‑1 (MQ=255)
                                                                                                                              |                                                                                                                             
CATGCGCGAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAAGCGGCGCGCGTCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGG  >  AM260479/628515‑628765

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
ACATGCGCGAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAAGCGGCGCGCGTCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAGGAAG  >  AM260479/628514‑628780
                                                                                                                                |                                                                                                                                           
ACATGCGCGAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTAC                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1216:20957:20791/150‑1 (MQ=60)
         AACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCtgt                                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1104:3654:53203/1‑147 (MQ=60)
             gctccgtctcgtgggctcggagatgtgtataagagacagAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGG                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1111:12883:41475/111‑1 (MQ=60)
              GTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGC                                                                                                         >  GWNJ‑0478:712:GW2002102894th:2:2104:2244:71478/1‑150 (MQ=60)
                          GGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCctgtctcttatacacatctccgagcccacgagacggagc                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:2107:7446:73007/1‑111 (MQ=60)
                          GGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGC                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1112:16238:41251/1‑149 (MQ=60)
                                                                           CGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTctgtctcttatacacatctccgagc                                            >  GWNJ‑0478:712:GW2002102894th:2:1212:12895:21101/1‑125 (MQ=60)
                                                                                CTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGC                                       >  GWNJ‑0478:712:GW2002102894th:2:1105:12104:28462/1‑150 (MQ=60)
                                                                                CTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGC                                       >  GWNJ‑0478:712:GW2002102894th:2:2102:9799:38803/1‑150 (MQ=60)
                                                                                CTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGC                                       >  GWNJ‑0478:712:GW2002102894th:2:2115:9950:99205/1‑150 (MQ=60)
                                                                                                                      AGGAGATGCAGCGGCGCGCGTCCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAGGAAG  <  GWNJ‑0478:712:GW2002102894th:2:2215:20313:67444/149‑1 (MQ=60)
                                                                                                                                |                                                                                                                                           
ACATGCGCGAACTTGTTCCACGGGCAGGCCAGCTGGCAGTCGTCGCAGCCGTAGACGCGGTTGCCCATGGGCGCGCGGAACTCCACCGGAATCGGGCCCTTGTGCTCGATGGTCAGGTAGGAGATGCAAGCGGCGCGCGTCAGCCGGTACGGCGCCACGATGGCGCGCGTGGGGCAGATGTCCAGGCAGCGGTGGCAGTTGCCGCAGTGGCCGGACTCGGGGGCGTCGGCCGGCAGCGGGATGTCGACCAGGATCTCGCCCAGGAAG  >  AM260479/628514‑628780

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