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 2,142,743 G→A P69L (CCC→CTC)  udk ← uridine/cytidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,142,7430GA97.6% 137.8 / NA 41P69L (CCC→CTC) udkuridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/17);  minor base T (1/0);  total (24/17)
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

ACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAAT  >  NC_000913/2142668‑2142819
                                                                           |                                                                            
aCCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGTTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATgg                                                                         >  2:363929/1‑81 (MQ=255)
  cGGCAGGTCAATTGCCGAGCCGCGTTTCAACGTTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTc                                                                       <  2:363892/81‑1 (MQ=255)
  cGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATTCATCGCGCTGAGATGGTc                                                                       <  2:363894/81‑1 (MQ=255)
  cGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTc                                                                       >  2:363928/1‑81 (MQ=255)
  cGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTc                                                                       <  1:363893/81‑1 (MQ=255)
   ggCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCg                                                                      <  2:363895/81‑1 (MQ=255)
     cAGGTCAATTGCCGAGCCGCGTTTCAACGTTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCg                                                                      >  1:363930/1‑79 (MQ=255)
     cAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCg                                                                      >  2:363943/1‑79 (MQ=255)
     cAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCg                                                                      >  1:363950/1‑79 (MQ=255)
      aGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAg                                                                   >  1:363939/1‑81 (MQ=255)
           aaTTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGt                                                              <  1:363896/81‑1 (MQ=255)
            aTTGCCGAGCCGCGTTTCAACGTTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTc                                                             <  1:363897/81‑1 (MQ=255)
              tGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCt                                                            >  1:363937/1‑80 (MQ=255)
               gCCGAGCCGCGTTTCAACGTTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCtt                                                           <  1:363898/80‑1 (MQ=255)
                  gAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAAcg                                                       >  2:363934/1‑81 (MQ=255)
                     ccGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGt                                                    >  1:363952/1‑81 (MQ=255)
                       gcgTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGttc                                                  >  2:363941/1‑81 (MQ=255)
                       gcgTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGttc                                                  >  1:363942/1‑81 (MQ=255)
                       gcgTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGttc                                                  <  2:363899/81‑1 (MQ=255)
                        cgTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGttct                                                 >  1:363932/1‑81 (MQ=255)
                         gTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGttctt                                                >  1:363983/1‑81 (MQ=255)
                          tttCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCTCTGTGATGGTCGTAGTTGGTCTTAACGCGttcttc                                               >  1:363931/1‑81 (MQ=255)
                            tCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTcc                                              <  1:363900/80‑1 (MQ=255)
                             cAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCa                                             <  2:363901/80‑1 (MQ=255)
                             cAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCAt                                            >  2:363933/1‑81 (MQ=255)
                              aaCGTTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATc                                           >  2:363936/1‑81 (MQ=255)
                               aCGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCg                                          <  2:363902/81‑1 (MQ=255)
                                 gCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGAc                                        <  1:363903/81‑1 (MQ=255)
                                     gTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAg                                      <  2:363904/79‑1 (MQ=255)
                                      tAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGAt                                    <  1:363905/80‑1 (MQ=255)
                                        aaTGCTCAAGCAGCAGACTGTGAACCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATgg                                  <  1:363906/80‑1 (MQ=255)
                                            cTCAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCttt                              >  1:363951/1‑80 (MQ=255)
                                              cAAGCAGCAGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGAt                           >  2:363956/1‑81 (MQ=255)
                                                 gcagcaGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGGTGGCTTTGATCtt                        <  1:363907/81‑1 (MQ=255)
                                                 gcagcaGACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCtt                        >  2:363935/1‑81 (MQ=255)
                                                     cagACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGt                     >  2:363947/1‑80 (MQ=255)
                                                       gACTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTaa                   >  1:363960/1‑80 (MQ=255)
                                                        aCTGTGATCCATCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTaa                   >  2:363958/1‑79 (MQ=255)
                                                                  aTCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCg       >  1:363945/1‑81 (MQ=255)
                                                                   tCGCGCTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCgg      <  1:363908/81‑1 (MQ=255)
                                                                       gcTGAGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAAt  >  2:363969/1‑81 (MQ=255)
                                                                           |                                                                            
ACCGGCAGGTCAATTGCCGAGCCGCGTTTCAACGCTTGTAAATGCTCAAGCAGCAGACTGTGATCCATCGCGCTGGGATGGTCGTAGTTGGTCTTAACGCGTTCTTCCATCGACAGATGGCTTTGATCTTTGTAATAGCAGTCTTCGGGAAT  >  NC_000913/2142668‑2142819

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

GATK/CNVnator alignment

N/A