BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,610,107	G→A	E541K (GAG→AAG)	hyfG →	hydrogenase 4, subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,610,107	0	G	A	99.7%	1356.4 / NA	356	E541K (GAG→AAG)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (134/221); minor base C (0/1); total (134/222)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 366 total aligned reads displayed.

GGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTCGCG  >  NC_000913/2609961‑2610254

ggCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAaaag                                                                                                                                                   <  5:99168/149‑1 (MQ=255)

ggCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAaaag                                                                                                                                                   <  7:371311/149‑1 (MQ=255)

ggCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAaaga                                                                                                                                                  <  2:816776/150‑1 (MQ=255)

  ccGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCg                                                                                                                                               <  6:1652376/151‑1 (MQ=255)

   gtgtgtTGCGTTACATGCTGCGGGGCAATACCGTTTCTGACGCACCGCTGATTATGGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGGCAGTCAATAACCGTGCCGTATAAAAAGATCGa                                                                                                                                              <  6:516853/149‑1 (MQ=37)

   cGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCACAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCg                                                                                                                                               >  2:1336145/1‑150 (MQ=255)

    ggTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGaa                                                                                                                                             <  5:1350135/151‑1 (MQ=255)

    ggTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGa                                                                                                                                              >  1:353621/1‑150 (MQ=255)

       gtTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGc                                                                                                                                          <  7:289861/151‑1 (MQ=255)

        ttGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCt                                                                                                                                         >  2:1087251/1‑151 (MQ=255)

          gCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTa                                                                                                                                        <  3:530459/150‑1 (MQ=255)

            gtcaCATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACgg                                                                                                                                     <  5:942903/148‑1 (MQ=255)

               aCATGCTGCGCGGCAATACCGTTTCTGACGCAGCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCAt                                                                                                                                  <  8:723712/151‑1 (MQ=255)

                cATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCAtt                                                                                                                                 <  1:608640/151‑1 (MQ=255)

                cATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCAtt                                                                                                                                 >  1:406680/1‑151 (MQ=255)

                   gctgcGCGGCAATACCGTTTCTGACGCACCGCTGATTATGGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGAt                                                                                                                              <  6:82876/151‑1 (MQ=255)

                      gcgcGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCg                                                                                                                            <  6:64835/150‑1 (MQ=255)

                      gcgcGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGt                                                                                                                           <  3:1254451/151‑1 (MQ=255)

                           gCAATACCGTTTTTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCg                                                                                                                      <  1:1271037/151‑1 (MQ=255)

                             aaTACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGtt                                                                                                                    >  2:429218/1‑151 (MQ=255)

                             aaTACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGtt                                                                                                                    >  1:1527113/1‑151 (MQ=255)

                                aCCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcgc                                                                                                                 >  8:1091029/1‑151 (MQ=255)

                                   gTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcgccgc                                                                                                              <  3:676739/151‑1 (MQ=255)

                                    tttCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATCGATCGTAACCGTTCGCCGCt                                                                                                             >  8:1559794/1‑151 (MQ=255)

                                      tCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTg                                                                                                            <  3:176367/150‑1 (MQ=255)

                                       cTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGCACGCTACGGCATTGATCGTAACCGTTCGCCGCTGa                                                                                                           >  8:359968/1‑150 (MQ=255)

                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGa                                                                                                           <  2:1372320/148‑1 (MQ=255)

                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGAACGTAACCGTTCGCCGCTGAAGt                                                                                                        <  5:184576/151‑1 (MQ=255)

                                          aCGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTa                                                                                                       <  2:311793/151‑1 (MQ=255)

                                          aCGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCATTCGCCGCTGAAGTa                                                                                                       <  2:1573741/151‑1 (MQ=255)

                                           cGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTa                                                                                                       >  3:1568803/1‑150 (MQ=255)

                                             cACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAgg                                                                                                    >  7:201412/1‑151 (MQ=255)

                                              aCCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAgg                                                                                                    >  2:536955/1‑150 (MQ=255)

                                               ccGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGGGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGAc                                                                                                  <  1:1375163/151‑1 (MQ=255)

                                                cGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACa                                                                                                 <  5:712485/151‑1 (MQ=255)

                                                cGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACa                                                                                                 <  7:1500206/151‑1 (MQ=255)

                                                   tGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACaga                                                                                               >  8:1614646/1‑150 (MQ=255)

                                                     aTTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACaga                                                                                               >  8:64890/1‑148 (MQ=255)

                                                     aTTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGAGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGACGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAACTAAGGACaga                                                                                               >  2:665951/1‑148 (MQ=255)

                                                      ttATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACaga                                                                                               >  1:666124/1‑147 (MQ=255)

                                                       tATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGAt                                                                                           <  5:1584090/150‑1 (MQ=255)

                                                        aTCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATg                                                                                          <  3:82301/150‑1 (MQ=255)

                                                        aTCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATg                                                                                          <  5:1066923/150‑1 (MQ=255)

                                                        aTCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGc                                                                                         >  2:278647/1‑151 (MQ=255)

                                                          gggTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTg                                                                                       <  2:1304159/150‑1 (MQ=255)

                                                          cGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCt                                                                                        <  1:912105/150‑1 (MQ=255)

                                                              aGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAg                                                                                    <  7:785776/150‑1 (MQ=255)

                                                              aGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAg                                                                                    <  4:812084/150‑1 (MQ=255)

                                                                 cTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTAc                                                                                >  5:1197226/1‑151 (MQ=255)

                                                                    gATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTg                                                                              >  7:1379159/1‑150 (MQ=255)

                                                                     aTCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGa                                                                             >  2:1433395/1‑150 (MQ=255)

                                                                          tGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACt                                                                        >  7:1325637/1‑150 (MQ=255)

                                                                           gCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtat                                                                      <  4:470964/151‑1 (MQ=255)

                                                                           gCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtat                                                                      >  3:968797/1‑151 (MQ=255)

                                                                           gCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACt                                                                        <  2:1035379/149‑1 (MQ=255)

                                                                            cTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtatt                                                                     >  6:186770/1‑151 (MQ=255)

                                                                              aCTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtattat                                                                   <  7:588327/151‑1 (MQ=255)

                                                                              aCTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtatta                                                                    <  8:506848/150‑1 (MQ=255)

                                                                                 ccTTTACCGACGGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCTAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATAGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcg                                                                <  1:533602/151‑1 (MQ=255)

                                                                                 ccTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcg                                                                <  2:971538/151‑1 (MQ=255)

                                                                                 ccTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcg                                                                <  6:1030302/151‑1 (MQ=255)

                                                                                    gTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcgcgc                                                             >  6:309962/1‑151 (MQ=255)

                                                                                    gTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcgcg                                                              <  1:1147029/150‑1 (MQ=255)

                                                                                      accgaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCg                                                           <  6:943976/151‑1 (MQ=255)

                                                                                        cgaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCgg                                                          >  5:65042/1‑150 (MQ=255)

                                                                                        cgaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGa                                                         <  5:332094/151‑1 (MQ=255)

                                                                                        cgaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGa                                                         <  4:1083104/151‑1 (MQ=255)

                                                                                         gaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGaa                                                        <  4:1331092/151‑1 (MQ=255)

                                                                                           ccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAAcc                                                      <  1:232647/151‑1 (MQ=255)

                                                                                           ccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAAcc                                                      <  4:1455333/151‑1 (MQ=255)

                                                                                             gtgtGACGCTGGTAGATGTGCGCAAGCTCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGGATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAATCTATTATGCGCGCCGGAACcgc                                                    <  7:594883/151‑1 (MQ=255)

                                                                                                 gACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACg                                                >  3:377802/1‑151 (MQ=255)

                                                                                                  aCGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACgg                                               <  1:178499/151‑1 (MQ=255)

                                                                                                       ggTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGaa                                           >  5:563048/1‑150 (MQ=255)

                                                                                                          aGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGCTGCTGAAGTTACTGAAAACGATTATGCGCGCCGGAACCGCGACGGTGAAatat                                       >  8:1633099/1‑151 (MQ=255)

                                                                                                          aGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAatat                                       <  6:1441373/151‑1 (MQ=255)

                                                                                                          aGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAata                                        <  1:309181/150‑1 (MQ=255)

                                                                                                           gATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATc                                      <  4:1138766/151‑1 (MQ=255)

                                                                                                              gtgCGCAAGCGCCAGTCAAAAACGGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCt                                   <  8:1043741/151‑1 (MQ=255)

                                                                                                              gtgCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCt                                   <  4:1667870/151‑1 (MQ=255)

                                                                                                                 cgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTcg                                >  4:961583/1‑151 (MQ=255)

                                                                                                                       cgcCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGAGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCAc                          >  3:880462/1‑151 (MQ=255)

                                                                                                                        gcCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACt                         >  8:1327935/1‑151 (MQ=255)

                                                                                                                        cccATTCAAAAACCGTGCAGTATAAAAAGATCGAACGGTACGGAATTGATCGTAACCGTTCGCCGCTGAAGTAAGCACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGCAACCGCGACGGTGAAATATCCCTTCGCGCCACt                         <  6:1218234/150‑1 (MQ=255)

                                                                                                                           aGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTgg                       >  6:257688/1‑150 (MQ=255)

                                                                                                                            gTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTgg                       >  4:56941/1‑149 (MQ=255)

                                                                                                                              cAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAgg                    <  3:1123740/150‑1 (MQ=255)

                                                                                                                              cAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAgg                    <  1:1017766/150‑1 (MQ=255)

                                                                                                                              cAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGt                   >  1:819707/1‑151 (MQ=255)

                                                                                                                              cAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGt                   >  8:822490/1‑151 (MQ=255)

                                                                                                                                 aaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTc                  >  2:1205576/1‑149 (MQ=255)

                                                                                                                                   aCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGcc              >  6:1360576/1‑151 (MQ=255)

                                                                                                                                        gCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGACCCCCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGc         >  8:18821/1‑151 (MQ=255)

                                                                                                                                        cccGTATAAAAAGAGCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGc         <  8:1458539/150‑1 (MQ=255)

                                                                                                                                         ccGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCt        <  4:126809/151‑1 (MQ=255)

                                                                                                                                          cGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGTGCCACTGGAGGTCAGCCCTGGCtt       >  7:1181235/1‑151 (MQ=255)

                                                                                                                                          cGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCtt       >  8:1220449/1‑151 (MQ=255)

                                                                                                                                          cGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCtt       <  4:31013/151‑1 (MQ=255)

                                                                                                                                             ataAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTcg    <  7:992223/151‑1 (MQ=255)

                                                                                                                                               aaaaAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGCCGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTcgcg  >  7:1043316/1‑151 (MQ=255)

GGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTCGCG  >  NC_000913/2609961‑2610254

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