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 F1 I1 R1	770	51.4	2439994	93.8%	2288714	110.3

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	812,947	G→A	226 (TAG→TAA)	bioD →	dethiobiotin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	812,947	0	G	A	97.1%	101.4 / ‑6.2	34	226 (TAG→TAA)	bioD	dethiobiotin synthetase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (17/16); total (17/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

GGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAT  >  NC_000913/812812‑813072

gggAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGTGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCa                                                                                                                            <  2:1074130/139‑1 (MQ=255)

   aaaaCGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTc                                                                                                                         <  2:370151/139‑1 (MQ=255)

                           gaCGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGAACCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAATCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTAGGTTAAATTgctagcga                                                                                                 >  2:750027/2‑139 (MQ=255)

                           cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgagcga                                                                                                 >  2:717456/1‑139 (MQ=255)

                           cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgagcga                                                                                                 >  1:844596/1‑139 (MQ=255)

                             cGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCgaga                                                                                               >  2:869832/1‑139 (MQ=255)

                                              cccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              <  2:844596/139‑1 (MQ=255)

                                              cccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              <  2:220617/139‑1 (MQ=255)

                                              cccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              >  2:839377/1‑139 (MQ=255)

                                                               agagaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTa                                                             <  1:839377/139‑1 (MQ=255)

                                                                                                 gCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTc                           >  2:337299/1‑139 (MQ=255)

                                                                                                  cGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgc                                                                                                       <  2:573092/62‑1 (MQ=255)

                                                                                                  cGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgc                                                                                                       >  1:573092/1‑62 (MQ=255)

                                                                                                        ccGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCa                    <  1:717456/139‑1 (MQ=255)

                                                                                                          ggAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             >  1:604991/1‑64 (MQ=255)

                                                                                                          ggAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             >  1:358041/1‑64 (MQ=255)

                                                                                                          ggAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             <  2:604991/64‑1 (MQ=255)

                                                                                                          ggAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             <  2:358041/64‑1 (MQ=255)

                                                                                                           gAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           >  1:370162/1‑97 (MQ=255)

                                                                                                           gAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           <  2:370162/97‑1 (MQ=255)

                                                                                                               gTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg             <  1:337299/139‑1 (MQ=255)

                                                                                                               gTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg             <  2:937585/139‑1 (MQ=255)

                                                                                                                   aTAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAAtg         >  2:618638/1‑139 (MQ=255)

                                                                                                                       aCCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGc                        <  1:237350/120‑1 (MQ=255)

                                                                                                                       aCCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGc                        >  2:237350/1‑120 (MQ=255)

                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           >  1:865309/1‑82 (MQ=255)

                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           <  2:865309/82‑1 (MQ=255)

                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt  >  2:275175/1‑139 (MQ=255)

                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt  >  1:262657/1‑139 (MQ=255)

                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt  <  1:243843/139‑1 (MQ=255)

                                                                                                                            gccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATcgc                                                                                           <  2:1094479/48‑1 (MQ=255)

                                                                                                                            gccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATcgc                                                                                           >  1:1094479/1‑48 (MQ=255)

                                                                                                                               ttgttgTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                                 >  1:966359/1‑55 (MQ=255)

                                                                                                                               ttgttgTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                                 <  2:966359/55‑1 (MQ=255)

GGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAT  >  NC_000913/812812‑813072

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

GATK/CNVnator alignment

N/A