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
A24 F999 I0 R1 116 27.9 2350853 86.4% 2031136 49.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,720,597 Δ1 bp coding (4/153 nt) hokA ← small toxic polypeptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,720,5960G.100.0% 209.8 / NA 44coding (5/153 nt)hokAsmall toxic polypeptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (24/20);  total (24/20)

AAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720548‑3720640
                                                |                                            
aaaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                              >  1:2138191/1‑49 (MQ=255)
 aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                            <  1:1878249/50‑3 (MQ=255)
     cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:1826353/1‑50 (MQ=255)
      gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                      <  1:1163294/50‑1 (MQ=255)
      gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:1801534/1‑49 (MQ=39)
            cAAATCACTATTGAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                >  1:2318958/1‑50 (MQ=255)
             aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                               >  1:770746/1‑50 (MQ=255)
              aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                              <  1:588146/50‑1 (MQ=255)
              aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                              <  1:443608/50‑1 (MQ=255)
                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            <  1:478101/50‑1 (MQ=255)
                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            <  1:2098859/50‑1 (MQ=255)
                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                          >  1:16285/1‑50 (MQ=255)
                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                          >  1:1857451/1‑50 (MQ=255)
                    tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:942085/1‑49 (MQ=255)
                    tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:969733/1‑49 (MQ=255)
                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      >  1:152335/1‑50 (MQ=255)
                        aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                    >  1:226723/1‑50 (MQ=255)
                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                    <  1:1905346/49‑1 (MQ=255)
                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  1:1785424/1‑50 (MQ=255)
                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  <  1:631532/50‑1 (MQ=255)
                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  <  1:2195512/50‑1 (MQ=255)
                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  1:1456312/1‑49 (MQ=255)
                           gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                 <  1:526325/50‑1 (MQ=255)
                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  1:1638850/1‑50 (MQ=255)
                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  1:406786/1‑50 (MQ=255)
                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               <  1:1455371/50‑1 (MQ=255)
                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               <  1:953837/50‑1 (MQ=255)
                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg              <  1:2116012/50‑1 (MQ=255)
                               gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGa             <  1:128459/50‑1 (MQ=255)
                                tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            <  1:1925247/50‑1 (MQ=255)
                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt           <  1:489863/49‑1 (MQ=255)
                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:613110/1‑50 (MQ=255)
                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:1490281/1‑50 (MQ=255)
                                   tCTATATTTCTGC‑GCATTACCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:413539/1‑49 (MQ=255)
                                    cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  1:125017/1‑49 (MQ=255)
                                    cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  1:319408/1‑49 (MQ=255)
                                    cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  1:1443582/1‑49 (MQ=255)
                                      atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac      >  1:890584‑M1/1‑47 (MQ=255)
                                      atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac      >  1:939746‑M1/1‑47 (MQ=255)
                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       <  1:933878‑M1/48‑3 (MQ=255)
                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt    >  1:210203‑M1/1‑45 (MQ=255)
                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   <  1:585513‑M1/50‑7 (MQ=255)
                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   <  1:2223688‑M1/50‑7 (MQ=255)
                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  1:1491990‑M1/1‑44 (MQ=255)
                                           tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg  <  1:429848‑M1/49‑8 (MQ=255)
                                           tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg  <  1:787128‑M1/49‑8 (MQ=255)
                                                |                                            
AAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720548‑3720640

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

N/A