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 I2 R1 524 31.4 1748232 97.8% 1709770 80.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,474,329 T→G Q75P (CAG→CCG)  rpsL ← 30S ribosomal subunit protein S12

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,474,3290TG97.4% 130.2 / ‑4.0 38Q75P (CAG→CCG) rpsL30S ribosomal subunit protein S12
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (18/19);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACG  >  NC_000913/3474253‑3474405
                                                                            |                                                                            
accacGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAg                                                                          <  2:614017/81‑1 (MQ=255)
            gtgGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCtt                                                               >  2:614064/1‑80 (MQ=255)
             tgGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTc                                                              <  2:614018/80‑1 (MQ=255)
              ggTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTcac                                                            >  2:614066/1‑81 (MQ=255)
               gTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTcac                                                            <  2:614019/80‑1 (MQ=255)
                  aCGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCaccacc                                                        <  2:614020/81‑1 (MQ=255)
                   cGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCaccacc                                                        >  2:614061/1‑80 (MQ=255)
                    gAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGa                                                      >  1:614051/1‑81 (MQ=255)
                    gAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGa                                                      >  1:614077/1‑81 (MQ=255)
                     aaCACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGAt                                                     >  1:614093/1‑81 (MQ=255)
                       cacCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGt                                                   <  2:614021/81‑1 (MQ=255)
                       cacCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATg                                                    >  1:614074/1‑80 (MQ=255)
                             ggAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAg                                             <  1:614022/81‑1 (MQ=255)
                              gAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGt                                            <  1:614023/81‑1 (MQ=255)
                              gAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGt                                            >  1:614069/1‑81 (MQ=255)
                               aGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTc                                           <  2:614024/81‑1 (MQ=255)
                                 gTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCAc                                         >  1:614079/1‑81 (MQ=255)
                                  tCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACt                                        >  2:614094/1‑81 (MQ=255)
                                     ttAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGT‑‑TGACCTTCACCACCGATGTAGGAAGTCACTTCGa                                    >  1:614072/1‑80 (MQ=255)
                                      tAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCg                                     >  1:614075/1‑80 (MQ=255)
                                      tAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCg                                     >  1:614073/1‑80 (MQ=255)
                                          acGACCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAAcc                                <  2:614025/81‑1 (MQ=255)
                                             aCCGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGtt                             <  1:614026/81‑1 (MQ=255)
                                               cGCCACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAg                           >  1:614080/1‑81 (MQ=255)
                                                 ccACGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTc                         <  2:614027/81‑1 (MQ=255)
                                                   aCGGATCAGGATCACGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAg                       <  1:614028/81‑1 (MQ=255)
                                                   aCGGATCAGGATCACGGAGTGCTCCGGCAGGATGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAg                       <  2:614029/81‑1 (MQ=255)
                                                     ggatcaggatcaCGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGAc                     >  1:614092/1‑81 (MQ=255)
                                                     ggatcaggatcaCGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGAc                     <  1:614030/81‑1 (MQ=255)
                                                          aggatcaCGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAaca                <  2:614031/81‑1 (MQ=255)
                                                                aCGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCa           <  2:614033/80‑1 (MQ=255)
                                                                aCGGAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCAt          <  2:614032/81‑1 (MQ=255)
                                                                  ggAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATAc        >  1:614081/1‑81 (MQ=255)
                                                                   gAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACt       <  2:614036/81‑1 (MQ=255)
                                                                   gAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACt       <  2:614034/81‑1 (MQ=255)
                                                                   gAGTGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACt       <  1:614035/81‑1 (MQ=255)
                                                                      tGCTCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTa    >  2:614097/1‑81 (MQ=255)
                                                                         tCCGGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACg  >  2:614102/1‑80 (MQ=255)
                                                                            |                                                                            
ACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACG  >  NC_000913/3474253‑3474405

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

GATK/CNVnator alignment

N/A