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	98.0%	1495.3 / 9.0	394	E541K (GAG→AAG)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand): ref base G (3/5); new base A (188/197); total (191/203)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 401 total aligned reads displayed.

GGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTCGCG  >  NC_000913/2609961‑2610254

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

  ccGGTGTTGCGTTACATGCTCCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATCTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATcc                                                                                                                                               >  7:278219/1‑150 (MQ=255)

   cGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCg                                                                                                                                               >  6:1027901/1‑150 (MQ=255)

    ggTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGaa                                                                                                                                             <  8:1578579/151‑1 (MQ=255)

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

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

          gCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGTCCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTAc                                                                                                                                       <  4:1714749/151‑1 (MQ=255)

          gCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTa                                                                                                                                        <  7:413107/150‑1 (MQ=255)

           ggTTAGATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGGGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACg                                                                                                                                      <  2:1179699/150‑1 (MQ=255)

           cGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACg                                                                                                                                      <  5:1279065/151‑1 (MQ=255)

           cGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACg                                                                                                                                      <  6:1220002/151‑1 (MQ=255)

            gTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACgg                                                                                                                                     >  7:918714/1‑151 (MQ=255)

             ttACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGc                                                                                                                                    >  2:864861/1‑151 (MQ=255)

                   gctgcGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGa                                                                                                                               <  8:220541/150‑1 (MQ=255)

                         cgGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAc                                                                                                                        <  3:333367/151‑1 (MQ=255)

                          ggCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAcc                                                                                                                       >  2:1376878/1‑151 (MQ=255)

                          ggCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAcc                                                                                                                       <  1:995935/151‑1 (MQ=255)

                           gCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCg                                                                                                                      >  1:1334677/1‑151 (MQ=255)

                            cAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTTATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTTCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGt                                                                                                                     <  5:373711/151‑1 (MQ=255)

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

                             aaTACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGtt                                                                                                                    <  2:1557213/151‑1 (MQ=255)

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

                               tACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcg                                                                                                                  >  8:294922/1‑151 (MQ=255)

                               tACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcg                                                                                                                  >  3:903433/1‑151 (MQ=255)

                                  cGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcgccg                                                                                                               >  6:292545/1‑151 (MQ=255)

                                     ttCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCt                                                                                                             <  7:1293188/150‑1 (MQ=255)

                                      tCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTg                                                                                                            <  4:1553945/150‑1 (MQ=255)

                                        tGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGa                                                                                                           >  5:754831/1‑149 (MQ=255)

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

                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGt                                                                                                        <  3:697150/151‑1 (MQ=255)

                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGACCGTAACCGTTCGCCGCTGa                                                                                                           <  8:71035/148‑1 (MQ=255)

                                          aagCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTTTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCTTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTa                                                                                                       <  1:1007432/149‑1 (MQ=255)

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

                                             cACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAgg                                                                                                    <  8:1364362/151‑1 (MQ=255)

                                             cACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCCTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTAGCCGCTGAAGTAAgg                                                                                                    >  3:1419903/1‑151 (MQ=255)

                                              aCCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAgg                                                                                                    >  8:358014/1‑150 (MQ=255)

                                               ccGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGAc                                                                                                  <  8:1481437/151‑1 (MQ=255)

                                                 gCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGAc                                                                                                  >  8:116790/1‑149 (MQ=255)

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

                                                    gATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACagaa                                                                                              <  7:510775/150‑1 (MQ=255)

                                                    gATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACaga                                                                                               <  3:980989/149‑1 (MQ=255)

                                                    gATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACaga                                                                                               <  2:290119/149‑1 (MQ=255)

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

                                                       tATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATg                                                                                          >  2:1315085/1‑151 (MQ=255)

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

                                                          cGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTg                                                                                       <  5:703279/151‑1 (MQ=255)

                                                                 gttGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTAc                                                                                <  4:624968/150‑1 (MQ=255)

                                                                 cTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTAc                                                                                >  3:1133725/1‑151 (MQ=255)

                                                                 cTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTAc                                                                                >  3:1532027/1‑151 (MQ=255)

                                                                  ttGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACt                                                                               >  3:1360400/1‑151 (MQ=255)

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

                                                                  ttGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCCTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACt                                                                               >  7:1374589/1‑151 (MQ=255)

                                                                   tGATCCCTGCTACTCCTGTACCGCCCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTAGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTg                                                                              >  3:465582/1‑151 (MQ=255)

                                                                    gATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTg                                                                              <  4:1389139/150‑1 (MQ=255)

                                                                    gATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTg                                                                              <  2:1748952/150‑1 (MQ=255)

                                                                       cccTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGaaa                                                                           <  6:1117134/150‑1 (MQ=255)

                                                                       cccTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGaaa                                                                           <  8:1690642/150‑1 (MQ=255)

                                                                        ccTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAAc                                                                         <  2:522478/151‑1 (MQ=255)

                                                                            cTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTCAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtat                                                                      >  1:1149398/1‑150 (MQ=255)

                                                                              aCTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtatta                                                                    <  4:275428/150‑1 (MQ=255)

                                                                              aCTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtatt                                                                     <  8:125063/149‑1 (MQ=255)

                                                                                   tGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcgcg                                                              >  4:45248/1‑151 (MQ=255)

                                                                                   tGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcgcg                                                              >  6:1325105/1‑151 (MQ=255)

                                                                                   tGTACCGACCGTGTGACGCTGGTAGATGTGAGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcgcg                                                              >  4:718478/1‑151 (MQ=255)

                                                                                     tACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGcc                                                            >  5:1222799/1‑151 (MQ=255)

                                                                                        cgaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCgg                                                          >  4:1151361/1‑150 (MQ=255)

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

                                                                                         gaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGaa                                                        <  5:226494/151‑1 (MQ=255)

                                                                                           ccgTGTTACTCTTGGAGATGTGCGCAAGCGCCAGTCAAAAACCTTGAATTATAAAAATATCTAAATATACGGCATTGATCGTAACCTTTATCATCTTAAGTAATTACATAAGATGCTTAAGTTACTGAAAAATATTATGCGATCCGGAAca                                                      <  1:453656/151‑2 (MQ=255)

                                                                                                 gACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTAAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGAc                                                 <  3:1715917/150‑1 (MQ=255)

                                                                                                    gCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCTAACTCTACGGCATTGATCGTAACCGTTCGGCGCTTAAGTTAGGACAGAAGATTCTGAAGTTACTTAAAACTATTATGCGCTACGGAACCTCGAGGGTg                                             <  5:166364/151‑1 (MQ=255)

                                                                                                        gTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAat                                         <  1:902203/151‑1 (MQ=255)

                                                                                                         tAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAat                                         >  7:610350/1‑150 (MQ=255)

                                                                                                          aGATGTGCGCAATCGCCAGTCAAAAACCGTGCCTTATAAAAAGATCTAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTTAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAatat                                       <  7:1556378/151‑1 (MQ=255)

                                                                                                          aGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGatatat                                       >  6:665125/1‑151 (MQ=255)

                                                                                                          aGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAatat                                       >  5:939220/1‑151 (MQ=255)

                                                                                                              gtgCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCCAACCCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGTACAGAAGATGCAGAAGTTACTGAAAACTATTATGCGAGCAGCAACAGCAACAGTGAAATATccc                                    >  6:91224/1‑150 (MQ=255)

                                                                                                               ggcgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCtt                                  <  2:69469/150‑1 (MQ=255)

                                                                                                                gcgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTc                                 >  6:338307/1‑151 (MQ=255)

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

                                                                                                                   cAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTcgcg                              >  7:1144069/1‑151 (MQ=255)

                                                                                                                      gcgcCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCa                           <  3:208967/151‑1 (MQ=255)

                                                                                                                        gcCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGcc                                                            <  7:1671984/116‑1 (MQ=255)

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

                                                                                                                        gcCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACt                         <  8:407824/151‑1 (MQ=255)

                                                                                                                        gcCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACt                         <  5:597637/151‑1 (MQ=255)

                                                                                                                           aGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTgg                       >  4:540825/1‑150 (MQ=255)

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

                                                                                                                               aaaaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTc                  >  8:1650296/1‑151 (MQ=255)

                                                                                                                                aaaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTGGGCGCGGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCGGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCa                 <  4:863173/151‑1 (MQ=255)

                                                                                                                                aaaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCa                 <  2:943427/151‑1 (MQ=255)

                                                                                                                                 aaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTc                  >  3:57211/1‑149 (MQ=255)

                                                                                                                                 aaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAg                <  1:1575957/151‑1 (MQ=255)

                                                                                                                                  aaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGc               >  2:831586/1‑151 (MQ=255)

                                                                                                                                   aCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGcc              <  1:1618737/151‑1 (MQ=255)

                                                                                                                                     cGTGCCGTATAAAAAGATCGCACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGCTGAAATATCCCTTCGCGCCACTGGAGCTCAGCCCt            >  7:1450818/1‑151 (MQ=255)

                                                                                                                                        gCCGTATAAAAAGATGGAACGCTACGTCATTGATCGTAACCGTTCGCCTCTGAAGTAAGGTCAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTgg          <  6:216041/150‑1 (MQ=255)

                                                                                                                                           gTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCttt      >  4:1455864/1‑151 (MQ=255)

                                                                                                                                             ataAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTc     <  6:1505727/150‑1 (MQ=255)

                                                                                                                                               aaaaAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGCACCGCGCCGGTGAAAGATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTcgcg  >  7:1088374/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: `‑`