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 150 26.0 2053212 90.0% 1847890 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,071,935 +G coding (626/627 nt) h16_A0979 → Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,071,9341.G100.0% 96.1 / NA 26*209G (TGA→GGA) h16_A0979Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (16/10);  total (16/10)

CGCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGG  >  AM260479/1071809‑1072031
                                                                                                                              |                                                                                                  
cgcgCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGccc                                                                                             >  2:1026489/1‑134 (MQ=255)
    cAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAgg                                                                                         >  1:89194/1‑134 (MQ=255)
                    gAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCt                                                                         <  1:333605/134‑1 (MQ=255)
                     aGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTc                                                                        >  2:50962/1‑134 (MQ=255)
                         ggggTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCg                                                                    >  1:431576/1‑134 (MQ=255)
                               gTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATc                                                               >  1:986358/1‑133 (MQ=255)
                                   gTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAgg                                                           >  2:796134/1‑133 (MQ=255)
                                    tCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAgg                                                           >  2:918298/1‑132 (MQ=255)
                                      gCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGc                                                       >  2:1006457/1‑134 (MQ=255)
                                          gCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTg                                                   <  2:986358/134‑1 (MQ=255)
                                                     gCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATc                                                               <  1:560868/111‑1 (MQ=255)
                                                     gCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATc                                                               >  2:560868/1‑111 (MQ=255)
                                                      cccTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAg                                      <  2:431576/134‑1 (MQ=255)
                                                               cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGa                                                                                           <  1:241654/73‑1 (MQ=255)
                                                               cgcTGTGGCCGGGCTGGACAGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGa                                                                                           >  2:241654/1‑73 (MQ=39)
                                                                    tgGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGa                         >  2:936360/1‑134 (MQ=255)
                                                                         gggCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcg                    >  2:886808/1‑134 (MQ=255)
                                                                           gCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcgcg                  <  1:886808/134‑1 (MQ=255)
                                                                                  cGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAg                                                            <  1:723374/85‑1 (MQ=255)
                                                                                  cGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAg                                                            >  2:723374/1‑85 (MQ=255)
                                                                                        cGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGc     >  1:761723/1‑134 (MQ=255)
                                                                                           cATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTgg  <  1:918298/134‑1 (MQ=255)
                                                                                                                gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGc                                                                              >  2:417144/1‑37 (MQ=255)
                                                                                                                gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGc                                                                              <  1:417144/37‑1 (MQ=255)
                                                                                                                gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcg                    >  2:208341/1‑95 (MQ=255)
                                                                                                                gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcg                    <  1:208341/95‑1 (MQ=255)
                                                                                                                              |                                                                                                  
CGCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGG  >  AM260479/1071809‑1072031

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
CGACCTGGCCACGCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGC  >  AM260479/1071798‑1072069
                                                                                                                                          |                                                                                                                                       
CGACCTGGCCACGCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGG                                                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1107:20465:100813/1‑150 (MQ=60)
  ttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGG                                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1213:14562:25787/78‑1 (MQ=60)
            gatatgtagctccgcctcgtgggctcggagatgtgtataagagacagaccCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCAC                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:1116:1404:24666/100‑1 (MQ=60)
                     CTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGT                                                                                                         >  GWNJ‑0478:712:GW2002102894th:2:1213:5189:33983/1‑150 (MQ=60)
                         gttttttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATC                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1210:2046:51881/67‑1 (MQ=60)
                           GGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCA                                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:2214:15338:68956/1‑149 (MQ=60)
                              CGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGA                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1207:2449:45681/150‑1 (MQ=60)
                                        aagagacagGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCC                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2105:4864:31392/141‑1 (MQ=60)
                                           cgtctcgtgggctcggagatgtgtataagagacagGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACC                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2115:5074:44406/115‑1 (MQ=60)
                                                 GCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGgacctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtct                                                                            >  GWNJ‑0478:712:GW2002102894th:2:1211:17185:92575/1‑89 (MQ=60)
                                                                                   gctcggagatgtgtataagagacagGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGC                                           <  GWNJ‑0478:712:GW2002102894th:2:1115:3319:28640/125‑1 (MQ=60)
                                                                                    GGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCT                                          >  GWNJ‑0478:712:GW2002102894th:2:2116:8459:70513/1‑150 (MQ=60)
                                                                                     GGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTG                                         <  GWNJ‑0478:712:GW2002102894th:2:2207:17241:84208/150‑1 (MQ=60)
                                                                                                     GCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGatgcggcccaggtcc                         <  GWNJ‑0478:712:GW2002102894th:2:2211:20466:76637/150‑16 (MQ=60)
                                                                                                                        GCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGActgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaaagaaac      >  GWNJ‑0478:712:GW2002102894th:2:2215:17822:65862/1‑60 (MQ=60)
                                                                                                                              GACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGC  <  GWNJ‑0478:712:GW2002102894th:2:2203:1450:14546/148‑1 (MQ=60)
                                                                                                                                          |                                                                                                                                       
CGACCTGGCCACGCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGC  >  AM260479/1071798‑1072069

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: