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 146 27.9 2327082 89.7% 2087392 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,602,642 +A intergenic (‑39/+87) h16_A3328 ← / ← h16_A3329 conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,602,6411.A100.0% 37.1 / NA 11intergenic (‑38/+88)h16_A3328/h16_A3329conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (6/5);  total (6/5)

GCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAG  >  AM260479/3602520‑3602727
                                                                                                                          |                                                                                      
gCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaatac                                                                                   >  1:1117189/1‑123 (MQ=255)
gCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaatac                                                                                   <  2:1117189/127‑5 (MQ=255)
        acaGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGa                                                                     >  2:361538/1‑134 (MQ=255)
                  cATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaatac                                                                                   >  1:175768/1‑105 (MQ=255)
                  cATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaatac                                                                                   <  2:175768/109‑5 (MQ=255)
                   aTCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTa                                                          >  1:1027395/1‑134 (MQ=255)
                          ccaccaTCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATtgcatg                                                   <  2:841418/134‑1 (MQ=255)
                            accaTCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATtgcatgca                                                 <  2:1027395/134‑1 (MQ=255)
                                          cAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGc                                  >  2:387279/1‑134 (MQ=255)
                                              tCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGc                              <  1:361538/134‑1 (MQ=255)
                                                             gcgATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTgg                >  2:90399/1‑133 (MQ=255)
                                                                                           gCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg  >  1:1154903/1‑117 (MQ=255)
                                                                                           gCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg  <  2:1154903/117‑1 (MQ=255)
                                                                                                     gCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGgt                    >  1:248326/1‑89 (MQ=255)
                                                                                                     gCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGgt                    <  2:248326/89‑1 (MQ=255)
                                                                                                                          |                                                                                      
GCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAG  >  AM260479/3602520‑3602727

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGGCTGGGCAGGTCCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTG  >  AM260479/3602500‑3602786
                                                                                                                                               |                                                                                                                                                
GGGCTGGGCAGGTCCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGG                                                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1203:9258:91101/150‑1 (MQ=60)
  GCTGGGCAGGTCCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAA                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2112:10904:6844/150‑1 (MQ=60)
     GGGCAGGTCCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGT                                                                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:2214:13788:83503/1‑150 (MQ=60)
             agattagcgagtgtctcgtgggctcggagatgtgtataagagacagGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGAT                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1113:1915:88586/104‑1 (MQ=60)
                       GGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGActgtctcttat                                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1110:5853:31219/1‑139 (MQ=60)
                        GTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAT                                                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:2211:7611:29697/1‑150 (MQ=60)
                                                         GTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAActgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccg                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:1213:17285:78620/1‑95 (MQ=60)
                                                                 GTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGC                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1204:12338:59310/150‑1 (MQ=60)
                                                                                                CCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCctg                                           >  GWNJ‑0478:712:GW2002102894th:2:2216:4294:59548/1‑147 (MQ=60)
                                                                                                          CGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACctgtctcttatacacatctccgagcccacg                                  >  GWNJ‑0478:712:GW2002102894th:2:1115:5847:101178/1‑119 (MQ=60)
                                                                                                                           CCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCctgtctcttatacacatctccgagcccacg                >  GWNJ‑0478:712:GW2002102894th:2:1216:5555:34808/1‑120 (MQ=60)
                                                                                                                              CGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCG              >  GWNJ‑0478:712:GW2002102894th:2:2215:13036:50006/1‑149 (MQ=60)
                                                                                                                               GGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCG              >  GWNJ‑0478:712:GW2002102894th:2:1211:5815:87755/1‑148 (MQ=60)
                                                                                                                                GTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaacagaagccaaataaaaccaaaaagc            >  GWNJ‑0478:712:GW2002102894th:2:1104:9206:46803/1‑49 (MQ=60)
                                                                                                                                         GTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTG  >  GWNJ‑0478:712:GW2002102894th:2:2205:16659:79423/1‑150 (MQ=60)
                                                                                                                                               |                                                                                                                                                
GGGCTGGGCAGGTCCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTG  >  AM260479/3602500‑3602786

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