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
A9 F48 I1 R1	82	50.4	1709896	99.0%	1692797	139.8

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	793,635	C→A	E52* (GAA→TAA)	modF ←	ABC family protein ModF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	793,635	0	C	A	90.0%	82.9 / ‑6.6	30	E52* (GAA→TAA)	modF	ABC family protein ModF
Reads supporting (aligned to +/- strand): ref base C (0/3); new base A (17/10); total (17/13)
Fisher's exact test for biased strand distribution p-value = 7.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCAA  >  NC_000913/793498‑793770

gtCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTacc                                                                                                                                       >  2:498221/1‑140 (MQ=255)

           cAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAgcgcg                                                                                                                            <  1:3772/140‑1 (MQ=255)

           cAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAgcgcg                                                                                                                            <  1:415713/140‑1 (MQ=255)

           aaGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAgcgcg                                                                                                                            <  2:93928/139‑1 (MQ=255)

             gggCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGcgcgcg                                                                                                                          >  1:74729/1‑140 (MQ=255)

                     cATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTNCCCCGCCAGCGCGCGGGCCAggg                                                                                                                  <  1:498221/140‑1 (MQ=255)

                          cGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGAt                                                                                                             >  1:849474/1‑140 (MQ=255)

                              gtTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTcc                                                                                                         >  2:136317/1‑140 (MQ=255)

                               ttATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTccc                                                                                                        >  2:357327/1‑140 (MQ=255)

                               ttATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTccc                                                                                                        >  1:254471/1‑140 (MQ=255)

                                  ttCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCt                                                                                                     >  2:562392/1‑140 (MQ=255)

                                  ttCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCt                                                                                                     >  2:827557/1‑140 (MQ=255)

                                       cTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCAt                                                                                                >  1:318540/1‑140 (MQ=255)

                                        agCCATTCGTCGCTGACGAGTTTTTGCAATTGCACGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCAtt                                                                                               <  2:168621/139‑1 (MQ=255)

                                           cATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGa                                                                                            >  2:167106/1‑140 (MQ=255)

                                            aTTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGaa                                                                                           <  2:820208/140‑1 (MQ=255)

                                                  cgCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACa                                                                                     <  2:254471/140‑1 (MQ=255)

                                                    cTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACaaa                                                                                   >  1:507781/1‑140 (MQ=255)

                                                        cGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGcc                                                                               <  1:167106/140‑1 (MQ=255)

                                                                     ttGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCcg                                                                  >  1:333650/1‑140 (MQ=255)

                                                                         acGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGtt                                                              <  2:647698/139‑1 (MQ=255)

                                                                          cGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGttt                                                             <  2:108951/140‑1 (MQ=255)

                                                                                    aCGAGTGATGTGGAAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCt                                                   <  2:724461/140‑1 (MQ=255)

                                                                                       aGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGaa                                                >  2:411575/1‑139 (MQ=255)

                                                                                       aGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAg                                                >  1:762392/1‑140 (MQ=255)

                                                                                       aGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAg                                                >  1:198140/1‑140 (MQ=255)

                                                                                                              gCCGTTCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTct                         >  2:336057/1‑139 (MQ=255)

                                                                                                                   tCACCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAAGCAATTGCAGCGTTTTTGTGTCGCTAAg                    <  1:411575/140‑1 (MQ=255)

                                                                                                                    aaCCTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGa                   <  2:29768/139‑1 (MQ=255)

                                                                                                                       cTTTCAAAAGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGa                >  1:469009/1‑140 (MQ=255)

                                                                                                                            aaaaGCGGAAGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGt           >  1:612545/1‑140 (MQ=255)

                                                                                                                                     aGTTACCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCaa  <  1:195038/140‑1 (MQ=255)

GTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCAA  >  NC_000913/793498‑793770

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

GATK/CNVnator alignment

N/A