BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,113,544	G→A	5.9%	intergenic (+138/‑35)	opgH → / → yceK	OPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,113,544	0	G	A	5.9%	707.2 / 21.5	221	intergenic (+138/‑35)	opgH/yceK	OPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein
Reads supporting (aligned to +/- strand): ref base G (113/92); new base A (5/8); total (120/101)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

Only 100 of 221 total aligned reads displayed.

GGCCCGTCTGCATTTCCGTGTCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATT  >  NC_000913/1113286‑1113808

                    tCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCACTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGAATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAAt                                                                                                                                                                                                                                      >  2:291995/1‑275 (MQ=255)

                    tCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAATACTTTGTGGGTCGGATTAGTCGTTCACGCCGCATCCGGCAAAAACCACGAGCCAATAAATATACCGGCCTTATGCCGGTATTTTTTTCCGAAAGAGTTATc                                                                                                                                                                                                                                          >  2:284282/1‑271 (MQ=255)

                                   tgggaTATTCTTTATGGGTGAGTAATTACGAAGGGAGAAAGTTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAATGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggt                                                                                                                                                                                                                       <  4:274981/272‑1 (MQ=255)

                                   gtggtttttcttaATGGGTGAGTTATTACGAAGGGATAAAGCTTAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAACCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggt                                                                                                                                                                                                                       <  4:291995/268‑1 (MQ=255)

                                                       aGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATAGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTAACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCACATGAGTTAAATTGTGTTGACCATCATGGTGACCTTGc                                                                                                                                                                                                   >  2:55405/1‑275 (MQ=255)

                                                                   gggATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGAGAAGGCGTTCACGACGCATCCGGCAAAAACAACGAGCCAATAAAAATACAGGCGTAATGCGGGTATTTTTTTACGCAAGAGGATTCCAATGTGTTTAATTGTGTTTAGCATGATGGTGACCTTCATGAGTGGAtgtg                                                                                                                                                                                       >  2:228100/1‑275 (MQ=38)

                                                                      aTAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTAACaaa                                                                                                                                                                                                                                       >  4:330848/1‑224 (MQ=255)

                                                                            cTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGtt                                                                                                                                                                                                                                                                        >  2:36681/1‑185 (MQ=255)

                                                                               aaTCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGc                                                                                                                                                                                                                                                                    >  2:345518/1‑186 (MQ=255)

                                                                                aTCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGc                                                                                                                                                                                                                                                                    <  4:36681/185‑1 (MQ=255)

                                                                                   cACTGGCATTGCGTAAACCGGATGCGGCTTCGANATAANAACGTAGTTGCCTGATGCGCTACNCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGt                                                                                                                                                                                                                                                                <  4:345518/186‑1 (MQ=255)

                                                                                                              cTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCa                                                                                                                                             <  2:169410/274‑1 (MQ=255)

                                                                                                               ttCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGt                                                                                                                                                                                                                                                                         <  7:1080888/149‑1 (MQ=255)

                                                                                                               ttCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGt                                                                                                                                                                                                                                                                         >  7:1392851/1‑149 (MQ=255)

                                                                                                                 cGCAATAAAAACGTAGTTGCCTGATGCGTTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTa                                                                                                                                                                                                                                                                       >  7:1253042/1‑149 (MQ=255)

                                                                                                                    aaTAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGc                                                                                                                                                                                                                                                                    >  7:865455/1‑149 (MQ=255)

                                                                                                                      tAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGtt                                                                                                                                                                                                                                                                        <  6:1175745/143‑1 (MQ=255)

                                                                                                                      tAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCg                                                                                                                                                                                                                                                                  >  5:276793/1‑149 (MQ=255)

                                                                                                                      tAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCg                                                                                                                                                                                                                                                                  >  6:278629/1‑149 (MQ=255)

                                                                                                                         aaaCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTa                                                                                                                                                                                                                                                               >  6:1153362/1‑149 (MQ=255)

                                                                                                                         aaaCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTa                                                                                                                                                                                                                                                               <  7:840141/149‑1 (MQ=255)

                                                                                                                           aCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAtt                                                                                                                                                                                                                                                             >  6:1328532/1‑149 (MQ=255)

                                                                                                                             gTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAtttt                                                                                                                                                                                                                                                           >  7:343338/1‑149 (MQ=255)

                                                                                                                              tAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAttttt                                                                                                                                                                                                                                                          <  5:575096/149‑1 (MQ=255)

                                                                                                                                 ttGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACACCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                                                                                                                                                                                                       >  5:1054499/1‑149 (MQ=255)

                                                                                                                                   gCCTGATTCGCTACGCTTATCAGCCCTAAATCGTTCCTGCAATTTATTGATTTTGCAATACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTGTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCa                                                                                                                       <  1:65148/275‑1 (MQ=255)

                                                                                                                                   gCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttt                                                                                                                                                                                                                                                            <  8:70193/142‑1 (MQ=255)

                                                                                                                                   gCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttt                                                                                                                                                                                                                                                            >  7:70193/1‑142 (MQ=255)

                                                                                                                                    ccTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGa                                                                                                                                                                                                                                                    >  8:248170/1‑149 (MQ=255)

                                                                                                                                      tGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGaaa                                                                                                                                                                                                                                                  <  7:185890/149‑1 (MQ=255)

                                                                                                                                        aTGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCg                                                                                                                                                                       <  4:123134/222‑1 (MQ=255)

                                                                                                                                         tGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                                                                                                                                                                                               <  8:804111/149‑1 (MQ=255)

                                                                                                                                             cTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAt                                                                                                                                                                                                                                                              >  5:1364985/1‑130 (MQ=255)

                                                                                                                                             cTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAt                                                                                                                                                                                                                                                              <  6:1364985/130‑1 (MQ=255)

                                                                                                                                              tACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATc                                                                                                                                                                                                                                          >  8:1336516/1‑149 (MQ=255)

                                                                                                                                              tACGCTTATCACGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCCGTATTTTTTTa                                                                                                                                                                                                                                                       <  7:961846/136‑1 (MQ=255)

                                                                                                                                                   ttATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATg                                                                                                                                                                                                                                     >  7:670251/1‑149 (MQ=255)

                                                                                                                                                    tATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGc                                                                                                                                                                                                                                    <  7:1427902/149‑1 (MQ=255)

                                                                                                                                                    tATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGc                                                                                                                                                                                                                                    <  6:1477482/149‑1 (MQ=255)

                                                                                                                                                     aTCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGtt                                                                                                                                                                                                                                                                        <  6:166127/112‑1 (MQ=255)

                                                                                                                                                     aTCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTCTAATTgtggt                                                                                                                                                                                                                       >  2:298272/1‑161 (MQ=255)

                                                                                                                                                      tCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGt                                                                                                                                                                                                                                  >  5:150669/1‑149 (MQ=255)

                                                                                                                                                       cAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGtt                                                                                                                                                                                                                                 <  7:614815/149‑1 (MQ=255)

                                                                                                                                                        aGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGttt                                                                                                                                                                                                                                >  7:94050/1‑149 (MQ=255)

                                                                                                                                                            cTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAGTGCGTTTAAtt                                                                                                                                                                                                                            <  5:1096900/149‑1 (MQ=255)

                                                                                                                                                            cTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAAtt                                                                                                                                                                                                                            <  6:17913/149‑1 (MQ=255)

                                                                                                                                                              aCATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgt                                                                                                                                                                                                                          <  5:1153362/149‑1 (MQ=255)

                                                                                                                                                              aCATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgt                                                                                                                                                                                                                          <  7:1355434/149‑1 (MQ=255)

                                                                                                                                                              aCATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAAtatta                                                                                                                                   <  4:172816/236‑1 (MQ=255)

                                                                                                                                                               cATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgtg                                                                                                                                                                                                                         <  7:493535/149‑1 (MQ=255)

                                                                                                                                                               cATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttt                                                                                                                                                                                                                                                            <  6:1083652/114‑1 (MQ=37)

                                                                                                                                                                aTCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtgg                                                                                                                                                                                                                        >  5:836462/1‑149 (MQ=255)

                                                                                                                                                                aTCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtgg                                                                                                                                                                                                                        >  6:5292/1‑149 (MQ=255)

                                                                                                                                                                aTCGTTCCTGCAATTTATTGATTTTGCAAGACTGTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTGTAATTgtgg                                                                                                                                                                                                                        >  6:972784/1‑149 (MQ=255)

                                                                                                                                                                    ttCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGt                                                                                                                                                                         >  1:52778/1‑192 (MQ=255)

                                                                                                                                                                        tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGcatc                                                                                                                                                                                                                >  5:239539/1‑149 (MQ=255)

                                                                                                                                                                         gCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTc                                                                                                                                                                        >  1:31547/1‑188 (MQ=255)

                                                                                                                                                                            aTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAaga                                                                                                                                                                                                                                                <  8:835418/113‑1 (MQ=255)

                                                                                                                                                                            aTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATgg                                                                                                                                                                                                            >  8:197171/1‑149 (MQ=255)

                                                                                                                                                                             tttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGt                                                                                                                                                                                                           <  8:1392851/149‑1 (MQ=255)

                                                                                                                                                                             tttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCAc                                                                                                                                                                    <  3:31547/188‑1 (MQ=255)

                                                                                                                                                                                aTTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAc                                                                                                                                                                                                        >  5:1193679/1‑149 (MQ=255)

                                                                                                                                                                                aTTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGAAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAc                                                                                                                                                                                                        <  5:1018899/149‑1 (MQ=255)

                                                                                                                                                                                  tGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCt                                                                                                                                                                                                      <  5:410673/149‑1 (MQ=255)

                                                                                                                                                                                    aTTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTg                                                                                                                                                                                                    >  5:335339/1‑149 (MQ=255)

                                                                                                                                                                                        tGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGa                                                                                                                                                                                                <  6:1040461/149‑1 (MQ=255)

                                                                                                                                                                                          cAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                                                                                                                                              >  8:750833/1‑149 (MQ=255)

                                                                                                                                                                                          cAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                                                                                                                                              <  8:1438864/149‑1 (MQ=255)

                                                                                                                                                                                           aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAAtt                                                                                                                                                                                                                            <  7:1234266/118‑1 (MQ=255)

                                                                                                                                                                                           aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTg                                                                                                                                                                                             >  7:1225434/1‑149 (MQ=255)

                                                                                                                                                                                            aGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGa                                                                                                                                                                                                <  8:877971/145‑1 (MQ=255)

                                                                                                                                                                                             gACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGc                                                                                                                                                                                           >  7:219220/1‑149 (MQ=255)

                                                                                                                                                                                              aCTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCt                                                                                                                                                                                          <  5:5292/149‑1 (MQ=255)

                                                                                                                                                                                                tttGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                                                                                                                                                                                               >  6:86378/1‑94 (MQ=255)

                                                                                                                                                                                                tttGGGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                                                                                                                                                                                               <  5:86378/94‑1 (MQ=255)

                                                                                                                                                                                                  tgtgGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTgg                                                                                                                                                                                      <  5:517010/149‑1 (MQ=255)

                                                                                                                                                                                                      ggTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCAc                                                                                   <  4:45289/244‑1 (MQ=255)

                                                                                                                                                                                                       gTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggt                                                                                                                                                                                                                       <  8:448647/111‑1 (MQ=255)

                                                                                                                                                                                                           gATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCattat                                                                                                                                                                             >  7:428086/1‑149 (MQ=255)

                                                                                                                                                                                                             tAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATTATGCCGCTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTTACGTTGCTGAGTGGCTGTGGCAGCAttatta                                                                                                                                                                           <  8:1211472/149‑1 (MQ=255)

                                                                                                                                                                                                              aaGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAg                                                                                                                                                                          <  8:885723/149‑1 (MQ=255)

                                                                                                                                                                                                                  agTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGc                                                                                                                                                                      >  6:466936/2‑149 (MQ=255)

                                                                                                                                                                                                                        cgccgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACtata                                                                                                                                                                <  6:559320/149‑1 (MQ=255)

                                                                                                                                                                                                                        cgccgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACtata                                                                                                                                                                <  5:305053/149‑1 (MQ=255)

                                                                                                                                                                                                                                     aaaaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcaggggc                                                                                                                                                   <  7:944839/149‑1 (MQ=255)

                                                                                                                                                                                                                                     aaaaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcaggggc                                                                                                                                                   >  6:598066/1‑149 (MQ=255)

                                                                                                                                                                                                                                     aaaaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcaggggc                                                                                                                                                   <  7:197171/149‑1 (MQ=255)

                                                                                                                                                                                                                                      aaaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcaggggca                                                                                                                                                  >  5:633239/1‑149 (MQ=255)

                                                                                                                                                                                                                                       aaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCAt                                                                                                                                                 <  7:1336516/149‑1 (MQ=255)

                                                                                                                                                                                                                                           aacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCa                                                                                                                                             >  6:225471/1‑149 (MQ=255)

                                                                                                                                                                                                                                           aacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCa                                                                                                                                             >  8:1132869/1‑149 (MQ=255)

                                                                                                                                                                                                                                                gCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCaa                                                                                                                                        >  7:107820/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcag                                                                                                                                                       >  6:1265379/1‑131 (MQ=255)

                                                                                                                                                                                                                                                   aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcag                                                                                                                                                       <  5:1265379/131‑1 (MQ=255)

                                                                                                                                                                                                                                                   aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATg                                                                                                                                                >  8:60324/1‑138 (MQ=255)

                                                                                                                                                                                                                                                     tAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAAtatta                                                                                                                                   >  5:679634/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       aaaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATc                                                                                                                                 <  6:426862/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        aaaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                                                                                                                                              >  7:531363/1‑87 (MQ=38)

                                                                                                                                                                                                                                                        aaaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                                                                                                                                              <  8:531363/87‑1 (MQ=38)

                                                                                                                                                                                                                                                          aTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCt                                                                                                                                                                                                  >  5:1411836/1‑81 (MQ=255)

GGCCCGTCTGCATTTCCGTGTCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATT  >  NC_000913/1113286‑1113808

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