BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,959,480	A→T	5.5%	intergenic (+35/+52)	ilvC → / ← ppiC	ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,959,480	0	A	T	5.5%	756.5 / 24.5	219	intergenic (+35/+52)	ilvC/ppiC	ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand): ref base A (92/115); new base T (7/5); total (99/120)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

Only 100 of 234 total aligned reads displayed.

CTGAGTACGGTAACTATCTGTTCTCTTACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAG  >  NC_000913/3959219‑3959750

               aTCTGTTCTCTTACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCCGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGGGGATTAAGTTCGCGCTGGTGCCCCCCACCCGACCCCCCCCCACACGGACAGTGAGAAAACCCTGACGGcc                                                                                                                                                                                                                                                    >  4:337427/1‑275 (MQ=255)

                           aCGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                                                                                                                                                                                                                        <  4:87547/275‑1 (MQ=255)

                                                aaCCGGTTAGGGCACGGCTGGAACCGGGCGAACTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGGGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAACCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                                                                                                                                    <  4:215561/274‑1 (MQ=255)

                                                        aTGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                                                                                                                                                                                            <  4:270263/274‑1 (MQ=255)

                                                                          ggCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGTGCAACTGCGTTATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                                                                                                                                                                          <  1:1657/274‑1 (MQ=255)

                                                                                                 ggAAGGCGCGGTAGATAACGGGAAACTGCGTGAGGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGcc                                                                                                                                                                                                                                                    <  2:160665/193‑1 (MQ=255)

                                                                                                                 aaCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCca                                                                                                                                                                                                                                                                                >  8:625012/1‑149 (MQ=255)

                                                                                                                    gggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                                                                                                                                                              <  7:839675/148‑1 (MQ=255)

                                                                                                                       cAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                                                                                                                                                          <  6:161632/149‑1 (MQ=255)

                                                                                                                        aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggag                                                                                                                                                                                                                                                                         >  6:177439/1‑149 (MQ=255)

                                                                                                                        aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGACCTGGGTGTGAAGCGGGCCGGTGGGCTCCAGTACCGGCCAagag                                                                                                                                           >  2:54864/1‑275 (MQ=255)

                                                                                                                            gCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagaggg                                                                                                                                                                                                                                                                     <  6:1042789/149‑1 (MQ=255)

                                                                                                                              gTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagggag                                                                                                                                                                                                                                                                   >  7:1076069/1‑149 (MQ=255)

                                                                                                                                gATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaa                                                                                                                                                                                                                                                                 >  6:759842/1‑149 (MQ=255)

                                                                                                                                         gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                                                                                                                                                        <  5:1286078/149‑1 (MQ=255)

                                                                                                                                            gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                                                                                                                                                      <  7:1012858/132‑1 (MQ=255)

                                                                                                                                            gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                                                                                                                                                     <  8:263594/149‑1 (MQ=255)

                                                                                                                                               aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagggtc                                                                                                                                                                                                                                                                   <  6:448832/132‑3 (MQ=255)

                                                                                                                                               aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagggtc                                                                                                                                                                                                                                                                   >  5:448832/1‑130 (MQ=255)

                                                                                                                                                ttCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                                                                                                                                                                                                 >  6:323893/1‑149 (MQ=255)

                                                                                                                                                ttCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                                                                                                                                                   <  2:27427/163‑1 (MQ=255)

                                                                                                                                                  cGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                                                                                                                                                                               <  7:625012/149‑1 (MQ=255)

                                                                                                                                                     aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCAGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                                                                                                                            <  7:994115/149‑1 (MQ=255)

                                                                                                                                                       ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                                                                                                                                          <  5:327219/149‑1 (MQ=255)

                                                                                                                                                        cATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                                                                                                                                                              <  8:1200248/112‑1 (MQ=255)

                                                                                                                                                              aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                                                                                                                                                   <  5:1057395/149‑1 (MQ=255)

                                                                                                                                                               ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                                                                                                                                                  <  8:545351/149‑1 (MQ=255)

                                                                                                                                                               ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                                                                                                                                                  <  7:1303695/149‑1 (MQ=255)

                                                                                                                                                               ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                                                                                                                                                  <  5:397905/149‑1 (MQ=255)

                                                                                                                                                               ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                                                                                                                                                  <  7:253313/149‑1 (MQ=255)

                                                                                                                                                                tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                                                                                                                                                                                                                                 <  8:810721/149‑1 (MQ=255)

                                                                                                                                                                tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCtt                                                                                                                                                                                                                                 <  8:954028/149‑1 (MQ=255)

                                                                                                                                                                tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCTGGAGAAGGCCtt                                                                                                                                                                                                                                 <  8:722803/149‑1 (MQ=255)

                                                                                                                                                                       gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCNTCTCCTGGAGAAGGCCTTGCTATNAGTTNCGGt                                                                                                                                                                                                                  >  4:347583/1‑157 (MQ=255)

                                                                                                                                                                            taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                                                                                                                          >  5:1079432/1‑144 (MQ=255)

                                                                                                                                                                             aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                                                                                                                                    >  7:31837/1‑149 (MQ=255)

                                                                                                                                                                              aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                                                                                                                   >  7:1223564/1‑149 (MQ=255)

                                                                                                                                                                              aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                                                                                                                   >  6:1232638/1‑149 (MQ=255)

                                                                                                                                                                               gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                                                                                                                  >  7:773248/1‑149 (MQ=255)

                                                                                                                                                                               gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                                                                                                                  >  8:19442/1‑149 (MQ=255)

                                                                                                                                                                                       gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAacac                                                                                                                                                                                                                                                            >  5:63335/1‑99 (MQ=255)

                                                                                                                                                                                       gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAacac                                                                                                                                                                                                                                                            <  6:63335/99‑1 (MQ=255)

                                                                                                                                                                                          gCTATTTGACAGAGTTGAAACGAATTGCGGTTGCGGGTTACGGGCGCGCGGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAAAACTCAAGGCCTTCTCCTGGGGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGCTATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCcgcg                                                                          <  2:340491/274‑1 (MQ=255)

                                                                                                                                                                                          gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                                                                                                                                       <  4:150881/213‑1 (MQ=255)

                                                                                                                                                                                              tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatga                                                                                                                                                                                                   >  5:407963/1‑149 (MQ=255)

                                                                                                                                                                                               atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                                                                                                                                  <  7:412550/149‑1 (MQ=255)

                                                                                                                                                                                               atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                                                                                                                                  <  5:1481583/149‑1 (MQ=255)

                                                                                                                                                                                               atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                                                                                                                                  <  5:1232638/149‑1 (MQ=255)

                                                                                                                                                                                                tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                                                                                                                                                              <  5:731835/72‑1 (MQ=255)

                                                                                                                                                                                                tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                                                                                                                                 >  5:790895/1‑149 (MQ=255)

                                                                                                                                                                                                     gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata                                                                                                                                                                                            >  8:204121/1‑149 (MQ=255)

                                                                                                                                                                                                        atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                                                                                                                                                                                         >  6:1074856/1‑149 (MQ=255)

                                                                                                                                                                                                          gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGa                                                                                                                                                                                       >  7:1332058/1‑149 (MQ=255)

                                                                                                                                                                                                           acaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCGCCCACAGGGAGAGGCAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTATTTGCAGTACAGCGCCTTAATGATGTGATATCCGag                                                                                                                                                                                      <  6:400817/147‑2 (MQ=38)

                                                                                                                                                                                                           aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTTATGCCCTCACCCCGACCCTCTACCACAGGGAGAGGGAGAATACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAGGATGTGATATCCGaa                                                                                                                                                                                      <  8:1076069/149‑1 (MQ=255)

                                                                                                                                                                                                           aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGTGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                                                                                                                                                                                      >  5:307741/1‑149 (MQ=255)

                                                                                                                                                                                                           aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                                                                                                                                                                                                       <  7:779156/116‑1 (MQ=255)

                                                                                                                                                                                                           aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                                                                                                                                                                                      >  6:1218195/1‑149 (MQ=255)

                                                                                                                                                                                                            aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGCGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                                                                                                                                                     >  8:839848/1‑149 (MQ=255)

                                                                                                                                                                                                            aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                                                                                                                                                                      <  4:132054/116‑1 (MQ=255)

                                                                                                                                                                                                            aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                                                                                                                                     <  8:678722/117‑1 (MQ=255)

                                                                                                                                                                                                             aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACt                                                                                                                                                                                    >  7:1435366/1‑149 (MQ=255)

                                                                                                                                                                                                                tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg                                                                                                                                                                                 <  8:121420/149‑1 (MQ=255)

                                                                                                                                                                                                                 aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                                                                                                                                <  5:555182/149‑1 (MQ=255)

                                                                                                                                                                                                                 aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                                                                                                                                <  6:737400/149‑1 (MQ=255)

                                                                                                                                                                                                                  ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtg                                                                                                                                                                               <  8:1435366/149‑1 (MQ=255)

                                                                                                                                                                                                                    gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGtgtg                                                                                                                                                                             <  5:1218195/149‑1 (MQ=255)

                                                                                                                                                                                                                      tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCa                                                                                                                                                                           >  6:149772/1‑149 (MQ=255)

                                                                                                                                                                                                                        ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGACTTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTttcttc                                           >  2:143539/1‑275 (MQ=255)

                                                                                                                                                                                                                              gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                                                                                                                                                >  8:800129/1‑88 (MQ=255)

                                                                                                                                                                                                                              gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                                                                                                                                                <  7:800129/88‑1 (MQ=255)

                                                                                                                                                                                                                                 aaGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAGCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                                                                                                                                                <  5:420705/149‑1 (MQ=255)

                                                                                                                                                                                                                                 aaGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                                                                                                                                                <  8:1332058/149‑1 (MQ=255)

                                                                                                                                                                                                                                   gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg                                                                                                                                                              <  8:393051/149‑1 (MQ=255)

                                                                                                                                                                                                                                        cgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                                                                                                                         <  5:124888/149‑1 (MQ=255)

                                                                                                                                                                                                                                        cgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                                                                                                                         <  6:790895/149‑1 (MQ=255)

                                                                                                                                                                                                                                             aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                                                                                                                                                    <  5:1125332/149‑1 (MQ=255)

                                                                                                                                                                                                                                              cgtccTCACCCCGACCATCTCCCACAGGGAGATGTAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACAACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGccc                     <  1:46719/272‑1 (MQ=255)

                                                                                                                                                                                                                                                 ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGAc                                                                                                                                                >  5:1327048/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                                                                                                                  <  8:1107132/114‑1 (MQ=255)

                                                                                                                                                                                                                                                  cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                                                                                                                  >  7:1107132/1‑114 (MQ=255)

                                                                                                                                                                                                                                                     aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCCCCTTAATGATGTGTTATCCGAACTGGGTGTGCCGCGGGCCGGTCGGNCCCAGTNCCGGACAAGAGACAACCACTTTATCGAACGCCGGACCCAACTGCCCCCGGCGGAATTCACCTAAATc                                                                                   >  3:21919/1‑206 (MQ=255)

                                                                                                                                                                                                                                                     aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                                                                                                                                            <  7:495205/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                                                                                                                                            >  6:1313467/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                                                                                                                                                                                      <  6:1099389/105‑1 (MQ=255)

                                                                                                                                                                                                                                                          gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                                                                                                                         >  6:529981/1‑131 (MQ=255)

                                                                                                                                                                                                                                                          gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                                                                                                                                       >  7:740184/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           tcccTCTCCCTGTGGGAGAGGGTGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                                                                                                                      >  6:853237/2‑149 (MQ=255)

                                                                                                                                                                                                                                                             ccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                                                                                                                                            <  7:585793/141‑1 (MQ=255)

                                                                                                                                                                                                                                                             ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa                                                                                                                                    >  6:562557/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                                                                                                                  >  8:775213/1‑101 (MQ=255)

                                                                                                                                                                                                                                                                 tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAATAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAACTAAATCACCGCCGCGTTTGCCTGATGGCCAAATGAAGTGTTTCTTCGCAAGCTTGCCGACATCGGCCCCGTTCTTAATCTGCTCgag  >  2:332376/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                 tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACttt                                                                                                                                >  7:1204320/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa                                                                                                                               >  8:605787/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt                                                                                                                              >  6:1086596/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAGTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt                                                                                                                              <  5:258346/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATGTCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                                                                                                                             <  7:910875/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGt                                                                                                                                                       <  6:288638/123‑1 (MQ=255)

                                                                                                                                                                                                                                                                    cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                                                                                                                             >  5:927351/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                                                                                                                             >  5:195548/1‑149 (MQ=255)

CTGAGTACGGTAACTATCTGTTCTCTTACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAG  >  NC_000913/3959219‑3959750

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