BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,406,125	C→G	5.2%	A70G (GCT→GGT)	accC →	biotin carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,406,125	0	C	G	5.2%	593.3 / 16.4	172	A70G (GCT→GGT)	accC	biotin carboxylase
Reads supporting (aligned to +/- strand): ref base C (95/68); new base G (7/2); total (102/70)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 182 total aligned reads displayed.

TGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCA  >  NC_000913/3405977‑3406273

         cTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCAcc                                                                                                                                            >  1:918619/1‑150 (MQ=255)

           gggCATCAAGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCgg                                                                                                                                          <  2:316618/150‑1 (MQ=255)

            ggCATCAAGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGgc                                                                                                                                         <  1:1722108/150‑1 (MQ=255)

             gCATCAAGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGgcg                                                                                                                                        <  2:1022514/150‑1 (MQ=255)

               aTCAAGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGGTGATTTCAGCGgcgct                                                                                                                                      <  1:400908/150‑2 (MQ=255)

               aTCAAGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGGTGATTTCAGCGgcgct                                                                                                                                      <  1:1680732/150‑2 (MQ=255)

                   aGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAg                                                                                                                                  <  2:1940248/150‑1 (MQ=255)

                   aGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAg                                                                                                                                  >  1:2358987/1‑150 (MQ=255)

                           gCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAatccat                                                                                                                          >  2:498900/1‑150 (MQ=255)

                             tgtgCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAatccatcc                                                                                                                        >  1:1396808/1‑150 (MQ=255)

                              gtgCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCg                                                                                                                       >  1:2437007/1‑150 (MQ=255)

                               tgCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCgg                                                                                                                      >  1:985048/1‑150 (MQ=255)

                                gCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCggg                                                                                                                     <  1:1106843/150‑1 (MQ=255)

                                 cACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGt                                                                                                                    >  2:275667/1‑150 (MQ=255)

                                     ccAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGCACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACg                                                                                                                >  1:894936/1‑150 (MQ=255)

                                     ccAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACg                                                                                                                >  1:374881/1‑150 (MQ=255)

                                        gcgcGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCt                                                                                                             <  2:1548145/150‑1 (MQ=255)

                                           cgGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCTGGTTACGGCTTcc                                                                                                          >  2:1739304/1‑150 (MQ=255)

                                           cgGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTcc                                                                                                          >  2:841818/1‑150 (MQ=255)

                                           cgGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTcc                                                                                                          <  2:1880515/150‑1 (MQ=255)

                                           cgGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTcc                                                                                                          >  1:1570711/1‑150 (MQ=255)

                                           cgGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTcc                                                                                                          <  1:1534569/150‑1 (MQ=255)

                                           cgGATCGCGAGCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTcc                                                                                                          <  2:1790418/150‑1 (MQ=255)

                                                 gcgATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCg                                                                                                    >  2:2008215/1‑150 (MQ=255)

                                                 gcgATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCg                                                                                                    <  1:469924/150‑1 (MQ=255)

                                                    aTCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCgaga                                                                                                 >  2:2375702/1‑150 (MQ=255)

                                                    aTCTAAAACACGTATTACTGGCAGATGAAACGGTCGGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCgaga                                                                                                 <  2:2321119/150‑1 (MQ=255)

                                                       tAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACg                                                                                              >  2:254310/1‑150 (MQ=255)

                                                         aaaCACGTATTACTGTCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGcc                                                                                            <  1:2474065/150‑1 (MQ=255)

                                                           acacGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCaa                                                                                          >  1:1812063/1‑150 (MQ=255)

                                                              cGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGgc                                                                                                                                         >  1:42048/1‑100 (MQ=255)

                                                              cGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACtt                                                                                       <  1:1636205/150‑1 (MQ=255)

                                                                   tACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCg                                                                                  <  2:1298780/150‑1 (MQ=255)

                                                                    aCTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGACCATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGACCGCCAACTTTGCCGa                                                                                 >  2:2434249/1‑150 (MQ=255)

                                                                    aCTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGa                                                                                 >  1:1023780/1‑150 (MQ=255)

                                                                     cTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAg                                                                                >  1:1667103/1‑150 (MQ=255)

                                                                       ggCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCa                                                                              >  2:381631/1‑150 (MQ=255)

                                                                       ggCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCAGGCAATCATCAGCGCCGCTGAAATCACCTGAGCAGTAGCAAACGATCCGGGAAACGGCATCCTCTCCGAGACCGGCAACTTTGCCGAGCa                                                                              >  2:1607194/1‑150 (MQ=255)

                                                                         cAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAgg                                                                            >  2:1372803/1‑150 (MQ=255)

                                                                         cAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAgg                                                                            >  1:851353/1‑150 (MQ=255)

                                                                         cAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAgg                                                                            >  1:652870/1‑150 (MQ=255)

                                                                           gATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGtt                                                                          <  2:583322/150‑1 (MQ=255)

                                                                                    gTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTcc                                                                 >  2:1192031/1‑150 (MQ=255)

                                                                                     tCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCg                                                                >  2:103461/1‑150 (MQ=255)

                                                                                     gcTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCg                                                                <  1:617897/149‑1 (MQ=255)

                                                                                        gTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCt                                                             >  2:1466771/1‑150 (MQ=255)

                                                                                         tATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCtt                                                            <  2:685370/150‑1 (MQ=255)

                                                                                         tATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCtt                                                            <  1:532496/150‑1 (MQ=255)

                                                                                          aTTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCttt                                                           >  1:1446812/1‑150 (MQ=255)

                                                                                             ggCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATc                                                        >  2:2206700/1‑150 (MQ=255)

                                                                                             ggCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATc                                                        <  2:6961/150‑1 (MQ=255)

                                                                                             ggCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATc                                                        >  1:2151709/1‑150 (MQ=255)

                                                                                             ggCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCGCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTAta                                                        >  1:1293909/1‑149 (MQ=255)

                                                                                               cccTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCtt                                                      >  1:1967597/1‑150 (MQ=255)

                                                                                               cccTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCtt                                                      >  2:1706636/1‑150 (MQ=255)

                                                                                                     tCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTgg                                                <  2:1396808/150‑1 (MQ=255)

                                                                                                      ccGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGc                                               >  1:2251722/1‑150 (MQ=255)

                                                                                                      ccGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGc                                               >  1:2260918/1‑150 (MQ=255)

                                                                                                       cGTCATTAAAAAGTTATGTGAACTTCCCGGCAATCATCAGCGCGGCTGAAATCACCGGCGCAGTAGCAATCCATGCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCGTTATCTTCAttttcc                                              <  2:818679/150‑5 (MQ=255)

                                                                                                        gTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGccc                                             <  1:85854/150‑1 (MQ=255)

                                                                                                        gTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGccc                                             <  2:785915/150‑1 (MQ=255)

                                                                                                         tCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCg                                            <  2:1023780/150‑1 (MQ=255)

                                                                                                            gTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGaaa                                         >  2:908496/1‑150 (MQ=255)

                                                                                                              aaaaaGTTATCTGAACATCACGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATACGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTACGGCTTTATCATCATTGGCACGAAAGc                                       >  2:2148316/1‑150 (MQ=255)

                                                                                                               aaaaGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCa                                      >  2:1952160/1‑150 (MQ=255)

                                                                                                                    ttATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAAc                                 <  1:1365120/150‑1 (MQ=255)

                                                                                                                    ttATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAAc                                 <  1:269756/150‑1 (MQ=255)

                                                                                                                    ttATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAAc                                 <  1:2008215/150‑1 (MQ=255)

                                                                                                                     tATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAAcc                                <  1:283117/150‑1 (MQ=255)

                                                                                                                       tCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCAt                              <  2:1671051/150‑1 (MQ=255)

                                                                                                                       tCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCAt                              >  1:813738/1‑150 (MQ=255)

                                                                                                                        cTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCAtt                             <  1:1452679/150‑1 (MQ=255)

                                                                                                                         tGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTc                            >  2:1099723/1‑150 (MQ=255)

                                                                                                                          gAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCg                           >  1:1525795/1‑150 (MQ=255)

                                                                                                                          gAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCg                           <  2:1727090/150‑1 (MQ=255)

                                                                                                                          gAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCg                           >  1:949741/1‑150 (MQ=255)

                                                                                                                             cATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCt                        >  2:427156/1‑150 (MQ=255)

                                                                                                                                cccGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGAt                     <  2:2060594/150‑1 (MQ=255)

                                                                                                                                 ccGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATg                    <  2:1436133/150‑1 (MQ=255)

                                                                                                                                 ccGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATg                    >  1:695287/1‑150 (MQ=255)

                                                                                                                                 ccGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATg                    >  1:1003655/1‑150 (MQ=255)

                                                                                                                                   ggCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATggg                  >  1:785746/1‑150 (MQ=255)

                                                                                                                                     cAATCATCAGCGCCGCTGAAATCAGCGGCGCAGTAGCAATCCATCCGTGTTACGGCTTACTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACAATTCGCCTGATGGGCg                <  1:725436/150‑1 (MQ=255)

                                                                                                                                     cAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATATTCATTGGCACGAAATCAGAATCCATTCGCCTGATGGGct                <  2:1262324/150‑2 (MQ=255)

                                                                                                                                      aaTCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGATCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGTGCGa               <  2:976774/150‑1 (MQ=255)

                                                                                                                                       tgctactgcgcCGGTGATTTCAGCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGAc              >  2:95108/8‑150 (MQ=255)

                                                                                                                                       atcatcAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGAc              >  2:2100721/1‑150 (MQ=255)

                                                                                                                                         ctactgcgcCGGTGATTTCAGCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACaa            >  1:2302185/6‑150 (MQ=255)

                                                                                                                                         catcaGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACaa            <  2:1354107/150‑1 (MQ=255)

                                                                                                                                          atcaGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACaaa           <  2:2449719/150‑1 (MQ=255)

                                                                                                                                           tcaGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCTGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAg          <  1:1739304/150‑1 (MQ=255)

                                                                                                                                           tcaGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGGTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAg          <  1:1742162/150‑1 (MQ=255)

                                                                                                                                           actgcgcCGGTGATTTCAGCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAg          >  2:511826/4‑150 (MQ=255)

                                                                                                                                            ctgcgcCGGTGATTTCAGCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGt         >  1:2060544/3‑150 (MQ=255)

                                                                                                                                             tgcgcCGGTGATTTCAGCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTa        >  1:1935803/2‑150 (MQ=255)

                                                                                                                                               cgccggTGATTTCAGCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATc      >  2:1215161/1‑150 (MQ=255)

                                                                                                                                               cgccgcTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATc      >  2:1205081/1‑150 (MQ=255)

                                                                                                                                               cgccgcTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATc      <  2:802233/150‑1 (MQ=255)

                                                                                                                                                 ccgcTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCg    <  1:1205081/150‑1 (MQ=255)

                                                                                                                                                   gcTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCa  >  1:2075631/1‑150 (MQ=255)

TGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCGGATCGCGATCTAAAACACGTATTACTGGCAGATGAAACGGTCTGTATTGGCCCTGCTCCGTCAGTAAAAAGTTATCTGAACATCCCGGCAATCATCAGCGCCGCTGAAATCACCGGCGCAGTAGCAATCCATCCGGGTTACGGCTTCCTCTCCGAGAACGCCAACTTTGCCGAGCAGGTTGAACGCTCCGGCTTTATCTTCATTGGCCCGAAAGCAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCA  >  NC_000913/3405977‑3406273

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