BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,392,157	A→T	5.2%	intergenic (‑129/+301)	tldD ← / ← yhdP	putative peptidase/DUF3971‑AsmA2 domains protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,392,157	0	A	T	5.2%	1065.8 / 3.6	367	intergenic (‑129/+301)	tldD/yhdP	putative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand): ref base A (197/136); new base T (2/28); total (199/168)
Fisher's exact test for biased strand distribution p-value = 1.39e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 382 total aligned reads displayed.

ACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTT  >  NC_000913/3392008‑3392306

aCTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                                                                                                                       <  1:4756263/150‑1 (MQ=255)

   tACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATGCGACGCTGTGCGTCTTATCAGGc                                                                                                                                                    <  2:1777736/150‑1 (MQ=37)

   tACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATGCGACGCTGTGCGTCTTATCAGGc                                                                                                                                                    <  2:466552/150‑1 (MQ=37)

     ccAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCaa                                                                                                                                                  >  1:357921/1‑150 (MQ=255)

     ccAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCaa                                                                                                                                                  >  1:2549838/1‑150 (MQ=255)

     ccAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAGCAGGCac                                                                                                                                                  >  1:5307508/1‑149 (MQ=255)

      cAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAAc                                                                                                                                                 >  1:3695519/1‑150 (MQ=255)

      cAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAAc                                                                                                                                                 <  2:5307682/150‑1 (MQ=255)

        ggTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg                                                                                                                                               <  2:3058644/150‑1 (MQ=255)

        ggTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg                                                                                                                                               <  2:2440640/150‑1 (MQ=255)

        ggTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg                                                                                                                                               <  2:1394154/150‑1 (MQ=255)

           tAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTg                                                                                                                                            >  2:5480863/1‑150 (MQ=255)

            aaGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGt                                                                                                                                           >  1:3094127/1‑150 (MQ=255)

              gACTCATGTTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTTTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCg                                                                                                                                         <  2:5085222/150‑1 (MQ=255)

                 tCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACTGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGAt                                                                                                                                      >  2:1927654/1‑150 (MQ=255)

                 tCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAt                                                                                                                                      >  2:842512/1‑150 (MQ=255)

                 tCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTACGGTTTTTAGGCCTGATAAGACGAACAGCGTCGCATCAGGCAACGGCTGACGGAt                                                                                                                                      >  1:565120/1‑150 (MQ=255)

                   aTCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACTCACAGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                                                                                                    >  2:472433/1‑150 (MQ=255)

                      gTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACCGCAGTCGGAAgcggc                                                                                                                                 >  1:4454124/1‑150 (MQ=255)

                         tttGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCAGGGTTGTAGGCCTGATAAGACGAACATCGTCGCATAATGCAACGCCTGACGGAAGAGGAGTa                                                                                                                              >  1:2894888/1‑150 (MQ=37)

                          ttGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTTTAGGCCTGATACAACGAACAGCGAGCCAGCAGGCAACGCCTGACGGATCCGCCGTaa                                                                                                                             >  2:5345093/1‑150 (MQ=21)

                            gCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGACGTAAAc                                                                                                                           >  1:440830/1‑150 (MQ=255)

                             cTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACg                                                                                                                          >  1:3757058/1‑150 (MQ=255)

                                 tCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCtt                                                                                                                      >  1:2136346/1‑150 (MQ=255)

                                  cGTTAGTTACTGCCGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTa                                                                                                                     <  2:2991555/150‑1 (MQ=255)

                                   gttagttaCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGGTAGTAGGCCTGATACGACGCAAAGCGTCGCATCAGGCACCGGCTGGCGGCAGCGGCGAAAACGccggag                                                                                                                    >  1:1551241/1‑146 (MQ=21)

                                    ttatttaCTGAAGTAGAGGTTTTTTTACGGCTGCCGGTTGCGGCGTGAACGCCCTATCCAGCCTACGGATATGTTCCGGTTTGTTGGCCTGATAAGACGCAAATCGTCGCATAAGGCAACGGCTGTCGTATGCGGCGTAAACGCCTTATc                                                                                                                   <  2:3239284/150‑1 (MQ=18)

                                    ttagttaCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATc                                                                                                                   <  2:2669537/150‑1 (MQ=255)

                                       gttaCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGa                                                                                                                >  2:5412063/1‑150 (MQ=255)

                                                gttgtgtaTTTTTTACGGCTGCCGGAGGCGGGGTGAACGCCCTACCCAGCCTACGGTTTTGTTCCGGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                       <  1:3133769/143‑1 (MQ=17)

                                                ggagagGATTTTTTCCGGCTGCCGGAGGCGGCGTGACCGCCCTCTCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                       <  2:372679/148‑1 (MQ=255)

                                                gTAGAGGATTTTTTATGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                       <  2:2269448/150‑1 (MQ=255)

                                                  agagTATTTTTTCCGGCTGCCGGAGGCGGCGTGAACGCCCTACCCAGCCTAGGGTTTTTGTCGGGTTTGTCGGCATTATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt                                                                                                     <  1:999616/150‑1 (MQ=17)

                                                  agagGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt                                                                                                     <  2:4573486/150‑1 (MQ=255)

                                                        ttttttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCg                                                                                               >  1:493260/1‑150 (MQ=255)

                                                        ttttttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCg                                                                                               >  1:1030506/1‑150 (MQ=255)

                                                         tttttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCCGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTGTTCCgg                                                                                              >  2:2827814/1‑149 (MQ=255)

                                                          ttttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGGTTGTAGGCCTGAGAAGACGCACAGCGTCGCAGCAGGCACCGGCTGACCGCTGGGGCGTAAACGCAGAGTCCGACCTACGGGGGTGTTCCggg                                                                                             >  2:3207213/1‑148 (MQ=17)

                                                             tACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGt                                                                                          <  2:3328797/150‑1 (MQ=255)

                                                             tACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCAACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGt                                                                                          <  2:3917451/150‑1 (MQ=255)

                                                              aCGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTa                                                                                         <  2:3614048/150‑1 (MQ=255)

                                                              aCGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTa                                                                                         <  2:899952/150‑1 (MQ=255)

                                                                 gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                      <  2:647087/150‑1 (MQ=255)

                                                                 gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                      <  2:3808181/150‑1 (MQ=255)

                                                                  cTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGcc                                                                                     >  1:428120/1‑150 (MQ=255)

                                                                   tGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                                                                    <  1:3355174/150‑1 (MQ=255)

                                                                   tGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                                                                    <  2:4002263/150‑1 (MQ=255)

                                                                   tGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGGTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGAGGCGAAAACGCCTTAACCGACCTAAGCTTAGGTTCCGATTGTAGGCCt                                                                                    >  1:267485/1‑150 (MQ=14)

                                                                     ccGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGGTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGa                                                                                  >  2:1647699/1‑150 (MQ=255)

                                                                        gATGCGGCGTGAATGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAATGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATaa                                                                               <  1:2773251/150‑1 (MQ=255)

                                                                         aTGCGGCGTGACCGCCCTATCCAGCCTCCGGTTTTGTTCCGTTTTGTAGGCCTGATAAGACGCACCGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAg                                                                              <  1:2294590/150‑1 (MQ=14)

                                                                             gggtttaaCGCCCTATCCCGACTTCGGTTTTTTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGc                                                                          <  1:1695389/144‑1 (MQ=11)

                                                                             ggcgTGACCGCCCTATCCAGCCTACGGTTATGTTCCGTTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGc                                                                          <  1:2178440/150‑1 (MQ=255)

                                                                             ggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGc                                                                          >  1:4993424/1‑150 (MQ=255)

                                                                              gcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTTGGCCTGATAAGACGca                                                                         >  1:5462980/1‑150 (MQ=255)

                                                                                gTGAACGCCCTATCCAGCCTACGTTTTTGTTCCGGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                       <  1:892469/150‑1 (MQ=255)

                                                                                gTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATGCGACGctgt                                                                       <  2:205629/150‑4 (MQ=255)

                                                                                gTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTGTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGCCGTAAACGCCATATCCGACCTACGGTTATGATCCGTTTGTAGGCCTGATAAGACGcaca                                                                       >  1:1900371/1‑150 (MQ=255)

                                                                                   aaCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCg                                                                    >  2:1373740/1‑150 (MQ=255)

                                                                                    aCGCCCTAGCCTGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGGTCCGTTTGTAGGCCTGATAAGACGCACAGCGt                                                                   <  2:1547868/150‑1 (MQ=255)

                                                                                      gCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCg                                                                 <  2:1333633/150‑1 (MQ=255)

                                                                                       cccTATCCAGCCTCCGGTTATGTTCCGGTTTGTCGGCCGGATAAGACGCACCGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGc                                                                <  2:4140973/150‑1 (MQ=14)

                                                                                       cccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGTTTGTAGGCCTGATAAGACGCACAGCGTCGc                                                                >  2:4646127/1‑150 (MQ=255)

                                                                                        ccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                                                               >  1:2957835/1‑150 (MQ=255)

                                                                                           aTCCAGCCTAGGGTTTTGTTCCGGTTTGTTGGCCTTATAAGACGCCCAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                            <  1:2416109/150‑1 (MQ=14)

                                                                                            tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                           >  1:758067/1‑150 (MQ=40)

                                                                                            tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                           <  1:4870274/150‑1 (MQ=40)

                                                                                                   tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg                                                    >  2:479333/1‑150 (MQ=37)

                                                                                                    aCGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGc                                                   >  1:1443630/1‑150 (MQ=37)

                                                                                                       tttttGTTCGGGTTTGTTGGCCTTATAAGACGCACTGGTTCGAATATGGAAACGGCTGTCGGTTGGGGGTTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGc                                                <  2:4405784/146‑1 (MQ=11)

                                                                                                       gTTATGTTCCGGTTTTTAGGCCTGATAAGACTCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGc                                                <  2:3302295/150‑1 (MQ=255)

                                                                                                        ttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGcc                                               >  1:536261/1‑150 (MQ=37)

                                                                                                          aTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCgg                                             <  2:450485/150‑1 (MQ=37)

                                                                                                             ttCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTGCGTTTGTTGGCGTGTTAATACTCACAGCGTCGAATCAGGCAGCGGCTGCCGGATg                                          >  1:3491156/1‑150 (MQ=11)

                                                                                                             ttCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCTTCAGGCAACGGCTGCCGGATg                                          >  1:3928379/1‑150 (MQ=25)

                                                                                                               ccGGTTTGTAGGCCTGATAAGACGCACCGCGGCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcg                                        <  2:111317/150‑1 (MQ=21)

                                                                                                               ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcg                                        >  1:3820014/1‑150 (MQ=34)

                                                                                                                cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCaa                                                       <  2:2895406/134‑1 (MQ=33)

                                                                                                                cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGCCGGATgcgg                                       >  1:2889377/1‑150 (MQ=18)

                                                                                                                cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcgg                                       >  2:3477323/1‑150 (MQ=33)

                                                                                                                cGGTTGGTCGGCCTGATACGACGCCCCGCGGCGCATCAGGCACCGGCTGGCGGATGCGGCGTAACCGCCTTAGCCGACCTCCGGTTAGGTTCCGTTTGTCGGCCTGATAAGACGCCCAGCGTCGCATCAGGCAACGGCTGCCGGATgcgg                                       <  2:4293636/150‑1 (MQ=11)

                                                                                                                 ggTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcggc                                      >  1:1289828/1‑150 (MQ=32)

                                                                                                                   tttGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGt                                    >  2:2167485/1‑150 (MQ=32)

                                                                                                                         ggCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGc                              >  1:3311481/1‑150 (MQ=34)

                                                                                                                          gCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGcc                             >  1:1204098/1‑150 (MQ=34)

                                                                                                                            cTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCtt                           >  1:4606233/1‑150 (MQ=35)

                                                                                                                              gATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCCTAt                         >  1:4061045/1‑150 (MQ=17)

                                                                                                                               aTAAGACGCACCGCGTCGCATCTGGCATCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTAGGGTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATc                        <  2:5372408/150‑1 (MQ=11)

                                                                                                                                  cgccGCACCGCGTCGCACCTGGAAACGGCTGTCGGATGCGGCGTAACCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGa                     <  2:2508480/147‑1 (MQ=11)

                                                                                                                                  aGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCCTATCCGa                     >  2:2324756/1‑150 (MQ=17)

                                                                                                                                    aCGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAcc                   >  1:3860142/1‑150 (MQ=37)

                                                                                                                                    aCGCACAGCGTCGCATCACGCAACGGCTGTCGGATCCGGCGTAAACGCCTTATCCGACATACGGATATCTTCCGTTTGTAGGCCTGATAAGTCGCACAGCGTCGCATCAGTCAACGGCTGCCGGATGCTGCGTAAACTCCTTATCCGAcc                   >  2:1859386/1‑150 (MQ=11)

                                                                                                                                           gCGTCGCATAAGTCAACGGCTGTCGTATGCGGCGTAAACGCCTTATCCGACCTAGGGTTATGTTCCGTTTGTTGGCCTGATAAGACGCACAGCTTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt            <  2:125659/150‑1 (MQ=11)

                                                                                                                                            cGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTa           <  2:4947790/150‑1 (MQ=37)

                                                                                                                                              tCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATg         >  1:4509637/1‑150 (MQ=37)

                                                                                                                                                 cATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGGTTGTAGGCCTGATAAGACGCACAGCGTCGCAGCAGGCAACGGCTGCCGGATGCGGCGGAAACGCCTTATCCGACCTACGGTTATGTTc      >  1:5236485/1‑150 (MQ=255)

                                                                                                                                                  ctCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTGCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTcc     <  2:3372133/149‑1 (MQ=38)

                                                                                                                                                     aGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGtt  >  1:4252983/1‑150 (MQ=40)

                                                                                                                                                     aGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGtt  >  1:2700732/1‑150 (MQ=40)

                                                                                                                                                     aGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGtt  >  1:1762327/1‑150 (MQ=40)

ACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTT  >  NC_000913/3392008‑3392306

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