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
A4 F60 I1 R1 17 33.2 1123130 98.8% 1109652 138.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,153,675 G→T intergenic (‑577/+21) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,153,6750GT93.3% 88.0 / 1.6 30intergenic (‑577/+21)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (14/14);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

GGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAG  >  NC_000913/4153545‑4153805
                                                                                                                                  |                                                                                                                                  
ggcgCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATc                                                                                                                           >  1:72896/1‑140 (MQ=255)
                 cGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATg                                                                                                          <  1:349609/140‑1 (MQ=255)
                            tGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAAt                                                                                               <  2:420131/140‑1 (MQ=255)
                             gggCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATg                                                                                              >  2:197277/1‑140 (MQ=255)
                                   cGTTTGAGTCACTGTAGGTCGGATAGGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAt                                                                                        >  2:384350/1‑140 (MQ=255)
                                    gTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAtt                                                                                       >  2:409704/1‑140 (MQ=255)
                                     tttGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTa                                                                                      >  1:480636/1‑140 (MQ=255)
                                       tGAGTCACTGTAGGTCGGATAGGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                    <  1:384350/140‑1 (MQ=255)
                                       tGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                    <  2:480636/140‑1 (MQ=255)
                                             aCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGtt                                                                              >  2:51699/1‑140 (MQ=255)
                                               tcgatgTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTAt                                                                            <  1:395918/135‑1 (MQ=255)
                                               tGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTAt                                                                            >  2:236895/1‑140 (MQ=255)
                                                   ggTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                    <  1:457494/128‑1 (MQ=255)
                                                   ggTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                    >  2:457494/1‑128 (MQ=255)
                                                   ggTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGt                                                                        >  2:175388/1‑140 (MQ=255)
                                                    gTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGtt                                                                       >  1:319919/1‑140 (MQ=255)
                                                    gTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGtt                                                                       >  1:531946/1‑140 (MQ=255)
                                                    gTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGtt                                                                       <  2:319919/140‑1 (MQ=255)
                                                         aTAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCgg                                                                  >  1:446190/1‑140 (MQ=255)
                                                           aaGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGt                                                                <  1:236895/140‑1 (MQ=255)
                                                             gATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTgg                                                              <  1:409704/140‑1 (MQ=255)
                                                              aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTggg                                                             >  2:317939/1‑140 (MQ=255)
                                                                            catcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGGTTTAATGCCAGCAAAAATGGTGAATTACCTGGGGTATCAGTTCGCGGGAGGGCTTGATAAACCgtg                                               >  1:171789/1‑140 (MQ=255)
                                                                                 gACATTATTGCGGCACTGGAGTTTGGCATCAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTcc                                          <  1:503773/140‑1 (MQ=255)
                                                                                 gACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTcc                                          <  2:72896/140‑1 (MQ=255)
                                                                                          gCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtca                                 <  1:481049/140‑1 (MQ=255)
                                                                                            ggCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtcatc                               >  1:284802/1‑140 (MQ=255)
                                                                                              cACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAg                             <  2:284802/140‑1 (MQ=255)
                                                                                                  ggAGTTTGGCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTtg                         <  2:531946/140‑1 (MQ=255)
                                                                                                      tttCGCAACAGCGGCGGATGCGTGCCTATCGCGGTAGCCGGTCTACGTTTCAATGCCAGCAGAAATGGTGAATTACCTGGGTTCTCAGTTCGCGTGTGGGCTTGATAAACCGTGTTTCCAGACATTCATCAGGTtgatga                     <  2:171789/140‑1 (MQ=255)
                                                                                                         ggCAACAGCGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGcc                  >  2:220294/1‑140 (MQ=255)
                                                                                                                         gcggcgCGATCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAg  <  1:197277/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                  
GGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAG  >  NC_000913/4153545‑4153805

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

N/A