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 F3 I160 R1 220 104.1 6051714 95.9% 5803593 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,807,028 T→C E121G (GAA→GGA)  BW25113_RS09050 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,807,0280TC100.0% 134.1 / NA 36E121G (GAA→GGA) BW25113_RS09050hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/17);  total (19/17)

GATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGATTATTA  >  NZ_CP009273/1806942‑1807114
                                                                                      |                                                                                      
gATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTg                                                                                     <  1:395896/90‑1 (MQ=255)
gATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTg                                                                                     <  1:704/90‑1 (MQ=255)
     atataATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                                <  2:2734530/90‑1 (MQ=255)
             ttGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGa                                                                        >  2:2941447/1‑90 (MQ=255)
               gTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAg                                                                      >  2:2026239/1‑90 (MQ=255)
               gTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAg                                                                      >  2:1116854/1‑90 (MQ=255)
                 tCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    <  2:2509588/90‑1 (MQ=255)
                 tCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    <  1:992380/90‑1 (MQ=255)
                 tCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    <  1:926837/90‑1 (MQ=255)
                  cATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGa                                                                   <  1:2928023/90‑1 (MQ=255)
                    tAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACa                                                                 >  2:2735565/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  2:2055300/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  2:2566603/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  2:2862925/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  1:1752829/1‑90 (MQ=255)
                                cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa                                                     >  2:3002263/1‑90 (MQ=255)
                                       aGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              <  2:1241337/90‑1 (MQ=255)
                                       aGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              <  2:861470/90‑1 (MQ=255)
                                        gTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              >  1:2448512/1‑89 (MQ=255)
                                        gTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              <  2:2448512/89‑1 (MQ=255)
                                          gAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGaaa                                                                              >  1:1790050/1‑55 (MQ=255)
                                          gAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGaaa                                                                              <  2:1790050/55‑1 (MQ=255)
                                               accaaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                                >  1:1199149/1‑48 (MQ=255)
                                               accaaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa                                                                                <  2:1199149/48‑1 (MQ=255)
                                                  aaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              <  2:2586994/79‑1 (MQ=255)
                                                  aaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCtt                                              >  1:2586994/1‑79 (MQ=255)
                                                  aaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACct                                   >  2:1282950/1‑90 (MQ=255)
                                                           tttATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGa                          <  2:1430890/90‑1 (MQ=255)
                                                                 cTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    <  2:2588371/42‑1 (MQ=255)
                                                                 cTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTg                                                                    >  1:2588371/1‑42 (MQ=255)
                                                                 cTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCa                    >  2:1300858/1‑90 (MQ=255)
                                                                  tGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCAc                   >  1:1602610/1‑90 (MQ=255)
                                                                   gATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCAt                                                            >  1:2920675/1‑48 (MQ=255)
                                                                   gATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCAt                                                            <  2:2920675/48‑1 (MQ=255)
                                                                   gATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACt                  >  2:2480338/1‑90 (MQ=255)
                                                                           tGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATaaa          <  1:1033961/90‑1 (MQ=255)
                                                                                   tttcCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGattatta  <  2:2388449/90‑1 (MQ=255)
                                                                                      |                                                                                      
GATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGATTATTA  >  NZ_CP009273/1806942‑1807114

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGCTCCCGACAGGATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGATTATTAAT  >  NZ_CP009273/1806930‑1807116
                                                                                                  |                                                                                        
TGCTCCTGACAGGATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCC                                                                                         <  SRR3722030.2001086/100‑1 (MQ=60)
  CTCCCGACAGGATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTG                                                                                       <  SRR3722030.333323/100‑1 (MQ=60)
            GATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATAC                                                                             <  SRR3722030.401836/100‑1 (MQ=60)
            GATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATAC                                                                             <  SRR3722030.716/100‑1 (MQ=60)
                             TCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTC                                                            <  SRR3722030.1007427/100‑1 (MQ=60)
                             TCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTC                                                            <  SRR3722030.941037/100‑1 (MQ=60)
                              CATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCC                                                           <  SRR3722030.2982002/100‑1 (MQ=60)
                                  GTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATA                                                       >  SRR3722030.1781885/1‑100 (MQ=60)
                                          TCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTT                                               >  SRR3722030.2493215/1‑100 (MQ=60)
                                            CTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGctgtctcttatacacatctgacgct                                             >  SRR3722030.1820092/1‑75 (MQ=60)
                                                 ATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGctgtctcttatacacatctgacgctgccgacg                                        >  SRR3722030.1216891/1‑68 (MQ=60)
                                                    GTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCC                                     >  SRR3722030.2634464/1‑100 (MQ=60)
                                                                   GTTTTTTATACTGATCAAATTGGTAATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCtgtctcttatacacatctgacgctgccgacgacaagg                      >  SRR3722030.2635870/1‑63 (MQ=60)
                                                                    TTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCAC                     >  SRR3722030.1628592/1‑100 (MQ=60)
                                                                     TTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAActgtctcttatacacatctgacgctgccgacg                    >  SRR3722030.2974532/1‑68 (MQ=60)
                                                                                       TGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGATTATTAAT  <  SRR3722030.1049551/100‑1 (MQ=60)
                                                                                                  |                                                                                        
TGCTCCCGACAGGATGTATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGATTATTAAT  >  NZ_CP009273/1806930‑1807116

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: