BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,229,322	C→G	5.8%	intergenic (+257/+116)	yohP → / ← dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,229,322	0	C	G	5.8%	468.2 / 17.3	137	intergenic (+257/+116)	yohP/dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand): ref base C (76/53); new base G (5/3); total (81/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

Only 100 of 138 total aligned reads displayed.

TAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAGCCCTTGCAATATCAGGGGAATTATTCAACACCCGAACATGCTGAAATAATTCCGTTGCTTCATCGTATTCTTTACGCAAATAACTCAACCACTGTTTAATCCGCGCAACG  >  NC_000913/2229069‑2229579

taattaatGTCATCAGGTCCGAAAATAACGAGAATCTTTCAGTCTCCAATCCGGTTGCGCTCCTGTCATGTGCATTGCTTCATATCATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCtt                                                                                                                                                                                                                                              <  4:272760/275‑1 (MQ=255)

                gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTCCCGGGTCGGGCGCGAGGGAACATCCCTGGAAAACGCGCCCTCAGCCCACACCCCATG‑GGGCTGCCCCGGCCTTCAGGGAACGCCCCgg                                                                                                                                                                                                                               >  2:332315/1‑274 (MQ=255)

                                          tctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCCCA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCCATTTTGACGCCCCCAAACAACCggt                                                                                                                                                                                                    >  2:195316/1‑273 (MQ=255)

                                                                               tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTg                                                                                                                                                                                                    <  2:52662/238‑1 (MQ=255)

                                                                                                          gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACAt                                                                                                                                                                         <  2:177665/238‑1 (MQ=255)

                                                                                                            aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCAt                                                                                                                                                                                                                                                                >  7:650334/1‑149 (MQ=32)

                                                                                                             gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCAtg                                                                                                                                                                                                                                                               >  8:158809/1‑149 (MQ=32)

                                                                                                                ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTgg                                                                                                                                                                                                                                                            <  5:105030/149‑1 (MQ=32)

                                                                                                                ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTgg                                                                                                                                                                                                                                                            <  7:1388327/149‑1 (MQ=32)

                                                                                                                ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTgg                                                                                                                                                                                                                                                            >  7:369952/1‑149 (MQ=32)

                                                                                                                 gcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TGGATGTggg                                                                                                                                                                                                                                                           <  7:961442/149‑1 (MQ=12)

                                                                                                                    gcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTg                                                                                                                                                                                                                                                        >  8:513774/1‑149 (MQ=32)

                                                                                                                     ccAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGc                                                                                                                                                                                                                                                       <  7:745058/149‑1 (MQ=32)

                                                                                                                      cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGc                                                                                                                                                                                                                                                          >  7:105920/1‑145 (MQ=32)

                                                                                                                      cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGcc                                                                                                                                                                                                                                                      >  6:653949/1‑149 (MQ=32)

                                                                                                                       aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TGGATGTGGGCTgagg                                                                                                                                                                                                                                                     <  5:178889/149‑4 (MQ=12)

                                                                                                                       aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGccc                                                                                                                                                                                                                                                     <  7:1340761/149‑1 (MQ=32)

                                                                                                                         tCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGc                                                                                                                                                                                                                                                          >  8:298484/1‑142 (MQ=32)

                                                                                                                            ccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGcc                                                                                                                                                                                                                                                >  8:916014/1‑149 (MQ=32)

                                                                                                                            ccgcccccCCCGGCACCCCCGGGTTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCTGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGaa                                                                                                                   <  4:195316/274‑1 (MQ=255)

                                                                                                                                  cccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAgg                                                                                                                                                                                                                                          <  7:513774/149‑1 (MQ=32)

                                                                                                                                  cccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAgg                                                                                                                                                                                                                                          <  5:248312/149‑1 (MQ=32)

                                                                                                                                   ccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAggg                                                                                                                                                                                                                                         >  6:746377/1‑149 (MQ=32)

                                                                                                                                    cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGa                                                                                                                                                                                                                                        >  7:482911/1‑149 (MQ=32)

                                                                                                                                    cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGa                                                                                                                                                                                                                                        >  6:112735/1‑149 (MQ=32)

                                                                                                                                      cTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAAc                                                                                                                                                                                                                                      <  8:369952/149‑1 (MQ=32)

                                                                                                                                            ccccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAg                                                                                                                                                                                                                                           >  5:1224661/1‑138 (MQ=32)

                                                                                                                                            ccccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCg                                                                                                                                                                                                                                >  8:1051478/1‑149 (MQ=32)

                                                                                                                                             cccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGCGGGCTGCCCCGGCCTTCAGGGAACGCCTCgg                                                                                                                                                                                                                               >  8:494799/1‑149 (MQ=17)

                                                                                                                                              ccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGc                                                                                                                                                                                                                              >  8:1250864/1‑149 (MQ=32)

                                                                                                                                              ccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGc                                                                                                                                                                                                                              >  7:904901/1‑149 (MQ=32)

                                                                                                                                                gggCTCTGGCGAACAAAAGCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAGCCCTTGCAATATCAgggg                                                                                              <  2:250003/275‑1 (MQ=255)

                                                                                                                                                  gCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAAtt                                                                                                                                                                                                                          <  5:833779/149‑1 (MQ=35)

                                                                                                                                                               aaaTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGccaccga                                                                                                                                                                                                             >  7:1189331/1‑147 (MQ=25)

                                                                                                                                                               aaaTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCaa                                                                                                                                                                                                             >  8:435712/1‑149 (MQ=35)

                                                                                                                                                                    gccgcTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAaacaac                                                                                                                                                                                                        >  6:471137/1‑149 (MQ=40)

                                                                                                                                                                            gcgGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGc                                                                                                                                                                                                <  6:457624/149‑1 (MQ=255)

                                                                                                                                                                              ggTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCt                                                                                                                                                                                              <  5:471137/149‑1 (MQ=255)

                                                                                                                                                                                    cTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAAc                                                                                                                                                                                                                                      >  5:119101/1‑103 (MQ=32)

                                                                                                                                                                                    cTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAAc                                                                                                                                                                                                                                      <  6:119101/103‑1 (MQ=32)

                                                                                                                                                                                    cTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGAt                                                                                                                                                                                        <  7:1051478/149‑1 (MQ=255)

                                                                                                                                                                                     tCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATa                                                                                                                                                                                       >  7:486178/1‑149 (MQ=255)

                                                                                                                                                                                     tCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATa                                                                                                                                                                                       >  8:1436088/1‑149 (MQ=255)

                                                                                                                                                                                       ggCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATaaa                                                                                                                                                                                     <  8:1152712/149‑1 (MQ=255)

                                                                                                                                                                                        gCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAAtttt                                                                                                                                                                                                                        >  7:363125/1‑113 (MQ=33)

                                                                                                                                                                                                cTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAaca                                                                                                                                                                            >  5:1428221/1‑149 (MQ=255)

                                                                                                                                                                                                cTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAaca                                                                                                                                                                            >  5:353565/1‑149 (MQ=255)

                                                                                                                                                                                                cTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAaca                                                                                                                                                                            >  5:534930/1‑149 (MQ=255)

                                                                                                                                                                                                 ttACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAacac                                                                                                                                                                           >  8:451990/1‑149 (MQ=255)

                                                                                                                                                                                                 ttACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAacac                                                                                                                                                                           >  8:1350956/1‑149 (MQ=255)

                                                                                                                                                                                                 ttACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAacac                                                                                                                                                                           >  7:1120422/1‑149 (MQ=255)

                                                                                                                                                                                                 ttACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAacac                                                                                                                                                                           >  7:1011305/1‑149 (MQ=255)

                                                                                                                                                                                                 ttACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACAt                                                                                                                                                            <  2:134751/164‑1 (MQ=255)

                                                                                                                                                                                                     ggCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTg                                                                                                                                                                       >  8:1441718/1‑149 (MQ=255)

                                                                                                                                                                                                     ggCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTAACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTg                                                                                                                                                                       >  7:1201381/1‑149 (MQ=255)

                                                                                                                                                                                                        tACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCa                                                                                                                                                                    >  5:450111/1‑149 (MQ=255)

                                                                                                                                                                                                                cGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACAt                                                                                                                                                            <  7:969394/149‑1 (MQ=255)

                                                                                                                                                                                                                 gggCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAAt                                                                                                                                                                                                                           >  6:815768/1‑85 (MQ=33)

                                                                                                                                                                                                                       gAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATg                                                                                                                                                     >  5:764471/1‑149 (MQ=255)

                                                                                                                                                                                                                         ggTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGt                                                                                                                                                   >  6:804327/1‑149 (MQ=255)

                                                                                                                                                                                                                          gTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtt                                                                                                                                                  >  6:954513/1‑149 (MQ=255)

                                                                                                                                                                                                                            aaCATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACAt                                                                                                                                                            <  7:192007/137‑1 (MQ=255)

                                                                                                                                                                                                                             aCATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGttttt                                                                                                                                               <  8:650334/149‑1 (MQ=255)

                                                                                                                                                                                                                               aTCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTa                                                                                                                                             >  5:925274/1‑149 (MQ=255)

                                                                                                                                                                                                                               aTCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTa                                                                                                                                             <  6:353565/149‑1 (MQ=255)

                                                                                                                                                                                                                                 cccTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTaga                                                                                                                                           <  5:1319285/149‑1 (MQ=255)

                                                                                                                                                                                                                                 cccTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTaga                                                                                                                                           >  5:1111557/1‑149 (MQ=255)

                                                                                                                                                                                                                                 cccTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTaga                                                                                                                                           <  8:874615/149‑1 (MQ=255)

                                                                                                                                                                                                                                 cccTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGGCaaaaaa                                                                                                                                                               >  8:1260381/1‑129 (MQ=255)

                                                                                                                                                                                                                                  ccTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTagag                                                                                                                                          >  8:413702/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGtt                                                                                                                                        >  6:77810/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGtt                                                                                                                                        >  6:171491/1‑149 (MQ=255)

                                                                                                                                                                                                                                     gTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGttt                                                                                                                                       <  8:1120422/149‑1 (MQ=255)

                                                                                                                                                                                                                                      tAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTtc                                                                                                                                      <  7:1350956/149‑1 (MQ=255)

                                                                                                                                                                                                                                      tAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTtc                                                                                                                                      >  6:156977/1‑149 (MQ=255)

                                                                                                                                                                                                                                       aaaaCGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAacac                                                                                                                                                                           >  7:238186/1‑111 (MQ=255)

                                                                                                                                                                                                                                       aaaaCGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAacac                                                                                                                                                                           <  8:238186/111‑1 (MQ=255)

                                                                                                                                                                                                                                        aaaCGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTtctc                                                                                                                                    >  6:142938/1‑149 (MQ=255)

                                                                                                                                                                                                                                         aaCGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCg                                                                                                                                   >  6:1198831/1‑149 (MQ=255)

                                                                                                                                                                                                                                            gcgcCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtttt                                                                                                                                                >  8:899260/1‑133 (MQ=255)

                                                                                                                                                                                                                                            gcgcCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtttt                                                                                                                                                <  7:899260/133‑1 (MQ=255)

                                                                                                                                                                                                                                               cccTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAGCCCTTGCAATATCAGGGGAATTATTc                                                                                      <  4:224548/188‑1 (MQ=255)

                                                                                                                                                                                                                                                ccTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCaa                                                                                                                            <  6:1308262/149‑1 (MQ=255)

                                                                                                                                                                                                                                                ccTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCaa                                                                                                                            >  6:440116/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 cTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAt                                                                                                                           >  6:31942/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  tCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAtt                                                                                                                          <  6:925274/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  tCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAtt                                                                                                                          <  7:158809/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aGCCCACA‑TGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGTCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGc                                                                                                                        >  6:381771/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGc                                                                                                                        >  7:1155159/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     gCCCACA‑TGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCt                                                                                                                                                                                <  8:919444/92‑1 (MQ=255)

                                                                                                                                                                                                                                                     gCCCACA‑TGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGcc                                                                                                                       >  6:203180/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     gCCCACA‑TGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGcc                                                                                                                       >  7:1025477/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     gCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATACCAATTGcc                                                                                                                       >  6:1381315/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      cccACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCt                                                                                                                      <  6:450111/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       ccACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTg                                                                                                                         >  8:462362/1‑145 (MQ=255)

                                                                                                                                                                                                                                                       ccACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTg                                                                                                                         <  7:462362/145‑1 (MQ=255)

                                                                                                                                                                                                                                                         aca‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGaa                                                                                                                   <  5:77810/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          ca‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAAt                                                                                                                  <  8:44797/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             tCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAtt                                                                                                                          <  5:1281884/139‑1 (MQ=255)

                                                                                                                                                                                                                                                             tCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAtt                                                                                                                          >  6:1281884/1‑139 (MQ=255)

TAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACA‑TCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAGCCCTTGCAATATCAGGGGAATTATTCAACACCCGAACATGCTGAAATAATTCCGTTGCTTCATCGTATTCTTTACGCAAATAACTCAACCACTGTTTAATCCGCGCAACG  >  NC_000913/2229069‑2229579

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