BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,561,972	A→T	5.4%	intergenic (‑149/‑41)	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,972	0	A	T	5.4%	413.6 / 11.2	111	intergenic (‑149/‑41)	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/47); new base T (3/3); total (61/50)
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.81e-03

Only 100 of 112 total aligned reads displayed.

AAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCC  >  NC_000913/3561702‑3562246

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

  aGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACa                                                                                                                                                                                                                                                                              >  2:259183/1‑275 (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                                                                                                                                                                                                                                                                                 >  5:1080526/1‑138 (MQ=255)

                                                                                                                                        ttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                                                                                                                                                                                 <  6:1080526/138‑1 (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)

                                                                                                                                                       cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACa                                                                                                                                                                                                                                                                              >  8:771178/1‑126 (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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                                                                                            aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaa                                                                                                                                                                                                                                                            >  6:592853/1‑91 (MQ=255)

                                                                                                                                                                                                              atatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaa                                                                                                                                                                                                                                                            <  7:9130/89‑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)

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

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

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

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

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

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

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

                                                                                                                                                                                                                      ctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAgg                                                                                                                                                                                                                                                 <  7:224180/92‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                                                                                                        aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggcgg                                                                                                                                                                                                         <  8:332300/130‑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                                                                                                                                                                                                                    <  5:396891/107‑1 (MQ=255)

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

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

                                                                                                                                                                                                                                     aaTGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                                                                                                                                                          <  8:888680/100‑1 (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)

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

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

                                                                                                                                                                                                                                                       aTTCTTTGCTAATATGTTCGATAACGAACATGTATGAGCTGTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAgggg                                                                                                                                                                                                             >  6:182097/1‑95 (MQ=255)

                                                                                                                                                                                                                                                         tCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCAt                                                                                                                                                                                                      <  3:144711/100‑1 (MQ=255)

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

                                                                                                                                                                                                                                                                tAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATa                                                                                                                                                                                                                 >  7:651870/1‑82 (MQ=255)

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

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

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

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

                                                                                                                                                                                                                                                                          gATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAAcc                                                                                                                                                                                                                                   >  5:647514/1‑54 (MQ=255)

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

                                                                                                                                                                                                                                                                           aTAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                                                                                                                                                                              <  5:239312/58‑1 (MQ=255)

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

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

                                                                                                                                                                                                                                                                              aCGAACATGTATGAGCTTTAACGAAATTGAATGAGGGCAGCTTTGATACCATAGTTCTGTTTGTGATAGGGGGCGCCATCAATGGTTCTGGTCTTGTGGCACACGCCGCTGGTAGCGGTTTATCCGGCCTGAGGCTGGGGGCGCAGTATCTCGCTTGCGCGACCTCTACCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTAGGAATTCCGCCGGGTCAGCGAGGGGCTGGCTGAACGTGACGTGCTGCTGAAAATGGcc  <  4:437552/275‑1 (MQ=255)

AAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCC  >  NC_000913/3561702‑3562246

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