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
A8 F1 I2 R1	97	53.3	2484976	94.5%	2348302	111.1

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,293,106	(G)_7→8	coding (1595/1659 nt)	nrfE →	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,293,099	1	.	G	100.0%	109.9 / NA	29	W530G (TGG→GGG)	nrfE	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (15/14); total (15/14)

CCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCT‑GGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCG  >  NC_000913/4292969‑4293223

ccGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCT‑gggggggg                                                                                                                      <  2:1098670/139‑2 (MQ=255)

  gTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCT‑ggggggggag                                                                                                                    >  2:458101/1‑136 (MQ=255)

    caacaaATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCT‑ggggggggagga                                                                                                                  <  1:579567/139‑6 (MQ=255)

          atgatgGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGAttgttg                                                                                                             <  1:375584/139‑1 (MQ=255)

                      tCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGgc                                                                                                 >  1:188296/1‑139 (MQ=255)

                      tCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGgc                                                                                                 >  2:395480/1‑139 (MQ=255)

                         aTTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCa                                                                                              <  2:357058/139‑1 (MQ=255)

                                                                         aCTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTtct                                              <  1:395480/139‑1 (MQ=255)

                                                                             ggCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACg                                          <  2:567174/139‑1 (MQ=255)

                                                                              gCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGc                                         <  1:458101/139‑1 (MQ=255)

                                                                                               ttCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTAcc                        >  2:747244/1‑139 (MQ=255)

                                                                                               ttCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCT‑gg                                                                                                                            >  1:327792/1‑38 (MQ=255)

                                                                                               ttCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCT‑gg                                                                                                                            <  2:327792/38‑1 (MQ=255)

                                                                                                     tgtatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGC‑GGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTa                                       >  2:963281/1‑117 (MQ=255)

                                                                                                     tgtatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGC‑GGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTa                                       <  1:963281/117‑1 (MQ=255)

                                                                                                      gtatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTActgc                                   >  1:1200735/1‑121 (MQ=255)

                                                                                                      gtatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTActgc                                   <  2:1200735/121‑1 (MQ=255)

                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTtgat                                                                                                           <  1:1001528/39‑1 (MQ=255)

                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTtgat                                                                                                           >  2:1001528/1‑39 (MQ=255)

                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCa                                                                                              >  2:357537/1‑52 (MQ=255)

                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCa                                                                                              <  1:357537/52‑1 (MQ=255)

                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTa                                       <  1:391049/107‑1 (MQ=255)

                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTa                                       >  2:391049/1‑107 (MQ=255)

                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTAtt                            <  2:626122/118‑1 (MQ=255)

                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTAtt                            >  1:626122/1‑118 (MQ=255)

                                                                                                                aGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCa       >  2:198650/1‑139 (MQ=255)

                                                                                                                    gtgtGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGtt                    >  1:94090/1‑122 (MQ=255)

                                                                                                                    gtgtGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGtt                    <  2:94090/122‑1 (MQ=255)

                                                                                                                     tgtgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCg  >  1:1074091/1‑139 (MQ=255)

                                                                                                                     tgtgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCg  >  1:250391/1‑139 (MQ=255)

                                                                                                                       tgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCt                                                 <  2:145297/90‑1 (MQ=255)

                                                                                                                       tgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCt                                                 >  1:145297/1‑90 (MQ=255)

                                                                                                                          gcTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGAt                                                                                 <  1:989312/55‑1 (MQ=255)

                                                                                                                          gcTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGAt                                                                                 >  2:989312/1‑55 (MQ=255)

CCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCT‑GGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCG  >  NC_000913/4292969‑4293223

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

GATK/CNVnator alignment

N/A