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
A5 F13 I0 R1 213 60.4 2554362 98.6% 2518600 147.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 4,158,471 A→C 36.1% S161R (AGC→CGC)  PP_3661 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,158,4710AC36.1% 20.6 / 3.2 26S161R (AGC→CGC) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base A (10/6);  new base C (6/3);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTG  >  NC_002947_CJ‑RC/4158325‑4158612
                                                                                                                                                   |                                                                                                                                             
tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                             <  3:44111/149‑1 (MQ=255)
  cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcgc                                                                                                                                           <  7:283168/149‑1 (MQ=255)
        gTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGGGGCGCGccccccc                                                                                                                                     >  3:74209/1‑148 (MQ=255)
         tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCCGccccccccc                                                                                                                                     >  1:229686/1‑147 (MQ=255)
                 ccccTGCTCATCGGCCTGGGGTGGGGCATCAAAGTGTTCACGCTCGCCCGCTATACCGCCTTGCTCACCGCCGCCGGCAGCGCCATGTGCGGCGCCGCCGCCGTGCGGGCATTCTAATCGGCGCGGCGCAGCGCCCCCCATAAAAgcgc                                                                                                                            <  8:239060/148‑1 (MQ=255)
                  cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCCGCCCCCCCCCAaaaaaaac                                                                                                                            >  3:252863/1‑144 (MQ=255)
                   ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCCGCGCCCCCCCAaaaaaaaac                                                                                                                           >  3:265702/1‑143 (MQ=255)
                   ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCACCCCCCCCCCAaaaaaaccc                                                                                                                           >  1:84118/1‑143 (MQ=255)
                   ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACAGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGCCGTGCCCGCCGCCCCCCAAAAACCCGcc                                                                                                                           >  4:123498/1‑149 (MQ=255)
                    ctgctgCGCGGCGTGGGGTGGGGCATGACTGTGTCCCCGCGCGGCCGCGAACCCGCCGTGCCCCCCGCCGCCGGCCGCGCCCGCGGCGGCGCCGCCGCGGTGCGGGCATCCGACGCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCAt                                                                                                                         <  8:199316/149‑1 (MQ=255)
                     tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAAAGCCaag                                                                                                                        >  1:203889/1‑147 (MQ=255)
                        tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCACCGCCCCCCAAAAAAACGCCATgg                                                                                                                       >  6:163248/1‑147 (MQ=255)
                              gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCAGCCCCCCCCAAAAAAACCCCATCGCAGTccg                                                                                                                >  6:249412/1‑146 (MQ=255)
                                                       cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGGGCAGCGCCCCCCCAAAAATCGCCAACTCAGTACGCAGCGGAGTGGGGTTGTGTAGCTTGt                                                                                      >  7:186258/1‑149 (MQ=255)
                                                       cAAGCGCGCCCGCTAACCCGCCTTGCTGCCCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCGGGGGGGGCATTCTAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGt                                                                                      <  8:271671/149‑1 (MQ=255)
                                                                ccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCCCCGCCCCCCCCAAAAAAGCCCAGGCCCGCGGGAGCGGAGTGGGGTTGGGTACCGTGTCggggttg                                                                              >  5:238680/1‑142 (MQ=255)
                                                                    gATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCAAAAAACGCCATGGCAGTCGGAAGCGGAGTGCGGTTCGGTACCATTTGGATGttgcgtt                                                                          >  8:61018/1‑143 (MQ=255)
                                                                                       cccccgGCAGCGCCTTCGGCGGCGCCGCCGCGGTGCTGGCATTCGAAGCGGCGGGCCTCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGcc                                                      <  8:37803/146‑1 (MQ=255)
                                                                                                   ccttCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCAt                                          <  2:84118/146‑1 (MQ=255)
                                                                                                         gcggcgCCGCCGCCGGGCTGGCATTCCACGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGAt                                    <  4:74209/149‑1 (MQ=255)
                                                                                                           ggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTGCGAAGCGTAGTGCTGTTAGGGACCATGTCGATGTTGCATTACCCACAGGCGATCAAATCCGGCGGGCTGCATCTGGGTaa                                  >  6:154716/1‑148 (MQ=255)
                                                                                                                            ggCATTCGAAGCGGCGCGGGGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgc                 <  2:229686/149‑1 (MQ=255)
                                                                                                                                          gcgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAc    >  2:181321/1‑149 (MQ=255)
                                                                                                                                          gcgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAc    >  1:135211/1‑149 (MQ=255)
                                                                                                                                            gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTg  >  2:38455/1‑149 (MQ=255)
                                                                                                                                            gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTg  >  2:33602/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTG  >  NC_002947_CJ‑RC/4158325‑4158612

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

GATK/CNVnator alignment

N/A