BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,616,692	C→T	G81S (GGC→AGC)	ydeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,692	0	C	T	98.2%	386.3 / ‑3.5	114	G81S (GGC→AGC)	ydeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (72/40); total (74/40)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 124 total aligned reads displayed.

GCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616602‑1616781

gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         >  1:490347/1‑93 (MQ=255)

gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         <  1:3409479/93‑1 (MQ=255)

gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         >  2:3136776/1‑93 (MQ=255)

gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         >  2:3354222/1‑93 (MQ=255)

 ctctTCCGACTTGCTTTTCGCCTCCGGTGAATCAGTCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgcc                                                                                        >  2:2975861/1‑93 (MQ=255)

 ctctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgcc                                                                                        <  2:1373103/93‑1 (MQ=255)

  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgccc                                                                                       >  2:3400718/1‑92 (MQ=255)

    ttCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGctgcctgc                                                                                     <  1:3370942/93‑1 (MQ=255)

     tCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCa                                                                                    <  1:1976475/93‑1 (MQ=255)

     tCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCa                                                                                    >  2:319530/1‑93 (MQ=255)

      ccGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCaa                                                                                   <  1:4285525/93‑1 (MQ=255)

       cGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAt                                                                                  <  1:1551515/93‑1 (MQ=255)

        gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 <  1:3354456/93‑1 (MQ=255)

        gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 <  1:1208990/93‑1 (MQ=255)

        gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 >  1:2101431/1‑93 (MQ=255)

         aCTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTc                                                                                >  1:4289246/1‑93 (MQ=255)

         aCTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTc                                                                                >  2:11535/1‑93 (MQ=255)

          cTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCg                                                                               >  1:2845716/1‑93 (MQ=255)

           ttGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGc                                                                              >  2:1647941/1‑93 (MQ=255)

           ttGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGc                                                                              >  2:1733680/1‑93 (MQ=255)

            tGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCa                                                                             <  1:703966/93‑1 (MQ=255)

            tGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCa                                                                             <  1:141918/93‑1 (MQ=255)

               ttttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGAc                                                                          >  2:937066/1‑93 (MQ=255)

               ttttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGAc                                                                          <  1:1432304/93‑1 (MQ=255)

                  tCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCgg                                                                       >  2:2681329/1‑93 (MQ=255)

                  tCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCgg                                                                       >  2:2681328/1‑93 (MQ=255)

                      cTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTTCCTGCAATTCGCAGACCGGGAAt                                                                   <  2:4747709/93‑1 (MQ=255)

                       tCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATc                                                                  <  1:2758174/93‑1 (MQ=255)

                          ggTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaa                                                               >  2:3785328/1‑93 (MQ=255)

                          ggTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaa                                                               >  2:2934684/1‑93 (MQ=255)

                           gTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaaa                                                              <  1:2495397/93‑1 (MQ=255)

                           gTGAATCAATCGCTATATGTTGCGGAAAGAGCATCCGAAAGCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaaa                                                              <  1:4225930/93‑1 (MQ=255)

                            tGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAAt                                                             <  2:1498195/93‑1 (MQ=255)

                             gAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATg                                                            >  2:219256/1‑93 (MQ=255)

                             gAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATg                                                            >  2:4037395/1‑93 (MQ=255)

                                  aatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCaaa                                                       >  1:1224316/1‑93 (MQ=255)

                                   atcGGTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAAt                                                      <  2:1289060/93‑1 (MQ=255)

                                    tcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATg                                                     <  2:3905093/93‑1 (MQ=255)

                                     cGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGt                                                    >  1:2165001/1‑93 (MQ=255)

                                     cGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGt                                                    <  1:4484820/93‑1 (MQ=255)

                                      gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   >  2:1184573/1‑93 (MQ=255)

                                      gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   >  2:4319241/1‑93 (MQ=255)

                                       cTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTAt                                                  <  2:974752/93‑1 (MQ=255)

                                       cTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTAt                                                  >  1:999483/1‑93 (MQ=255)

                                        tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATc                                                 >  1:4709158/1‑93 (MQ=255)

                                        tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATc                                                 >  1:1459943/1‑93 (MQ=255)

                                        tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAACGTagc                                                 >  2:4336987/1‑91 (MQ=255)

                                         ttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATcc                                                <  1:3974071/93‑1 (MQ=255)

                                          tCTGTTGCGGAAAGAGCATCCGAAACCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCa                                               <  1:4327318/93‑1 (MQ=255)

                                          tCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCa                                               >  1:44880/1‑93 (MQ=255)

                                            tGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATc                                             <  2:3213299/93‑1 (MQ=255)

                                            tGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATc                                             >  2:2719118/1‑93 (MQ=255)

                                             gTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCa                                            >  2:1144593/1‑93 (MQ=255)

                                             gTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCa                                            >  1:333915/1‑93 (MQ=255)

                                                 cGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAg                                        >  2:971885/1‑93 (MQ=255)

                                                 cGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAg                                        <  2:3112262/93‑1 (MQ=255)

                                                  ggAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGt                                       >  2:1399641/1‑93 (MQ=255)

                                                   gAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCCCCAGtt                                      >  1:406749/1‑93 (MQ=255)

                                                   gAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGtt                                      >  1:406750/1‑93 (MQ=255)

                                                   gAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGtt                                      >  2:874472/1‑93 (MQ=255)

                                                    aaaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGttc                                     >  2:148145/1‑92 (MQ=255)

                                                     aaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGCATCCATCACCAGTTGc                                    >  1:2716557/1‑93 (MQ=255)

                                                     aaGAGCAGCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGc                                    <  2:1883527/93‑1 (MQ=255)

                                                       gagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGccc                                  >  2:2998659/1‑93 (MQ=255)

                                                       gagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGccc                                  >  1:3523757/1‑93 (MQ=255)

                                                        agCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCg                                 <  1:3548021/93‑1 (MQ=255)

                                                          cATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGTGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGc                               <  1:3077797/93‑1 (MQ=255)

                                                           aTCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCg                              >  1:1644831/1‑93 (MQ=255)

                                                               gAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAAt                          >  2:1052142/1‑93 (MQ=255)

                                                                aaaaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATa                         >  2:4570469/1‑93 (MQ=255)

                                                                aaaaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATa                         >  1:3838845/1‑93 (MQ=255)

                                                                aaaaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATa                         >  2:1189548/1‑93 (MQ=255)

                                                                 aaaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATAc                        >  1:2651169/1‑93 (MQ=255)

                                                                 aaaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATAc                        <  1:2363884/93‑1 (MQ=255)

                                                                  aaCCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACg                       >  1:4138424/1‑93 (MQ=255)

                                                                   aCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGACCGGCGTAATACGc                      <  2:4696337/93‑1 (MQ=255)

                                                                    ccGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGcc                     >  1:347744/1‑93 (MQ=255)

                                                                     cGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGcca                    >  1:3507930/1‑93 (MQ=255)

                                                                      gATAAACGCGCCGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGCCGTCATAAGccac                   >  2:643139/1‑93 (MQ=255)

                                                                      gATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccac                   >  1:2820901/1‑93 (MQ=255)

                                                                      gATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCAATCACCAGTTGCCCGGCGTAATACGccac                   >  1:473596/1‑93 (MQ=255)

                                                                         aaaCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcat                >  2:2823859/1‑93 (MQ=255)

                                                                           aCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatca              >  1:4700082/1‑93 (MQ=255)

                                                                           aCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatca              <  2:410965/93‑1 (MQ=255)

                                                                           aCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcataa              >  1:1659329/1‑91 (MQ=255)

                                                                           aCGCGAAGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGAATCCATCACCAGTTGCCCGGCGTAATACGCCACcatca              >  2:782323/1‑93 (MQ=255)

                                                                              cgacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatga           <  2:1930250/93‑1 (MQ=255)

                                                                               gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          >  2:2986966/1‑93 (MQ=255)

                                                                               gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          <  1:2787417/93‑1 (MQ=255)

                                                                               gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          <  2:2628238/93‑1 (MQ=255)

                                                                               gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          >  1:3656973/1‑93 (MQ=255)

                                                                               gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          <  2:2007271/93‑1 (MQ=255)

                                                                                acgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgatg         >  1:4339338/1‑93 (MQ=255)

                                                                                 cgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATgg        <  2:4596428/93‑1 (MQ=255)

                                                                                 cgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATgg        >  2:3282704/1‑93 (MQ=255)

                                                                                  gaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCCAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGc       >  2:2094295/1‑93 (MQ=255)

                                                                                   aTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCa      >  2:3767661/1‑93 (MQ=255)

                                                                                     tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    <  2:3852995/93‑1 (MQ=255)

                                                                                     tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    >  1:3000718/1‑93 (MQ=255)

                                                                                       aGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATa  <  1:752336/93‑1 (MQ=255)

GCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616602‑1616781

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