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
A6 F1 I1 R1	750	51.0	2595496	93.6%	2429384	106.4

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,389	A→G	K163R (AAG→AGG)	cybB →	cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,389	0	A	G	93.5%	88.7 / ‑2.2	31	K163R (AAG→AGG)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (12/17); total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

CGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATT  >  NC_000913/1491254‑1491514

cGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                            <  1:642752/139‑1 (MQ=255)

                               aCGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTaaa                                                                                             <  2:1023721/139‑1 (MQ=255)

                                           cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGTACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                            <  2:323871/96‑1 (MQ=255)

                                           cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGTACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                            >  1:323871/1‑96 (MQ=255)

                                           cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGt                                                                                          <  2:7820/130‑1 (MQ=255)

                                           cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGt                                                                                          >  1:7820/1‑130 (MQ=255)

                                           cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGAGATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAgg                                                                                 >  2:849320/1‑139 (MQ=255)

                                             tAAAGTCGTTGCATGAAACGCTGGCGAATCTGGGATCTTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGCCGACACACTTATACGCATGATGCCGCGTAAACGTCCCTGAAGGat                                                                               >  2:720898/1‑139 (MQ=255)

                                             tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACACTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGat                                                                               >  2:857427/1‑139 (MQ=255)

                                                aGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATAtt                                                                            >  1:550277/1‑139 (MQ=255)

                                                     tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCATACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAg                                                                       >  2:40741/1‑139 (MQ=255)

                                                     tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAg                                                                       >  2:58257/1‑139 (MQ=255)

                                                     tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAg                                                                       >  2:838449/1‑139 (MQ=255)

                                                     tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAg                                                                       >  2:100750/1‑139 (MQ=255)

                                                         aTGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaaaa                                                                   <  1:849320/139‑1 (MQ=255)

                                                                          cTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCg                                                  <  1:100750/139‑1 (MQ=255)

                                                                                tatTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGt                                            <  2:101414/139‑1 (MQ=255)

                                                                                    ttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCATACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                        <  1:40741/139‑1 (MQ=255)

                                                                                    ttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                        <  2:1050956/139‑1 (MQ=255)

                                                                                       tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacaca                                                                                                                     >  2:1080161/1‑59 (MQ=255)

                                                                                       tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacaca                                                                                                                     <  1:1080161/59‑1 (MQ=255)

                                                                                       tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAg                                     <  1:58257/139‑1 (MQ=255)

                                                                                                       tGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATAc                          <  2:574656/134‑1 (MQ=255)

                                                                                                       tGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATAc                          >  1:574656/1‑134 (MQ=255)

                                                                                                         cGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                        >  2:208231/1‑118 (MQ=255)

                                                                                                         cGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                        <  1:208231/118‑1 (MQ=255)

                                                                                                              gCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCg              <  1:857427/139‑1 (MQ=255)

                                                                                                                aCTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGa            <  1:113760/139‑1 (MQ=255)

                                                                                                                   ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgc         <  2:776273/139‑1 (MQ=255)

                                                                                                                      acacCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGcc      <  2:550277/139‑1 (MQ=255)

                                                                                                                          cacGATTTCTGGAGGGACAACACACCTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGACTACGGGTCTAAATCGGAAGCGCCtatt  <  1:836728/139‑1 (MQ=255)

CGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATT  >  NC_000913/1491254‑1491514

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

GATK/CNVnator alignment

N/A