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
A14 F1 I2 R1 774 49.2 2112398 96.9% 2046913 118.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,210,944 G→A pseudogene (1845/2001 nt) yehQ → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,210,9440GA97.9% 150.4 / ‑6.3 48pseudogene (1845/2001 nt)yehQpseudogene
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/17);  minor base C (0/1);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTC  >  NC_000913/2210809‑2211076
                                                                                                                                       |                                                                                                                                    
aCCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTg                                                                                                                                   >  2:957382/1‑139 (MQ=255)
              ccAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCt                                                                                                                     >  1:585791/1‑139 (MQ=255)
                                 gTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTAt                                                                                                  <  2:103874/139‑1 (MQ=255)
                                           catcaGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcac                                                                                        <  2:7602/139‑1 (MQ=255)
                                                 ttGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCGGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATgtgt                                                                                  >  2:639968/1‑139 (MQ=255)
                                                         aCGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                          >  2:648332/1‑139 (MQ=255)
                                                                  ggCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAAt                                                                 >  2:713878/1‑139 (MQ=255)
                                                                                 aCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCa                                                  >  2:296179/1‑139 (MQ=255)
                                                                                  cTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcac                                                                                        >  2:756895/1‑100 (MQ=255)
                                                                                  cTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcac                                                                                        <  1:756895/100‑1 (MQ=255)
                                                                                  cTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCAc                                                 >  1:101850/1‑139 (MQ=255)
                                                                                    gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACgg                                               >  2:1043610/1‑139 (MQ=255)
                                                                                      tCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGt                                             <  2:585791/139‑1 (MQ=255)
                                                                                             cGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCa                                      <  2:312778/139‑1 (MQ=255)
                                                                                               gCTGGTGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAAt                                    >  2:43991/1‑139 (MQ=255)
                                                                                               gCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGCGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAAt                                    >  2:230193/1‑139 (MQ=255)
                                                                                                cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGc                                                                                              <  1:987067/80‑1 (MQ=255)
                                                                                                cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGc                                                                                              >  2:987067/1‑80 (MQ=255)
                                                                                                cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATg                                   >  1:578519/1‑139 (MQ=255)
                                                                                                  ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGaa                                 >  1:533427/1‑139 (MQ=255)
                                                                                                  ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGaa                                 >  1:513177/1‑139 (MQ=255)
                                                                                                  ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGaa                                 >  2:741122/1‑139 (MQ=255)
                                                                                                         tCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTgg                                                                                           >  1:976733/1‑74 (MQ=255)
                                                                                                         tCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTgg                                                                                           <  2:976733/74‑1 (MQ=255)
                                                                                                         tCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGgt                          <  1:957382/139‑1 (MQ=255)
                                                                                                                   agTGCTATTGGTCAGGCATAATTGCTGGCGAATGTCCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTgg                                                                               <  2:651250/76‑1 (MQ=255)
                                                                                                                   agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTgg                                                                               >  1:651250/1‑76 (MQ=255)
                                                                                                                   agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                >  2:572614/1‑139 (MQ=255)
                                                                                                                   agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                >  2:1048369/1‑139 (MQ=255)
                                                                                                                   agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                >  1:358167/1‑139 (MQ=255)
                                                                                                                   agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                >  1:416375/1‑139 (MQ=255)
                                                                                                                   agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                >  1:535278/1‑139 (MQ=255)
                                                                                                                     tGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                       <  2:739512/66‑1 (MQ=255)
                                                                                                                     tGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                       >  1:739512/1‑66 (MQ=255)
                                                                                                                     tGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgct                  <  1:1018627/135‑1 (MQ=255)
                                                                                                                     tGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgct                  >  2:1018627/1‑135 (MQ=255)
                                                                                                                      gCTATTGGTCAGGCATACTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgctttgcc             <  1:648332/139‑2 (MQ=255)
                                                                                                                      gCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCg             <  2:101850/139‑1 (MQ=255)
                                                                                                                      gCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCg             <  2:535278/139‑1 (MQ=255)
                                                                                                                      gCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCg             <  1:713878/139‑1 (MQ=255)
                                                                                                                      gCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCg             <  1:1043610/139‑1 (MQ=255)
                                                                                                                       cTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGa            >  2:245115/1‑139 (MQ=255)
                                                                                                                          ttGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACa         <  2:533427/139‑1 (MQ=255)
                                                                                                                             gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCGCATGTGTTGGGACaa                                                                          <  2:166336/71‑1 (MQ=255)
                                                                                                                             gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCGCATGTGTTGGGACaa                                                                          >  1:166336/1‑71 (MQ=255)
                                                                                                                             gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAatt      >  2:146345/1‑139 (MQ=255)
                                                                                                                                 ggCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTc  >  1:782053/1‑139 (MQ=255)
                                                                                                                                 ggCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTc  >  2:331216/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                    
ACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTC  >  NC_000913/2210809‑2211076

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

GATK/CNVnator alignment

N/A