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.3%	1528.0 / ‑1.1	405	E541K (GAG→AAG)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand): ref base G (1/2); new base A (173/228); total (174/231)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 411 total aligned reads displayed.

GGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTCGCG  >  NC_000913/2609961‑2610254

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

ggCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGAt                                                                                                                                                 <  7:1202930/151‑1 (MQ=255)

ggCCGGGGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAaaag                                                                                                                                                   <  6:215302/149‑1 (MQ=255)

  ccGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCg                                                                                                                                               <  4:551782/151‑1 (MQ=255)

    ggTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGaa                                                                                                                                             <  7:547962/151‑1 (MQ=255)

     gtgtTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAAc                                                                                                                                            <  1:1180492/151‑1 (MQ=255)

      tgtTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACg                                                                                                                                           >  6:1629674/1‑151 (MQ=255)

       gtTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACg                                                                                                                                           <  3:729455/150‑1 (MQ=255)

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

           cGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACg                                                                                                                                      >  3:1769286/1‑151 (MQ=255)

               aCATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCAt                                                                                                                                  <  4:1137376/151‑1 (MQ=255)

               aCATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGGACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCAt                                                                                                                                  <  4:415180/151‑1 (MQ=255)

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

                 aTGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGCCCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTg                                                                                                                                >  1:636442/1‑151 (MQ=255)

                        gcgGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTa                                                                                                                          >  7:870018/1‑150 (MQ=255)

                        gcgGCAATACCGTTTCTGACGCACAGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTa                                                                                                                          >  8:1408226/1‑150 (MQ=255)

                          ggCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAcc                                                                                                                       >  5:1720941/1‑151 (MQ=255)

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

                              aTACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTc                                                                                                                   >  4:635848/1‑151 (MQ=255)

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

                                 ccGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCCAACGCTACGGCATTGATCGTAACCGTTcgcc                                                                                                                >  5:766681/1‑151 (MQ=255)

                                 ccGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCGGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcgcc                                                                                                                <  5:1864247/151‑1 (MQ=255)

                                  cGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcgccg                                                                                                               >  2:748751/1‑151 (MQ=255)

                                  cGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcgccg                                                                                                               >  2:1056456/1‑151 (MQ=255)

                                    tttCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCt                                                                                                             >  7:1073190/1‑151 (MQ=255)

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

                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGa                                                                                                           <  7:743850/148‑1 (MQ=255)

                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGt                                                                                                        <  7:1243/151‑1 (MQ=255)

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

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

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

                                             cACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAgg                                                                                                    <  7:1517691/151‑1 (MQ=255)

                                               ccGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGAc                                                                                                  <  2:494258/151‑1 (MQ=255)

                                                 gCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGAc                                                                                                  >  6:323383/1‑149 (MQ=255)

                                                    gATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACagaag                                                                                             <  2:1549972/151‑1 (MQ=255)

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

                                                       tATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATGAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGAt                                                                                           <  2:837766/150‑1 (MQ=255)

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

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

                                                           ggTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGa                                                                                      <  6:1396917/151‑1 (MQ=255)

                                                           ggTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGa                                                                                      <  3:1673040/151‑1 (MQ=255)

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

                                                            gTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGaa                                                                                     <  3:1293428/151‑1 (MQ=255)

                                                              aGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAg                                                                                    <  5:285028/150‑1 (MQ=255)

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

                                                                   tGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACt                                                                               <  6:520854/150‑1 (MQ=255)

                                                                    gATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTg                                                                              >  5:1473662/1‑150 (MQ=255)

                                                                    gATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGa                                                                             <  3:528393/151‑1 (MQ=255)

                                                                     aTCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGaa                                                                            <  3:1554261/151‑1 (MQ=255)

                                                                       cccTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGaaaa                                                                          <  8:1044931/151‑1 (MQ=255)

                                                                        ccTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAAc                                                                         >  4:32729/1‑151 (MQ=255)

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

                                                                        ccTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAAc                                                                         <  3:934095/151‑1 (MQ=255)

                                                                         cTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACt                                                                        >  7:65220/1‑151 (MQ=255)

                                                                          tGCTACTCCTGTACCGCCCGTGTGACGCTGGTACATGTGCGCCAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAACTTACTGAAAACt                                                                        >  4:1386688/1‑150 (MQ=255)

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

                                                                           gCTAGTCCTGTACCGACGGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGGTGAAGTTACTGAAAACta                                                                       <  2:1759329/150‑1 (MQ=255)

                                                                           gCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACta                                                                       <  2:986162/150‑1 (MQ=255)

                                                                            cTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtatt                                                                     >  3:1235880/1‑151 (MQ=255)

                                                                            cTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtatt                                                                     >  3:633539/1‑151 (MQ=255)

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

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

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

                                                                                tCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgc                                                                 >  8:327164/1‑151 (MQ=255)

                                                                                 ccTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcg                                                                <  7:1729436/151‑1 (MQ=255)

                                                                                 ccTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgc                                                                 <  7:425263/150‑1 (MQ=255)

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

                                                                                    gTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGATCGCTACGGGATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcgc                                                               <  2:285133/149‑1 (MQ=255)

                                                                                    gTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcgc                                                               <  1:1325628/149‑1 (MQ=255)

                                                                                         gaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAGCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGATGTTACTGAAAACTTTTATGCGCGCCGGac                                                        >  8:1195549/1‑150 (MQ=255)

                                                                                         gaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGGTCGCCGCTGAAGTAAGGACAGAAGATGATGAAGATACTGAAAACTATTATTCGCGCCGGa                                                         >  8:1704746/1‑150 (MQ=255)

                                                                                           ccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAAcc                                                      >  4:1841842/1‑151 (MQ=255)

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

                                                                                            cgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAGCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAGCTATTATGCGCGCCGGACCcg                                                     >  2:1197154/1‑151 (MQ=255)

                                                                                            cgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACcg                                                     >  2:198506/1‑151 (MQ=255)

                                                                                             gtgtGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGACCGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACcgc                                                    >  2:137596/1‑151 (MQ=255)

                                                                                             gtgtGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACcgc                                                    >  6:1726779/1‑151 (MQ=255)

                                                                                              tgtgACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACcgcg                                                   >  5:1087385/1‑151 (MQ=255)

                                                                                                  aCGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACgg                                               >  1:1055648/1‑151 (MQ=255)

                                                                                                   cGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACggg                                              >  2:588908/1‑150 (MQ=255)

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

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

                                                                                                     cTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGCAAACTATTATGCGCGCCGGAACCGCGACGGTgc                                            >  6:1249514/1‑150 (MQ=255)

                                                                                                       ggTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGaaa                                          <  3:1071103/151‑1 (MQ=255)

                                                                                                        gTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGATGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAat                                         <  3:438222/151‑1 (MQ=255)

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

                                                                                                              gtgCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCt                                   >  4:1522430/1‑151 (MQ=255)

                                                                                                               tgCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGCTCGAACGCTACGGCATTGATCGTACCCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAACACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCtt                                  >  5:586903/1‑151 (MQ=255)

                                                                                                               tgCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGCACCGCGACGGTGAAATATCCCtt                                  >  8:1653912/1‑151 (MQ=255)

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

                                                                                                                 cgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTcg                                <  3:27046/151‑1 (MQ=255)

                                                                                                                    aaGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCgcgc                             >  6:116798/1‑151 (MQ=255)

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

                                                                                                                           aGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGa                      <  4:737816/151‑1 (MQ=255)

                                                                                                                              cAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGt                   <  3:779629/151‑1 (MQ=255)

                                                                                                                              aaaaaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAATTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGt                   <  5:1669533/150‑1 (MQ=255)

                                                                                                                                 aaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTc                  >  8:1859335/1‑149 (MQ=255)

                                                                                                                                        gCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGc         <  8:549957/151‑1 (MQ=255)

                                                                                                                                         ccGTATAAAAAGATCGAACGCTACGGCATTGGTCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCt        <  8:1552053/151‑1 (MQ=255)

                                                                                                                                             ataAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTcg    <  3:321875/151‑1 (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: `‑`