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	93.9%	1255.8 / 45.3	363	E541K (GAG→AAG)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand): ref base G (11/11); new base A (145/195); total (156/207)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 369 total aligned reads displayed.

GGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTCGCG  >  NC_000913/2609961‑2610254

ggCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGAt                                                                                                                                                 <  4:327569/151‑1 (MQ=255)

ggCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGAt                                                                                                                                                 <  6:215893/151‑1 (MQ=255)

    ggTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTTCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGaa                                                                                                                                             <  5:1051859/151‑1 (MQ=255)

    ggTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGaa                                                                                                                                             <  2:544137/151‑1 (MQ=255)

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

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

     gtgtTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAAc                                                                                                                                            <  3:1129254/151‑1 (MQ=255)

       gtTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACg                                                                                                                                           <  8:823611/150‑1 (MQ=255)

         tGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCt                                                                                                                                         <  5:172162/150‑1 (MQ=255)

          gCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTa                                                                                                                                        <  2:818985/150‑1 (MQ=255)

           cGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACg                                                                                                                                      <  8:870924/151‑1 (MQ=255)

            gTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTTAAAAACCGTGCCGTATAAAAAGATCGAACGCTACg                                                                                                                                      <  5:1284826/150‑1 (MQ=255)

              tACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGc                                                                                                                                    <  8:721872/150‑1 (MQ=255)

               aCATGCTGCGCGGCAATACCGTTTCTGACGCACGGCTGAATATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCAGGTAGATGAGCGCAAGCGACACTCAAAAACCGTCCCGTATAAAAAGATCGAACGCTACGCCa                                                                                                                                   >  1:545514/1‑150 (MQ=255)

               aCATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCAt                                                                                                                                  >  1:1247320/1‑151 (MQ=255)

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

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

                   gctgcGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGAt                                                                                                                              >  3:1005020/1‑151 (MQ=255)

                   gctgcGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGAt                                                                                                                              >  6:589894/1‑151 (MQ=255)

                      gcgcGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGt                                                                                                                           <  4:926494/151‑1 (MQ=255)

                      gcgcGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGt                                                                                                                           <  5:557433/151‑1 (MQ=255)

                         cgGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAc                                                                                                                        >  8:1137556/1‑151 (MQ=255)

                            cAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGt                                                                                                                     >  3:795647/1‑151 (MQ=255)

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

                              aTACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAACGCTACGGCATTGATCGTAACCGTTc                                                                                                                   <  4:1212322/151‑1 (MQ=255)

                              aTACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTc                                                                                                                   <  2:1417629/151‑1 (MQ=255)

                                 ccGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcgcc                                                                                                                <  1:942392/151‑1 (MQ=255)

                                    tttCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCt                                                                                                             >  2:1265565/1‑151 (MQ=255)

                                     ttCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTg                                                                                                            >  6:392304/1‑151 (MQ=255)

                                     ttCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGGGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCt                                                                                                             <  1:1445839/150‑1 (MQ=255)

                                        tGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGa                                                                                                           >  2:929093/1‑149 (MQ=255)

                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGt                                                                                                        <  5:1391207/151‑1 (MQ=255)

                                         gACGCACCGCTGAGTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGt                                                                                                        <  8:345761/151‑1 (MQ=255)

                                          aCGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTa                                                                                                       <  3:75675/151‑1 (MQ=255)

                                           cGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTaa                                                                                                      <  8:375201/151‑1 (MQ=255)

                                           cGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTaa                                                                                                      <  1:764604/151‑1 (MQ=255)

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

                                              aCCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAgg                                                                                                    >  3:1044702/1‑150 (MQ=255)

                                              aCCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGa                                                                                                   <  2:277735/151‑1 (MQ=255)

                                              aCCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCCTTCGCCGCTGAAGTAAgg                                                                                                    >  7:439251/1‑150 (MQ=255)

                                                     aTTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACaga                                                                                               >  4:878223/1‑148 (MQ=255)

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

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

                                                           ggTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGa                                                                                      <  7:583577/151‑1 (MQ=255)

                                                           ggTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGa                                                                                      <  7:865390/151‑1 (MQ=255)

                                                            gTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGaa                                                                                     <  4:510101/151‑1 (MQ=255)

                                                            gTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGaa                                                                                     <  6:1285692/151‑1 (MQ=255)

                                                            gTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGa                                                                                      >  4:579080/1‑150 (MQ=255)

                                                              aGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGt                                                                                   <  6:1223675/151‑1 (MQ=255)

                                                                  ttGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACt                                                                               >  5:803604/1‑151 (MQ=255)

                                                                    gATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTg                                                                              <  6:1165242/150‑1 (MQ=255)

                                                                    gATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACAGa                                                                             <  1:219423/151‑1 (MQ=255)

                                                                          tGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACt                                                                        >  5:694459/1‑150 (MQ=255)

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

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

                                                                           gCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACt                                                                        <  6:1012702/149‑1 (MQ=255)

                                                                           gCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACt                                                                        <  3:449115/149‑1 (MQ=255)

                                                                             tACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtatt                                                                     >  7:1002743/1‑150 (MQ=255)

                                                                              aCTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGCTGCTGAAGTTACTGAAAAGtattat                                                                   >  7:1056646/1‑151 (MQ=255)

                                                                              aCTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtatt                                                                     <  4:636909/149‑1 (MQ=255)

                                                                                tCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgc                                                                 >  7:703882/1‑151 (MQ=255)

                                                                                 ccTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATg                                                                  <  1:1097929/149‑1 (MQ=255)

                                                                                 ccTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATg                                                                  <  4:815281/149‑1 (MQ=255)

                                                                                      accgaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAATGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCTCCg                                                           <  4:764364/151‑1 (MQ=255)

                                                                                           ccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAAcc                                                      >  8:1389463/1‑151 (MQ=255)

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

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

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

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

                                                                                               gtgACTCTGGTAGATGTGCGCAAGCGCCAGTC‑AAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAACTTACTGAAAACTATTATGCGCGCCGGAACCGCGAc                                                 <  8:51294/151‑1 (MQ=255)

                                                                                               gtgACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACcgcg                                                   <  7:1151664/150‑1 (MQ=255)

                                                                                                   cGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGt                                              <  2:1062862/151‑1 (MQ=255)

                                                                                                   cGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGt                                              <  2:113464/151‑1 (MQ=255)

                                                                                                   cGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGt                                              <  7:262302/151‑1 (MQ=255)

                                                                                                    gCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTg                                             <  1:710106/151‑1 (MQ=255)

                                                                                                    gCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTg                                             >  1:844834/1‑151 (MQ=255)

                                                                                                    gCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTg                                             >  4:975787/1‑151 (MQ=255)

                                                                                                    gCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTg                                             >  7:501129/1‑151 (MQ=255)

                                                                                                    gCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTg                                             <  2:718759/151‑1 (MQ=255)

                                                                                                          aGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAatat                                       >  4:162946/1‑151 (MQ=255)

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

                                                                                                          aGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCCAACGCTACGGCATTGAGCGTAACCGTTCGCCGCGGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAata                                        <  2:34005/150‑1 (MQ=255)

                                                                                                            aTGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATcc                                     <  3:1182307/151‑1 (MQ=255)

                                                                                                                gcgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCtt                                  <  5:419469/150‑1 (MQ=255)

                                                                                                                gcgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTc                                 <  7:146642/151‑1 (MQ=255)

                                                                                                                gcgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTc                                 <  4:1040211/151‑1 (MQ=255)

                                                                                                                gcgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTc                                 <  4:1112217/151‑1 (MQ=255)

                                                                                                                    aaGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTCACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCgcgc                             >  4:220118/1‑151 (MQ=255)

                                                                                                                       cgcCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCAc                          >  8:436939/1‑151 (MQ=255)

                                                                                                                            gTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGTATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTgg                       >  2:581665/1‑149 (MQ=255)

                                                                                                                               aaaaaCCGAGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATCCTGAACTTACTGAAAACTAATATGCGCGCCGCAACCGCGACGGTGAACTATCCCTTCGCGCCACAGGAGGt                   >  1:329639/1‑150 (MQ=255)

                                                                                                                                aaaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGGCa                 <  4:1377410/151‑1 (MQ=255)

                                                                                                                                  aaCCGTGACGTATAAAAAGATCGAACGCTACCGCATTGATCGTAACCGTTCGCCGCTAAAGTAAGGACAGAAGATGCTCAAGTTACTGAAAACTATTATGCGAGACGGAACCGAGACGGTCAAAAATCACTTAGCGCCACTGGATGTCAGc               >  8:283884/1‑151 (MQ=255)

                                                                                                                                    ccGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCTCTGGAGGTCAGccc             <  8:1256795/151‑1 (MQ=255)

                                                                                                                                      gTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAATTATCCCTTCGCGCCACTGGAGGTCAGCCCTg           >  6:89436/1‑151 (MQ=255)

                                                                                                                                      gTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTg           <  3:1304957/151‑1 (MQ=255)

                                                                                                                                             ataAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGGCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTc     <  6:75125/150‑1 (MQ=255)

                                                                                                                                             ataAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTcg    >  2:1343670/1‑151 (MQ=255)

                                                                                                                                             ataAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTc     <  5:114846/150‑1 (MQ=255)

                                                                                                                                              taAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTcgc   >  6:1078235/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: `‑`