BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,561,969	A→T	5.5%	intergenic (‑146/‑44)	glpE ← / → glpD	thiosulfate:cyanide sulfurtransferase (rhodanese)/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,561,969	0	A	T	5.5%	398.4 / 11.1	110	intergenic (‑146/‑44)	glpE/glpD	thiosulfate:cyanide sulfurtransferase (rhodanese)/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand): ref base A (58/46); new base T (3/3); total (61/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-02

Only 100 of 113 total aligned reads displayed.

AAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAAC  >  NC_000913/3561702‑3562223

aaaGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                                                                                                                                                         <  2:257722/275‑1 (MQ=255)

                                                                                                            tCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGGCGCCGCTGGa                                                                                                                                             >  2:238646/1‑275 (MQ=255)

                                                                                                            tCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCAAAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGa                                                                                                                                             >  2:367005/1‑275 (MQ=255)

                                                                                                                         tGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCa                                                                                                                                                                                >  4:128724/1‑227 (MQ=255)

                                                                                                                             tctctcTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCaa                                                                                                                                                                            <  2:128725/227‑1 (MQ=255)

                                                                                                                              ctctctTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAt                                                                                                                                                                                                                                                      >  2:452713/1‑152 (MQ=255)

                                                                                                                               tctctTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGAAGACGCCGCTGGACGCGGTTTATCCGTGCTGa                                                                                                                          >  2:512276/1‑275 (MQ=255)

                                                                                                                                 tctttcttTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                                                                                                                                                          <  7:624068/145‑1 (MQ=255)

                                                                                                                                 tctttcttTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                                                                                                                                                          >  8:624068/1‑145 (MQ=255)

                                                                                                                                  ctttctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTAt                                                                                                                                                                                                                                                  <  4:452712/152‑1 (MQ=255)

                                                                                                                                     tctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                                                                                                                                                   <  8:517172/148‑1 (MQ=255)

                                                                                                                                        ttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                                                                                                                                                          <  6:1080526/138‑1 (MQ=255)

                                                                                                                                        ttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                                                                                                                                                          >  5:1080526/1‑138 (MQ=255)

                                                                                                                                                 aGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACg                                                                                                                                                                                                                                       >  6:570167/1‑148 (MQ=255)

                                                                                                                                                   tGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                                                                                                                                                    <  6:1027280/133‑1 (MQ=255)

                                                                                                                                                   tGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                                                                                                                                                    >  5:1027280/1‑133 (MQ=255)

                                                                                                                                                       cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACa                                                                                                                                                                                                                                                       <  7:771178/126‑1 (MQ=255)

                                                                                                                                                       cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                                                                                                                                                    <  8:938994/129‑1 (MQ=255)

                                                                                                                                                       cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                                                                                                                                                    >  7:938994/1‑129 (MQ=255)

                                                                                                                                                       cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAtt                                                                                                                                                                                                                                                     <  7:469169/128‑1 (MQ=255)

                                                                                                                                                       cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAtt                                                                                                                                                                                                                                                     >  8:469169/1‑128 (MQ=255)

                                                                                                                                                       cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTg                                                                                                                                                                                                                                 <  7:160527/148‑1 (MQ=255)

                                                                                                                                                       cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTg                                                                                                                                                                                                                                 >  8:160527/1‑148 (MQ=255)

                                                                                                                                                         aaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGaa                                                                                                                                                                                                                               <  5:570167/148‑1 (MQ=255)

                                                                                                                                                          aTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTgc                                                                                                >  2:339882/1‑274 (MQ=255)

                                                                                                                                                                   acaGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgcagc                                                                                                                                                                                                                     >  7:553794/1‑148 (MQ=255)

                                                                                                                                                                        ttACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGa                                                                                                                                                                                                                <  8:553794/148‑1 (MQ=255)

                                                                                                                                                                             tCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATa                                                                                                                                                                                                                                                              >  7:304032/1‑97 (MQ=255)

                                                                                                                                                                             tCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATa                                                                                                                                                                                                                                                              <  8:304032/97‑1 (MQ=255)

                                                                                                                                                                                aaGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAt                                                                                                                                                                                                                                                      >  7:966165/1‑102 (MQ=255)

                                                                                                                                                                                 aGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCTGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCaac                                                                         >  4:388235/1‑273 (MQ=255)

                                                                                                                                                                                 aGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGCTCTCGCTTGCGCGACCGCTTCCGCCAGTTCaaa                                                                         >  2:7436/1‑274 (MQ=255)

                                                                                                                                                                                              cGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAt                                                                                                                                                                                                                                                      >  6:266785/1‑88 (MQ=255)

                                                                                                                                                                                              cGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAt                                                                                                                                                                                                                                                      <  5:266785/88‑1 (MQ=255)

                                                                                                                                                                                                 aagagacagataAATATAACTCTAAAATTTTTTTTCAAGTTTACCGAAAGCGCGATTCTTTGCTAATAGGTTCGATAACGAACATTTATGAGCTTTAACGAAGGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGGGCTGGTATCGCGGCAGACGCCGGTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGc                                                         <  4:330958/265‑1 (MQ=255)

                                                                                                                                                                                                   tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAg                                                                                                                                                                                                                                   <  7:695998/102‑1 (MQ=255)

                                                                                                                                                                                                   tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAg                                                                                                                                                                                                                                   >  8:695998/1‑102 (MQ=255)

                                                                                                                                                                                                   tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                                                                                                                                   >  8:782590/1‑134 (MQ=255)

                                                                                                                                                                                                   tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                                                                                                                                   <  7:782590/134‑1 (MQ=255)

                                                                                                                                                                                                    aaTGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAt                                                                                                                                                                                                                                                      >  8:3009/1‑82 (MQ=255)

                                                                                                                                                                                                    aaTGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAt                                                                                                                                                                                                                                                      <  7:3009/82‑1 (MQ=255)

                                                                                                                                                                                                         tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                                                                                                                                                         <  7:398154/74‑1 (MQ=255)

                                                                                                                                                                                                         tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                                                                                                                                                         >  8:398154/1‑74 (MQ=255)

                                                                                                                                                                                                         tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                                                                                                                                                         >  5:63027/1‑74 (MQ=255)

                                                                                                                                                                                                         tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                                                                                                                                                         <  6:63027/74‑1 (MQ=255)

                                                                                                                                                                                                         tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                                                                                                                                                    >  5:953570/1‑79 (MQ=255)

                                                                                                                                                                                                         tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                                                                                                                                                    <  6:953570/79‑1 (MQ=255)

                                                                                                                                                                                                          ataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTTCGCTAcc                                               >  2:475575/1‑275 (MQ=255)

                                                                                                                                                                                                            aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                                                                                                                                                   <  6:166563/77‑1 (MQ=255)

                                                                                                                                                                                                            aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                                                                                                                                                   >  5:166563/1‑77 (MQ=255)

                                                                                                                                                                                                            aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaa                                                                                                                                                                                                                                     <  5:592853/91‑1 (MQ=255)

                                                                                                                                                                                                              atatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaa                                                                                                                                                                                                                                     >  8:9130/1‑89 (MQ=255)

                                                                                                                                                                                                               tatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                                                                                                                                                       <  5:323796/118‑1 (MQ=255)

                                                                                                                                                                                                               tatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                                                                                                                                                       >  6:323796/1‑118 (MQ=255)

                                                                                                                                                                                                               tatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCTCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCCAAACTCATTCCCGTTGGCTTGCCCTTCCTTGa                                           >  4:416774/1‑274 (MQ=255)

                                                                                                                                                                                                                atCACTCTAAAATGTTTTTTCAATGTTAGCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgcag                                                                                                                                                                                                                      >  3:24059/1‑102 (MQ=255)

                                                                                                                                                                                                                atCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGGAGGAGCTCGCTTGCGCGACCTCTTCCGCAAGTTCAAAACTCATTCCCGGGGGCCTGCGCTACCTTGAg                                          >  4:18673/1‑274 (MQ=255)

                                                                                                                                                                                                                 tCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                                                                                                                                   >  6:976638/1‑120 (MQ=255)

                                                                                                                                                                                                                 tCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                                                                                                                                   <  5:976638/120‑1 (MQ=255)

                                                                                                                                                                                                                  cACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTTTCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCCGTTCAAAACTCATTCACGTTGGCCTGCGCTACCTTTAg                                          >  1:143977/1‑272 (MQ=255)

                                                                                                                                                                                                                  cACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTTCTGATGCTGGAGGCGCAGGCTCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAg                                          >  2:544378/1‑272 (MQ=255)

                                                                                                                                                                                                                    ctctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATg                                                                                                                                                                                                                  <  1:24059/102‑1 (MQ=255)

                                                                                                                                                                                                                     tctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGCATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACgg                                                             >  2:330959/1‑250 (MQ=255)

                                                                                                                                                                                                                      ctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                                                                                                                                                   <  8:1111225/67‑1 (MQ=255)

                                                                                                                                                                                                                      ctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                                                                                                                                                   >  7:1111225/1‑67 (MQ=255)

                                                                                                                                                                                                                      ctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAgg                                                                                                                                                                                                                          >  8:224180/1‑92 (MQ=255)

                                                                                                                                                                                                                        aaaaTGTTTTTTCAATGTTACCTAAAGCGTGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaa                                                                                                                                                                                                          <  7:772900/106‑1 (MQ=255)

                                                                                                                                                                                                                        aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaa                                                                                                                                                                                                          >  8:772900/1‑106 (MQ=255)

                                                                                                                                                                                                                        aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAggggg                                                                                                                                                                                     >  7:343603/1‑127 (MQ=255)

                                                                                                                                                                                                                        aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAggggg                                                                                                                                                                                     <  8:343603/127‑1 (MQ=255)

                                                                                                                                                                                                                        aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggcgg                                                                                                                                                                                  >  7:332300/1‑130 (MQ=255)

                                                                                                                                                                                                                          aaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                                                                                                                                   >  6:298512/1‑111 (MQ=255)

                                                                                                                                                                                                                          aaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                                                                                                                                   <  5:298512/111‑1 (MQ=255)

                                                                                                                                                                                                                                   tCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATc                                                                                                                                                                                                     >  5:1075895/1‑100 (MQ=255)

                                                                                                                                                                                                                                   tCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATc                                                                                                                                                                                                     <  6:1075895/100‑1 (MQ=255)

                                                                                                                                                                                                                                    cAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgtg                                                                                                                                                                                             >  6:396891/1‑107 (MQ=255)

                                                                                                                                                                                                                                    cAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgtg                                                                                                                                                                                             <  5:396891/107‑1 (MQ=255)

                                                                                                                                                                                                                                     aaTGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                                                                                                                                   <  8:888680/100‑1 (MQ=255)

                                                                                                                                                                                                                                     aaTGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                                                                                                                                   >  7:888680/1‑100 (MQ=255)

                                                                                                                                                                                                                                        gTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCt                                                                                                                                                >  7:1112901/1‑148 (MQ=255)

                                                                                                                                                                                                                                            ccTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCCATg                                                                                                                                                                          <  6:1014166/118‑1 (MQ=255)

                                                                                                                                                                                                                                            ccTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATg                                                                                                                                                                          >  5:1014166/1‑118 (MQ=255)

                                                                                                                                                                                                                                               aaaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGAcgcg                                                                                                                                         <  8:1112901/148‑1 (MQ=255)

                                                                                                                                                                                                                                                aaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                                                                                                                                                      >  5:1099468/1‑86 (MQ=255)

                                                                                                                                                                                                                                                aaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                                                                                                                                                      <  6:1099468/86‑1 (MQ=255)

                                                                                                                                                                                                                                                     cgATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggcg                                                                                                                                                                                   >  1:144711/1‑100 (MQ=255)

                                                                                                                                                                                                                                                       aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAgggg                                                                                                                                                                                      <  5:182097/95‑1 (MQ=255)

                                                                                                                                                                                                                                                       aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCt                                                                                                                                                                     >  8:85565/1‑112 (MQ=255)

                                                                                                                                                                                                                                                       aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCt                                                                                                                                                                     <  7:85565/112‑1 (MQ=255)

                                                                                                                                                                                                                                                       aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATcc                                                                                                                                 >  6:808641/1‑148 (MQ=255)

                                                                                                                                                                                                                                                       aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCATGCTGGCATGCGCTACCTTGAGCACTATGAATTCCGCCTGTTCACCGAGGCGCTGGCTGaaa  >  1:54368/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                tAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATa                                                                                                                                                                                          <  8:651870/82‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aaTATGTTCGATAACGAACCTTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGaaa                                                                                                                                                                                                              <  7:908231/61‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aaTATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGaaa                                                                                                                                                                                                              >  8:908231/1‑61 (MQ=255)

                                                                                                                                                                                                                                                                   aatGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                                                                                                                                                      >  5:266904/2‑67 (MQ=39)

                                                                                                                                                                                                                                                                   aatGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                                                                                                                                                      <  6:266904/66‑1 (MQ=39)

                                                                                                                                                                                                                                                                       ttCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaa                                                                                                                                                                                                          >  2:209142/1‑59 (MQ=255)

                                                                                                                                                                                                                                                                          gATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAAcc                                                                                                                                                                                                            <  6:647514/54‑1 (MQ=255)

                                                                                                                                                                                                                                                                           aTAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                                                                                                                                                       >  6:239312/1‑58 (MQ=255)

                                                                                                                                                                                                                                                                           aTAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                                                                                                                                                      <  4:209141/59‑1 (MQ=255)

AAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAAC  >  NC_000913/3561702‑3562223

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: `‑`