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 152 24.9 2083062 91.6% 1908084 109.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,602,688 (G)6→5 intergenic (‑85/+41) 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,6830G.100.0% 85.6 / NA 20intergenic (‑80/+46)h16_A3328/h16_A3329conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (11/9);  total (11/9)

GTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGC  >  AM260479/3602574‑3602809
                                                                                                              |                                                                                                                              
gTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGc                                                                                                                    >  1:241133/1‑122 (MQ=255)
gTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGc                                                                                                                    <  2:241133/122‑1 (MQ=255)
                   gcgcCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATg                                                                                                                      <  1:355198/101‑1 (MQ=255)
                   gcgcCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATg                                                                                                                      >  2:355198/1‑101 (MQ=255)
                    cgcCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg                                                                                    <  2:791302/134‑1 (MQ=255)
                                         aGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCagcgcg                                                               <  2:678590/134‑1 (MQ=255)
                                                 ccccGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGcgcc                                                       >  1:297227/1‑134 (MQ=255)
                                                   ccGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATg                                                                                                                      >  2:599620/1‑69 (MQ=255)
                                                   ccGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATg                                                                                                                      <  1:599620/69‑1 (MQ=255)
                                                                  cgaaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCtctc                                      >  1:243726/4‑134 (MQ=255)
                                                                     aTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAGAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAg                                   >  1:275275/1‑134 (MQ=255)
                                                                                     ggATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTgc                   >  1:727523/1‑134 (MQ=255)
                                                                                        tCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGcgcc                <  2:727523/134‑1 (MQ=255)
                                                                                              gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                                     >  1:792877/1‑107 (MQ=255)
                                                                                              gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                                     <  2:42455/107‑1 (MQ=255)
                                                                                              gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                                     >  1:42455/1‑107 (MQ=255)
                                                                                              gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                                     <  2:792877/107‑1 (MQ=255)
                                                                                              gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCAGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                                     <  2:363726/107‑1 (MQ=255)
                                                                                              gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCAGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                                     >  1:363726/1‑107 (MQ=255)
                                                                                                        tgcaGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGc  >  2:780922/1‑132 (MQ=255)
                                                                                                              |                                                                                                                              
GTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGC  >  AM260479/3602574‑3602809

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
GGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGCAG  >  AM260479/3602559‑3602824
                                                                                                                                  |                                                                                                                                        
atacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGG                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2107:2071:40476/99‑1 (MQ=60)
agatgtgtataagagacagAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGG                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1206:2474:57218/131‑1 (MQ=60)
GGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGG                                                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:1114:8816:98047/1‑150 (MQ=60)
                                            agcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCG                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1214:7441:99349/87‑1 (MQ=60)
                                                 GTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCG                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1202:17279:78588/1‑150 (MQ=60)
                                                                  CCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCG                                                    >  GWNJ‑0478:712:GW2002102894th:2:1116:18012:12818/1‑150 (MQ=60)
                                                                     GTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAGAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGC                                                 >  GWNJ‑0478:712:GW2002102894th:2:1203:20830:28303/1‑150 (MQ=60)
                                                                                 gggctcggagatgtgtataagagacaGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAAC                                     <  GWNJ‑0478:712:GW2002102894th:2:1214:13890:88903/124‑1 (MQ=60)
                                                                                     TACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCG                                 >  GWNJ‑0478:712:GW2002102894th:2:2115:6468:16920/1‑150 (MQ=60)
                                                                                              GTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACctgtctcttatac                        >  GWNJ‑0478:712:GW2002102894th:2:1105:16749:93193/1‑137 (MQ=60)
                                                                                              GTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACctgtctcttatac                        >  GWNJ‑0478:712:GW2002102894th:2:2203:17471:7519/1‑137 (MQ=60)
                                                                                              GTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCAGGTGTGGCAGCGGCAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGGCCGCTCTCGAGCTTCACCTTGAACctgtctcttatac                        >  GWNJ‑0478:712:GW2002102894th:2:1208:7980:10170/1‑137 (MQ=60)
                                                                                                                 ATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGG     >  GWNJ‑0478:712:GW2002102894th:2:1116:5073:53550/1‑150 (MQ=60)
                                                                                                                    GCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACctgtctcttatacacatctccgagcccacgagacg  >  GWNJ‑0478:712:GW2002102894th:2:2209:13666:12362/1‑115 (MQ=60)
                                                                                                                                  |                                                                                                                                        
GGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGCAG  >  AM260479/3602559‑3602824

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: