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 154 26.6 2069536 90.2% 1866721 108.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 582,759 G→T C352F (TGT→TTT)  h16_A0555 → zinc metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479582,7590GT100.0% 98.6 / NA 27C352F (TGT→TTT) h16_A0555zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/14);  total (13/14)

GACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCC  >  AM260479/582631‑582884
                                                                                                                                |                                                                                                                             
gACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTgcgc                                                                                                                          >  2:593433/1‑134 (MQ=255)
        aaaaTTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCgg                                                                                                                  >  2:474911/1‑134 (MQ=255)
           aTTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGcc                                                                                                               <  1:237673/134‑1 (MQ=255)
           aTTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGcc                                                                                                               <  1:474911/134‑1 (MQ=255)
           aTTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGcc                                                                                                               <  1:482480/134‑1 (MQ=255)
                      tCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATc                                                                                                      >  1:953835/1‑132 (MQ=255)
                          aTGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGaa                                                                                                >  2:744162/1‑134 (MQ=255)
                                  aCTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGcc                                                                                        <  1:920131/134‑1 (MQ=255)
                                    tCACCCGCTGACCTCCGGGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATg                                                                                                                     <  2:515135/103‑1 (MQ=255)
                                    tCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATg                                                                                                                     >  1:515135/1‑103 (MQ=255)
                                      aCCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCa                                                                                      <  1:587136/132‑1 (MQ=255)
                                          gCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGAcc                                                                                >  1:818592/1‑134 (MQ=255)
                                                  ccGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGcgcg                                                                                                                         <  2:75763/85‑1 (MQ=255)
                                                  ccGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGcgcg                                                                                                                         >  1:75763/1‑85 (MQ=255)
                                                      gCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTgc                                                                    >  1:368555/1‑134 (MQ=255)
                                                             gCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCgcgg                                                              <  2:953835/133‑1 (MQ=255)
                                                                 aCATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTg                                                         <  2:368555/134‑1 (MQ=255)
                                                                        ggACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCCGCGCTGcgcacg                                                   >  1:836068/1‑133 (MQ=255)
                                                                                  aGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCgg                                         >  2:425675/1‑133 (MQ=255)
                                                                                             cAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCAc                             <  1:425675/134‑1 (MQ=255)
                                                                                                     aGGTCGCCATCACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCt                     >  1:738078/1‑134 (MQ=255)
                                                                                                               cacGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTgcgc                        <  2:229122/121‑1 (MQ=255)
                                                                                                               cacGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTgcgc                        >  1:229122/1‑121 (MQ=255)
                                                                                                                   ccTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTa       <  1:744162/134‑1 (MQ=255)
                                                                                                                   ccTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTa       <  2:738078/134‑1 (MQ=255)
                                                                                                                      gAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACgg    >  1:511237/1‑134 (MQ=255)
                                                                                                                        aTACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGgcc  <  2:511237/134‑1 (MQ=255)
                                                                                                                                |                                                                                                                             
GACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCC  >  AM260479/582631‑582884

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCCTTGCCTTCAC  >  AM260479/582627‑582894
                                                                                                                                    |                                                                                                                                       
tccatccaagctcgttcctttttctgtgtttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCC                                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1103:12677:46728/44‑1 (MQ=60)
           GAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCC                                                                                                               >  GWNJ‑0478:712:GW2002102894th:2:2212:7970:54639/1‑148 (MQ=60)
              AATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAA                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1114:13770:92386/150‑1 (MQ=60)
              AATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAA                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1215:18004:98485/150‑1 (MQ=60)
              AATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAA                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1215:12031:59383/150‑1 (MQ=60)
               tcttatttctatcttttatttgtacttttttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAA                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1208:7226:73054/42‑1 (MQ=60)
                         GTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCtgtctcttatacacat                                                                                               >  GWNJ‑0478:712:GW2002102894th:2:1216:11051:77331/1‑134 (MQ=60)
                               TGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCT                                                                                         >  GWNJ‑0478:712:GW2002102894th:2:2204:18080:16524/1‑150 (MQ=60)
                                     GACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAA                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2211:19291:57667/150‑1 (MQ=60)
                                       CTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTctgtctcttatacacatctccgagcccacgagacg                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1107:15071:21360/1‑115 (MQ=60)
                                         CACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACC                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2104:11298:83755/148‑1 (MQ=60)
                                           CCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCG                                                                             >  GWNJ‑0478:712:GW2002102894th:2:1208:15584:54917/1‑150 (MQ=60)
                                                            agcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGC                                                            <  GWNJ‑0478:712:GW2002102894th:2:2105:1914:38533/86‑1 (MQ=60)
                                                             GATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCCGCGCTGCGCACGC                                                            >  GWNJ‑0478:712:GW2002102894th:2:2205:11782:30599/1‑149 (MQ=60)
                                                                                          ATGCCAGCAAGCGCAAGGTCGCCATCACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTC                              >  GWNJ‑0478:712:GW2002102894th:2:2115:2921:98335/1‑150 (MQ=60)
                                                                                                GCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGC                        <  GWNJ‑0478:712:GW2002102894th:2:1210:8648:88827/150‑1 (MQ=60)
                                                                                                    GCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGC                     >  GWNJ‑0478:712:GW2002102894th:2:1114:5239:43938/1‑149 (MQ=60)
                                                                                                           GTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGC             >  GWNJ‑0478:712:GW2002102894th:2:1216:20579:54509/1‑150 (MQ=60)
                                                                                                                      acagGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCCTTGCCTTCAC  <  GWNJ‑0478:712:GW2002102894th:2:1212:17433:58172/146‑1 (MQ=60)
                                                                                                                      GCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCCTTGCCTTCAC  <  GWNJ‑0478:712:GW2002102894th:2:2114:17562:32708/150‑1 (MQ=60)
                                                                                                                       CCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCCTctgtctctt  >  GWNJ‑0478:712:GW2002102894th:2:2106:11605:101231/1‑140 (MQ=60)
                                                                                                                                    |                                                                                                                                       
CGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCCTTGCCTTCAC  >  AM260479/582627‑582894

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: