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 149 26.7 2337422 87.5% 2045244 105.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,253,193 +C coding (243/258 nt) h16_A2071 → Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,253,1911.C100.0% 80.1 / NA 22T81P (ACC→CCC) h16_A2071Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (10/12);  total (10/12)

GCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACACCCGCCTCCCC  >  AM260479/2253081‑2253292
                                                                                                               |                                                                                                     
gCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGcc                                                                                 <  1:679569/134‑1 (MQ=255)
        cACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                           <  1:622274/132‑1 (MQ=255)
                        ccgccgGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGAt                                                                            <  1:345548/115‑1 (MQ=255)
                        ccgccgGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGAt                                                                            >  2:345548/1‑115 (MQ=255)
                          gccgGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTAcgc                                                       >  2:78372/1‑134 (MQ=255)
                                   gTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                           >  1:288067/1‑105 (MQ=255)
                                   gTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                           <  2:288067/105‑1 (MQ=255)
                                                gcCTGTGCATTTCCTCCGATAACACCAACGGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAAccc                                                                                                    <  1:557605/66‑2 (MQ=255)
                                                gcCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAAccc                                                                                                    >  2:557605/1‑65 (MQ=255)
                                                       cATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGt                           <  1:78372/133‑1 (MQ=255)
                                                                cGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAg                                               <  1:1157527/104‑1 (MQ=255)
                                                                cGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAg                                               >  2:1157527/1‑104 (MQ=255)
                                                                   tAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAAccc                                                                                                    <  2:881788/47‑2 (MQ=255)
                                                                   tAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAAccc                                                                                                    >  1:881788/1‑46 (MQ=255)
                                                                              gCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCaa               <  2:433295/122‑1 (MQ=255)
                                                                              gCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCaa               >  1:433295/1‑122 (MQ=255)
                                                                                  gtTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAg                                                                                      >  1:812647/1‑47 (MQ=38)
                                                                                  gtTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAg                                                                                      <  2:812647/47‑1 (MQ=38)
                                                                                   ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                           >  1:86616/1‑57 (MQ=255)
                                                                                   ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                           >  2:418477/1‑57 (MQ=255)
                                                                                   ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                           <  1:418477/57‑1 (MQ=255)
                                                                                   ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATc                                                                           <  2:86616/57‑1 (MQ=255)
                                                                                   ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACACCCGCCTcccc  <  2:1067407/130‑1 (MQ=255)
                                                                                   ttCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACACCCGCCTcccc  >  1:1067407/1‑130 (MQ=255)
                                                                                                     cAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGaa                                <  1:635834/82‑1 (MQ=255)
                                                                                                     cAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGaa                                >  2:635834/1‑82 (MQ=255)
                                                                                                               |                                                                                                     
GCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACACCCGCCTCCCC  >  AM260479/2253081‑2253292

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ 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
TCACCTGGCAGCCACCGCGGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACACCCGCCTCCCCGCATCGCCCGACGCCCTGCCCGCCCATGCCCCGGCAGCGCAC  >  AM260479/2253054‑2253334
                                                                                                                                            |                                                                                                                                             
agattgaggttcttgctgtttttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGC                                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1103:2300:92892/51‑1 (MQ=60)
      ggcatacgagtttagcgagtgtctcgtgggctcggagatgtgtataagagacagACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACACCAACGGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAG                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1216:4618:20327/96‑1 (MQ=60)
                          CGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTAC                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2106:4783:40494/150‑1 (MQ=60)
                               gacagACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGA                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1205:6216:91729/145‑1 (MQ=60)
                                gcagaatacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGAC                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1209:7761:22397/87‑1 (MQ=60)
                                  GCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGAC                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2105:2602:4938/148‑1 (MQ=60)
                                               GATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACctgtctcttatacac                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:1201:15425:25635/1‑135 (MQ=60)
                                                             gtgtataagagacaGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAA                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2216:7627:51930/135‑1 (MQ=60)
                                                                           gagtgtctcgtgggctcggagatgtgtataagagacagGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGC                                                           <  GWNJ‑0478:712:GW2002102894th:2:2105:4649:83322/112‑1 (MQ=60)
                                                                               GTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaa                                                        >  GWNJ‑0478:712:GW2002102894th:2:2201:10454:79854/1‑78 (MQ=60)
                                                                                 GCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACAC                                                      <  GWNJ‑0478:712:GW2002102894th:2:1104:6235:97748/149‑1 (MQ=60)
                                                                                          CCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACACCCGCCTCCCC                                            >  GWNJ‑0478:712:GW2002102894th:2:1211:5378:6414/1‑150 (MQ=60)
                                                                                              TAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaa                                         >  GWNJ‑0478:712:GW2002102894th:2:2115:11921:35373/1‑77 (MQ=60)
                                                                                               AACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgc                                       >  GWNJ‑0478:712:GW2002102894th:2:1107:3249:46123/1‑87 (MQ=60)
                                                                                               AACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACACCCGCCTCCCCGCATC                                       >  GWNJ‑0478:712:GW2002102894th:2:2212:7250:36759/1‑150 (MQ=60)
                                                                                                                                    GGCAAACCCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaatcagccactactagtcctgccc  >  GWNJ‑0478:712:GW2002102894th:2:1108:15898:2640/1‑50 (MQ=60)
                                                                                                                                            |                                                                                                                                             
TCACCTGGCAGCCACCGCGGTCGACCCGCATCCTGCACGTACGCACCGATGCCGCCGGCCACGTCGCGCAGGTGCGCCTGTGCATTTCCTCCGATAACAACAACAGCGTGTTCGCTCCGCAGCCGGTGCAGTGGCAAACCTGCGCCCGCTGGTAGCGGCCGAGATCGCGCAGTACCGAGACTACGCGGCCCGAGCAGGCCGCTGTTGAAGCCGTGACCCCCAGCAACACCCGCCTCCCCGCATCGCCCGACGCCCTGCCCGCCCATGCCCCGGCAGCGCAC  >  AM260479/2253054‑2253334

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