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 1,231,177 Δ1 bp intergenic (+14/‑201) glnE → / → h16_A1128 glutamine synthetase adenylyltransferase/Membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,231,1740G.100.0% 81.2 / NA 19intergenic (+11/‑204)glnE/h16_A1128glutamine synthetase adenylyltransferase/Membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (9/10);  total (9/10)

CGGCGGCATGGTCGATATCGAATTCACGGGTGCAGTACCTGGTGCTGCTGCACAACCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCACCGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGCGCGCGCCAGCACCAGCT  >  AM260479/1231088‑1231285
                                                                                      |                                                                                                               
cggcggCATGGTCGATATCGAATTCACGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGc                                                                                  >  2:774383/1‑115 (MQ=255)
cggcggCATGGTCGATATCGAATTCACGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGc                                                                                  <  1:774383/115‑1 (MQ=255)
 ggcggcATGGTCGATATCGAATTCACGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGc                                                              <  2:646285/134‑1 (MQ=255)
   cggcATGGTCGATATCGAATTCACGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGcca                                                            <  2:227082/134‑3 (MQ=255)
                 tCGAATTCACGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCg                                              >  1:97127/1‑134 (MQ=255)
                         cacggTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGcc                                                                        <  1:299873/98‑1 (MQ=255)
                         cacggTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGcc                                                                        >  2:299873/4‑101 (MQ=255)
                           cggTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGAt                                     <  1:354377/133‑1 (MQ=255)
                             gTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGc                                   <  2:97127/134‑1 (MQ=255)
                                          tgctgctgcACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTc                      <  2:223529/134‑1 (MQ=255)
                                            ctgctgcACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCcg                    >  1:536645/1‑134 (MQ=255)
                                                  gcacagccGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCcgcgcgc               >  1:476511/7‑134 (MQ=255)
                                                    acagccGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCgg                                             <  1:945747/98‑1 (MQ=255)
                                                    acagccGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCgg                                             >  2:945747/5‑102 (MQ=255)
                                                               cACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGCGCGCGCCAGCACCAGCt  >  1:690779/1‑134 (MQ=255)
                                                                 cccGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGc                                                     >  2:125360/1‑81 (MQ=255)
                                                                 cccGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGc                                                     <  1:125360/81‑1 (MQ=255)
                                                                         tGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAg                                                                                       <  1:524726/39‑1 (MQ=37)
                                                                         tGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAg                                                                                       >  2:524726/1‑39 (MQ=37)
                                                                                      |                                                                                                               
CGGCGGCATGGTCGATATCGAATTCACGGGTGCAGTACCTGGTGCTGCTGCACAACCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCACCGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGCGCGCGCCAGCACCAGCT  >  AM260479/1231088‑1231285

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ 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
AAGGTCGCCGAAGGCCATCCCAACCCGACCGACCTGTTCGACCTGAAGCACGATCGCGGCGGCATGGTCGATATCGAATTCACGGGTGCAGTACCTGGTGCTGCTGCACAACCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCACCGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGCGCGCGCCAGCACCAGCTGCGGCTCGAGGGGCAGGCGCATGCGCGCGTCGCGC  >  AM260479/1231032‑1231320
                                                                                                                                                 |                                                                                                                                               
ttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGATATCGAATTCACGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATAT                                                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2203:5242:30570/80‑1 (MQ=60)
              ttttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGC                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:2109:18319:10886/69‑1 (MQ=60)
                                  ttttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCC                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2101:6668:10277/69‑1 (MQ=60)
                                    gacagCCTGAAGCACGATCGCGGCGGCATGGTCGATATCGAATTCACGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATC                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2116:2007:77140/145‑1 (MQ=60)
                                              agatgtgtataagagacagGGTCGATATCGAATTCACGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCcacgc                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1202:19285:96368/131‑6 (MQ=60)
                                                          GCGGCATGGTCGATATCGAATTCACGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGG                                                                                >  GWNJ‑0478:712:GW2002102894th:2:1106:9037:61153/1‑150 (MQ=60)
                                                                   gctccgtctcgtgggctcggagatgtgtataagagacaggcacagCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTaggcat                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1109:17858:38812/105‑7 (MQ=60)
                                                                          gatgtgtataagagacagACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACC                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2210:18969:79339/132‑1 (MQ=60)
                                                                                  acGGTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTC                                                         <  GWNJ‑0478:712:GW2002102894th:2:1206:20406:51894/148‑1 (MQ=60)
                                                                                     GTGCAGTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGC                                                      >  GWNJ‑0478:712:GW2002102894th:2:2101:6963:97154/1‑150 (MQ=60)
                                                                                          GTACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGCGCGCG                                                 >  GWNJ‑0478:712:GW2002102894th:2:1215:16735:47668/1‑150 (MQ=60)
                                                                                            ACCTGGTGCTGCTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttc                                               >  GWNJ‑0478:712:GW2002102894th:2:1113:2933:13636/1‑90 (MQ=60)
                                                                                                       CTGCACAGCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGCGCGCGCCAGCACCAGCTG                                    >  GWNJ‑0478:712:GW2002102894th:2:2113:12849:5448/1‑150 (MQ=60)
                                                                                                                           CGCAGCTGACGCGCAACGCGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGCGCGCGCCAGCACCAGCctgtctcttatacacatctcc                 >  GWNJ‑0478:712:GW2002102894th:2:2212:18323:86167/1‑129 (MQ=60)
                                                                                                                                             gcGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCCACGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGCGCGCGCCAGCACCAGCTGCGGCctgtctcttatacacctctccgagcccacg  >  GWNJ‑0478:712:GW2002102894th:2:1101:1953:11826/3‑118 (MQ=60)
                                                                                                                                                 |                                                                                                                                               
AAGGTCGCCGAAGGCCATCCCAACCCGACCGACCTGTTCGACCTGAAGCACGATCGCGGCGGCATGGTCGATATCGAATTCACGGGTGCAGTACCTGGTGCTGCTGCACAACCCGGGCGCACCCGCAGCTGACGCGCAACGCGGGGCAATATCGCGCTGCTGCGCGCAGCGGGCGAGCTGGGCCTGATCGACGCACCGCTGGCACTGGCGGTGGGCGATGCCTACCGGCTGTTCCGCGCGCGCCAGCACCAGCTGCGGCTCGAGGGGCAGGCGCATGCGCGCGTCGCGC  >  AM260479/1231032‑1231320

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