BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	805,312	G→T	6.8%	E604* (GAG→TAG) ‡	ybhJ →	putative hydratase YbhJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	805,312	0	G	T	6.8%	349.5 / 6.9	105	E604* (GAG→TAG) ‡	ybhJ	putative hydratase YbhJ
Reads supporting (aligned to +/- strand): ref base G (48/49); new base T (5/3); total (53/52)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGA  >  NC_000913/805165‑805459

cTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAg                                                                                                                                                   >  1:5664/1‑150 (MQ=255)

  gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCTGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCt                                                                                                                                                 >  1:1374012/1‑150 (MQ=255)

  gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCt                                                                                                                                                 >  2:136058/1‑150 (MQ=255)

    ttccttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGa                                                                                                                                               <  1:968093/150‑1 (MQ=255)

      ccttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACa                                                                                                                                             >  1:299006/1‑150 (MQ=255)

      ccttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACa                                                                                                                                             >  2:1119462/1‑150 (MQ=255)

       cttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAAATGTCAGCGAGCTGACa                                                                                                                                             >  2:194596/1‑150 (MQ=255)

       cttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACag                                                                                                                                            >  1:718605/1‑150 (MQ=255)

        ttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACaga                                                                                                                                           >  1:1283112/1‑150 (MQ=255)

         tccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCTCGCTGACattc                                                                                                                                          <  1:483039/150‑4 (MQ=255)

            ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCTCGCTGACattcccc                                                                                                                                       <  1:422753/150‑7 (MQ=255)

            ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCTCGCTGACattcccc                                                                                                                                       <  1:1209040/150‑7 (MQ=255)

            ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGgtg                                                                                                                                       >  1:1140826/1‑150 (MQ=255)

            ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGgtg                                                                                                                                       >  1:1109613/1‑150 (MQ=255)

               gAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGttt                                                                                                                                    <  1:90424/150‑1 (MQ=255)

                      cttcttATCGTTCAAATCCAATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                                                                                                                             >  1:507370/1‑150 (MQ=255)

                         cttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTa                                                                                                                          >  2:1435831/1‑150 (MQ=255)

                         cttATCGTTCAAATCCGATTGGGCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTa                                                                                                                          <  1:476473/150‑1 (MQ=255)

                          ttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTaa                                                                                                                         >  1:518392/1‑150 (MQ=255)

                           tATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAg                                                                                                                        <  1:1080661/150‑1 (MQ=255)

                              cGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAg                                                                                                                     >  1:1297434/1‑150 (MQ=255)

                              cGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAg                                                                                                                     >  1:601042/1‑150 (MQ=255)

                               gTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGTGCTGACAGAGGTGTTTGCGCGCATTAAGCAGa                                                                                                                    >  1:382100/1‑150 (MQ=255)

                                  cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTg                                                                                                                 <  1:882421/150‑1 (MQ=255)

                                  cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGAATGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTg                                                                                                                 >  2:1065317/1‑150 (MQ=255)

                                   aaaTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGc                                                                                                                <  2:718605/150‑1 (MQ=255)

                                   aaaTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGc                                                                                                                >  1:2632/1‑150 (MQ=255)

                                    aaTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGa                                                                                                                                               >  2:447505/1‑118 (MQ=255)

                                    aaTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGa                                                                                                                                               <  1:447505/118‑1 (MQ=255)

                                     aTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTATCAGAAGTAAAGCGACTGGTTATCTGGAAAATCATCGTCTGGCGGGT‑AATGTCATAGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTg                                                                                                              <  1:91416/150‑1 (MQ=255)

                                     aTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTg                                                                                                              <  1:136058/150‑1 (MQ=255)

                                      tCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTgg                                                                                                             <  1:1249559/150‑1 (MQ=255)

                                        cGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGTGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTc                                                                                                           <  2:382100/150‑1 (MQ=255)

                                         gATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                                                                          <  1:206576/150‑1 (MQ=255)

                                         gATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                                                                          <  2:1504648/150‑1 (MQ=255)

                                             ggTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                      <  2:1261048/150‑1 (MQ=255)

                                             ggTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                      <  1:1119462/150‑1 (MQ=255)

                                               tCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa                                                                                                    <  2:1486491/150‑1 (MQ=255)

                                                   gCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATAtt                                                                                                <  2:174541/150‑1 (MQ=255)

                                                     ggAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGa                                                                                              <  2:148107/150‑1 (MQ=255)

                                                       aGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                                                                            >  1:1107657/1‑150 (MQ=255)

                                                        gTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                                                                           >  2:1535723/1‑150 (MQ=255)

                                                            aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTg                                                                                       <  1:916946/150‑1 (MQ=255)

                                                            aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGGTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTg                                                                                       <  1:663085/150‑1 (MQ=255)

                                                             cccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGc                                                                                      >  2:67168/1‑150 (MQ=255)

                                                             cccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGc                                                                                      >  1:6402/1‑150 (MQ=255)

                                                                   ctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                      >  1:206431/1‑128 (MQ=255)

                                                                   ctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTg                                                                                >  1:1055069/1‑150 (MQ=255)

                                                                    tcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGa                                                                               >  2:220800/1‑150 (MQ=255)

                                                                                cGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                                                                   <  2:505136/150‑1 (MQ=255)

                                                                                  gTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATgg                                                                 <  1:1504242/150‑1 (MQ=255)

                                                                                   ttATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGt                                                                <  1:1505737/150‑1 (MQ=255)

                                                                                    tATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                               <  2:67211/150‑1 (MQ=255)

                                                                                      tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                                                                  <  1:518843/129‑1 (MQ=255)

                                                                                      tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                                                                  >  2:518843/1‑129 (MQ=255)

                                                                                      tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa                                                             <  2:838411/150‑1 (MQ=255)

                                                                                      tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa                                                             >  2:674899/1‑150 (MQ=255)

                                                                                      tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa                                                             >  1:741977/1‑150 (MQ=255)

                                                                                         tAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc                                                          >  1:823624/1‑150 (MQ=255)

                                                                                           gCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCgg                                                        >  1:1016086/1‑150 (MQ=255)

                                                                                           gCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCgg                                                        <  2:577390/150‑1 (MQ=255)

                                                                                           gCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGGCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCgg                                                        <  1:1355025/150‑1 (MQ=255)

                                                                                             aGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTg                                                                          >  2:148574/1‑130 (MQ=255)

                                                                                             aGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTg                                                                          <  1:148574/130‑1 (MQ=255)

                                                                                             aGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg                                                      <  2:113058/150‑1 (MQ=255)

                                                                                             aGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg                                                      >  1:113058/1‑150 (MQ=255)

                                                                                             aGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg                                                      <  1:159965/150‑1 (MQ=255)

                                                                                              gAAGTATAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGTCAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGa                                                     <  1:1506946/150‑1 (MQ=255)

                                                                                                aGTAAAGCGCCTCCTGAGCTGGCTCATCAGCGTCTGGCGGGG‑AATGTCACCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCGGCCAACTGAAATTGGCAGAATGGTCTAGGCGGTGacc                                                   >  2:1067245/1‑148 (MQ=255)

                                                                                                 gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAc                                                  <  2:1140826/150‑1 (MQ=255)

                                                                                                 gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAc                                                  <  2:126507/150‑1 (MQ=255)

                                                                                                  tAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAcc                                                 <  2:601042/150‑1 (MQ=255)

                                                                                                   aaaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                <  1:1435831/150‑1 (MQ=255)

                                                                                                      gCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGc                                             <  2:350149/150‑1 (MQ=255)

                                                                                                      gCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGc                                             <  1:643399/150‑1 (MQ=255)

                                                                                                       cGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCg                                            <  1:1056306/150‑1 (MQ=255)

                                                                                                       cGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCg                                            >  2:890402/1‑150 (MQ=255)

                                                                                                         aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGAt                                          <  2:518392/150‑1 (MQ=255)

                                                                                                           tgctgAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATgg                                        >  1:777063/1‑150 (MQ=255)

                                                                                                            gctgAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGt                                       <  2:1393105/150‑1 (MQ=255)

                                                                                                            gctgAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGt                                       >  1:755160/1‑150 (MQ=255)

                                                                                                             ctgAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGtt                                      <  1:746007/150‑1 (MQ=255)

                                                                                                                 gCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGAAGGTTCTgc                                  >  1:407686/1‑150 (MQ=255)

                                                                                                                          tCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACa                         >  1:104503/1‑150 (MQ=255)

                                                                                                                           cAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAg                        >  1:935227/1‑150 (MQ=255)

                                                                                                                                 cTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGgcggcg                  <  2:823624/150‑1 (MQ=255)

                                                                                                                                   ggCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAg                >  2:1461316/1‑150 (MQ=255)

                                                                                                                                      gggg‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg             <  1:930132/150‑1 (MQ=255)

                                                                                                                                       cacctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg             >  1:81000/7‑150 (MQ=255)

                                                                                                                                       cacctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg             >  2:1234185/7‑150 (MQ=255)

                                                                                                                                        acctctGTCAGCTAGCTGACAGAGGTGATTGCGCGCATTAAGCAGATTGCTGGTCAGCAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGc            >  2:37940/6‑150 (MQ=255)

                                                                                                                                         cctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTGCTGCGCGTGAACAGGCGGCGAGCTGcc           >  1:607333/5‑150 (MQ=255)

                                                                                                                                            aTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGc        >  1:704821/1‑150 (MQ=255)

                                                                                                                                             tGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg       <  2:1107657/150‑1 (MQ=255)

                                                                                                                                             tGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg       <  2:5664/150‑1 (MQ=255)

                                                                                                                                             tGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg       <  1:1065317/150‑1 (MQ=255)

                                                                                                                                              gTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt      >  2:538284/1‑150 (MQ=255)

                                                                                                                                                cAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtgt    >  2:680460/1‑150 (MQ=255)

                                                                                                                                                 ggCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGtgtg   >  1:379015/2‑150 (MQ=255)

                                                                                                                                                  gCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGa  >  1:1384909/1‑150 (MQ=255)

CTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGG‑AATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGA  >  NC_000913/805165‑805459

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