BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,375	+A	coding (458/516 nt)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,372	1	.	A	99.2%	481.4 / ‑3.8	132	coding (455/516 nt)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base A (35/96); total (36/96)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

GCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGC‑‑‑G‑‑‑‑AAC‑T‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCC  >  NC_000913/4008280‑4008449

gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTc                                                                                              <  1:4755172/93‑1 (MQ=255)

gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTc                                                                                              >  1:3964422/1‑93 (MQ=255)

gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTc                                                                                              <  1:3045384/93‑1 (MQ=255)

gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTc                                                                                              >  2:4412301/1‑93 (MQ=255)

gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTc                                                                                              <  1:5617646/93‑1 (MQ=255)

 cTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGTTTC‑a                                                                                            >  1:3694321/1‑93 (MQ=255)

  tGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aa                                                                                           >  2:6985593/1‑93 (MQ=255)

    gaccAGGCGATTAACGCTCGTGGGTGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaa                                                                                         <  1:53737/90‑2 (MQ=255)

    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaa                                                                                         <  2:1091025/93‑2 (MQ=255)

    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaa                                                                                         <  1:3323321/93‑2 (MQ=255)

    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaa                                                                                         <  2:4294971/93‑2 (MQ=255)

    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTC‑aaaa                                                                                         >  1:5685688/1‑92 (MQ=255)

     aGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGTTTC‑aaaag                                                                                        >  1:3471991/1‑91 (MQ=255)

     aGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaag                                                                                        <  1:5196705/93‑3 (MQ=255)

     aGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaag                                                                                        <  2:464164/93‑3 (MQ=255)

     aGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGATTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaag                                                                                        >  1:3256397/1‑91 (MQ=255)

      gCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGTTTC‑aaaaag                                                                                       >  2:5254769/1‑90 (MQ=255)

      gCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaagg                                                                                       <  2:3998768/93‑4 (MQ=255)

       cAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaaggc                                                                                      >  1:1483609/1‑89 (MQ=255)

       cAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑aaaaggc                                                                                      <  2:1641142/93‑5 (MQ=255)

         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGTTTCAAAAGGCCa                                                                                     >  2:6456540/1‑93 (MQ=255)

         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCa                                                                                     >  2:6162742/1‑93 (MQ=255)

         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCa                                                                                     <  1:2973541/93‑1 (MQ=255)

         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCa                                                                                     <  2:2277139/93‑1 (MQ=255)

         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCa                                                                                     <  1:5038271/93‑1 (MQ=255)

          gCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAt                                                                                    <  1:5231259/93‑1 (MQ=255)

           cGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATg                                                                                   >  1:5129671/1‑93 (MQ=255)

            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgc                                                                                  <  1:2662152/93‑1 (MQ=255)

            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgc                                                                                  <  2:270367/93‑1 (MQ=255)

            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgc                                                                                  <  1:1176917/93‑1 (MQ=255)

             aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATgcg                                                                                 >  2:7102674/1‑93 (MQ=255)

             aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATgag                                                                                 >  2:3594123/1‑91 (MQ=255)

              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgcgc                                                                                <  2:3022571/93‑1 (MQ=255)

              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgcgc                                                                                <  2:6206847/93‑1 (MQ=255)

              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgcgc                                                                                <  1:1756611/93‑1 (MQ=255)

              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgcgc                                                                                <  1:821666/93‑1 (MQ=255)

              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgcgc                                                                                <  1:7789852/93‑1 (MQ=255)

               tAACGCTCGTTGGGGCTATCGCGTAAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCt                                                                               <  2:415083/93‑1 (MQ=255)

               tAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCt                                                                               <  2:7935202/93‑1 (MQ=255)

                aaCGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCCTTTCAAAAGGCCATtcgccg                                                                              >  1:1632266/1‑91 (MQ=255)

                 aCGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTgg                                                                             >  2:5359641/1‑93 (MQ=255)

                    cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAAGCCATGCGCTGGAAg                                                                          >  2:2825328/1‑93 (MQ=255)

                     tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGGAgc                                                                         >  2:6624528/1‑92 (MQ=255)

                     tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                         <  2:6408742/93‑2 (MQ=255)

                     tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                         <  1:590139/93‑2 (MQ=255)

                      cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcg                                                                        <  1:6264619/93‑3 (MQ=255)

                      cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcg                                                                        <  1:3847176/93‑3 (MQ=255)

                      cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcg                                                                        <  2:7316603/93‑3 (MQ=255)

                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                       <  2:3276488/93‑4 (MQ=255)

                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                       <  1:1678421/93‑4 (MQ=255)

                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAgcgg                                                                       >  1:825140/1‑90 (MQ=255)

                        ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                      <  2:5206941/93‑5 (MQ=255)

                        ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                      <  1:2023289/93‑5 (MQ=255)

                        ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                      <  2:6272242/93‑5 (MQ=255)

                        ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                      <  1:7582962/93‑5 (MQ=255)

                        ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                      <  1:6940488/93‑5 (MQ=255)

                        ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggg                                                                      >  1:6786195/1‑89 (MQ=255)

                              cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aa                                                         <  2:5876197/93‑1 (MQ=255)

                              cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aa                                                         <  2:3476810/93‑1 (MQ=255)

                              cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aa                                                         >  2:7825112/1‑93 (MQ=255)

                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGGAGCGGTGGC‑‑‑G‑‑‑‑aaga                                                       >  1:3816411/1‑91 (MQ=255)

                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aaga                                                       <  2:1742079/93‑3 (MQ=255)

                                  cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aagaaa                                                     <  1:6133066/93‑5 (MQ=255)

                                  cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aagaaa                                                     <  1:7277068/93‑5 (MQ=255)

                                  cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aagaaa                                                     <  2:7537857/93‑5 (MQ=255)

                                  cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aagaaa                                                     <  2:5096779/93‑5 (MQ=255)

                                                                  gcTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGtt                        <  1:2469917/92‑1 (MQ=255)

                                                                   cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGttt                       <  1:1059237/93‑2 (MQ=255)

                                                                    tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGtttg                      <  1:6350927/93‑3 (MQ=255)

                                                                     ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                     >  2:137114/1‑90 (MQ=255)

                                                                     ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                     <  1:251791/93‑4 (MQ=255)

                                                                     ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                     >  2:1261370/1‑90 (MQ=255)

                                                                      gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGtttggg                    >  2:5683019/1‑89 (MQ=255)

                                                                      gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGtttggg                    >  1:2317758/1‑89 (MQ=255)

                                                                            aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg              <  2:3772951/93‑1 (MQ=255)

                                                                            aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg              >  2:502144/1‑93 (MQ=255)

                                                                             tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa             <  2:4595843/93‑1 (MQ=255)

                                                                             tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa             <  1:3150267/93‑1 (MQ=255)

                                                                             tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa             <  1:7051633/93‑1 (MQ=255)

                                                                              ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt            <  2:5824838/93‑1 (MQ=255)

                                                                              ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt            <  2:4703577/93‑1 (MQ=255)

                                                                               gcgggcggGGTTTCAAAACGCCATGCGCTGGAAGCGGTGGCTAAGA‑‑‑AAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg           <  2:2141291/93‑1 (MQ=255)

                                                                                cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa          <  1:2852182/93‑1 (MQ=255)

                                                                                 gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa         <  1:6448749/93‑1 (MQ=255)

                                                                                 gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa         <  1:297305/93‑1 (MQ=255)

                                                                                 gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa         <  1:6162743/93‑1 (MQ=255)

                                                                                  ggcggcGTTTCAAAAGGCCATGCGGTGGAAGCGGTGGC‑‑‑GAAGAAACGG‑GGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac        <  1:7078414/93‑1 (MQ=255)

                                                                                  ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac        <  1:2454381/93‑1 (MQ=255)

                                                                                  ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac        <  1:7078413/93‑1 (MQ=255)

                                                                                  ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac        >  1:2114222/1‑93 (MQ=255)

                                                                                  ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAGGAac        >  2:2476240/1‑93 (MQ=255)

                                                                                   gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg       >  2:2387781/1‑93 (MQ=255)

                                                                                    cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga      <  2:184452/93‑1 (MQ=255)

                                                                                    cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga      >  1:4484584/1‑93 (MQ=255)

                                                                                      gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacc    >  1:5719364/1‑92 (MQ=255)

                                                                                       cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc   >  2:3399617/1‑93 (MQ=255)

                                                                                       cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc   <  2:5618425/93‑1 (MQ=255)

                                                                                        gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcg  >  1:1056564/1‑92 (MQ=255)

                                                                                        gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc  <  2:6601284/93‑1 (MQ=255)

                                                                                        gTTTCAAAAGGCCATGCGCTGGAAGCGGCGGC‑‑‑GAAGAAAC‑TGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcg  >  2:5160899/1‑92 (MQ=255)

GCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGC‑‑‑G‑‑‑‑AAC‑T‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCC  >  NC_000913/4008280‑4008449

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`