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	99.0%	358.7 / ‑5.2	101	G81S (GGC→AGC)	ydeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (50/50); total (50/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 105 total aligned reads displayed.

GCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616602‑1616781

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

 ctctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgcc                                                                                        >  1:2133651/1‑93 (MQ=255)

  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       <  1:1161280/93‑1 (MQ=255)

  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       >  2:513249/1‑93 (MQ=255)

  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       >  1:4444226/1‑93 (MQ=255)

   ctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTg                                                                                      >  1:1778170/1‑93 (MQ=255)

   ctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTg                                                                                      >  1:4195338/1‑93 (MQ=255)

   ctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTg                                                                                      >  2:1530234/1‑93 (MQ=255)

    ttCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGctgcctgc                                                                                     >  1:363661/1‑93 (MQ=255)

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

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

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

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

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

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

           ttGCTTTTCGCCTCCGGTGAATCAATCGGTTTCTGTTGCGGAAAGAGGATGCGAAAACGGATAAACGCAACGATTAAGCTGCCTGCAATTCGc                                                                              >  2:6344/1‑93 (MQ=255)

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

              cTTTTCGCCTCCAGTGAATCAATCGCTTTCTGTTGCGGAAAGTGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGa                                                                           <  2:2263147/93‑1 (MQ=255)

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

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

                  tCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCgg                                                                       <  2:356502/93‑1 (MQ=255)

                    gCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGa                                                                     <  2:4396724/93‑1 (MQ=255)

                     ccTCCGGTGAATCAATCTCTTTCTGTTGCGGAAAGAGCATCCGAAAGCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGaa                                                                    <  2:360489/93‑1 (MQ=255)

                     ccTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGaa                                                                    >  2:2806568/1‑93 (MQ=255)

                       tCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATc                                                                  >  1:2872479/1‑93 (MQ=255)

                         cGGTGAATCAATCGCTTTCTGTTGCGGAAATAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGa                                                                <  2:1568811/93‑1 (MQ=255)

                           gTGAATCATTCGCTTTCTGTGGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaaa                                                              <  2:3959628/93‑1 (MQ=255)

                           gTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGcaa                                                              >  1:3376471/1‑93 (MQ=255)

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

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

                           gTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaaa                                                              >  1:2014377/1‑93 (MQ=255)

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

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

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

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

                               atcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGcc                                                          <  1:327726/93‑1 (MQ=255)

                               atcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGcc                                                          >  2:3983590/1‑93 (MQ=255)

                                tcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCa                                                         <  2:416788/93‑1 (MQ=255)

                                 caatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCaa                                                        >  2:491214/1‑93 (MQ=255)

                                  aatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCaaa                                                       <  2:4113456/93‑1 (MQ=255)

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

                                    tcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATg                                                     >  2:53431/1‑93 (MQ=255)

                                      gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   <  2:4057887/93‑1 (MQ=255)

                                      gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   <  1:1111092/93‑1 (MQ=255)

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

                                        tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCCCAGACCGGGAATCGAAATGCCAAATGTATc                                                 >  2:3337494/1‑93 (MQ=255)

                                         ttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATcc                                                >  2:481329/1‑93 (MQ=255)

                                           cTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCAt                                              <  1:9992/93‑1 (MQ=255)

                                           cTGTTGCGGAAAGAGCACCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCAt                                              <  1:1798431/93‑1 (MQ=255)

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

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

                                               tGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCAcc                                          <  1:3819833/93‑1 (MQ=255)

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

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

                                                  ggAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGt                                       <  2:34446/93‑1 (MQ=255)

                                                  ggAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGt                                       <  1:292664/93‑1 (MQ=255)

                                                   gAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGtt                                      <  1:2811426/93‑1 (MQ=255)

                                                    aaaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTg                                     >  2:2055117/1‑93 (MQ=255)

                                                    aaaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTg                                     <  1:2587161/93‑1 (MQ=255)

                                                    aaaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTg                                     <  1:1601286/93‑1 (MQ=255)

                                                     aaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGc                                    >  1:4232924/1‑93 (MQ=255)

                                                      agagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGcc                                   >  2:1536873/1‑93 (MQ=255)

                                                      agagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGcc                                   >  2:4410394/1‑93 (MQ=255)

                                                      agagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGcc                                   <  1:544993/93‑1 (MQ=255)

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

                                                         gCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCgg                                <  1:1339619/93‑1 (MQ=255)

                                                         gCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCgg                                >  2:1004632/1‑93 (MQ=255)

                                                            tCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCga                             >  2:3033733/1‑92 (MQ=255)

                                                            tCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCga                             >  1:399633/1‑92 (MQ=255)

                                                              cGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTaa                           <  2:3019098/93‑1 (MQ=255)

                                                               gAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAAt                          <  1:2043260/93‑1 (MQ=255)

                                                                aaaaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATa                         <  1:1509330/93‑1 (MQ=255)

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

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

                                                                  aaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACg                       >  1:3878471/1‑93 (MQ=255)

                                                                  aaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACg                       >  1:2015523/1‑93 (MQ=255)

                                                                   aCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGc                      >  1:903363/1‑93 (MQ=255)

                                                                     cGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGCAATACCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccc                    >  2:609715/1‑92 (MQ=255)

                                                                     cGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGcca                    <  2:277550/93‑1 (MQ=255)

                                                                      ggtAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccac                   >  1:3667818/3‑93 (MQ=255)

                                                                      gATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccac                   <  1:2233021/93‑1 (MQ=255)

                                                                       aTAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccacc                  >  2:3850367/1‑93 (MQ=255)

                                                                          aaCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatc               >  2:507209/1‑93 (MQ=255)

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

                                                                            cgcgACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatcat             >  1:1114207/1‑93 (MQ=255)

                                                                             gcgACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatg            >  1:2169557/1‑93 (MQ=255)

                                                                              cgacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatga           >  2:4673571/1‑93 (MQ=255)

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

                                                                                 cgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGCAATACGCCACCATCATGATgg        >  2:55840/1‑93 (MQ=255)

                                                                                  gaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGc       >  1:1862456/1‑93 (MQ=255)

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

                                                                                   aTTAAGCTGCCTGCAATTCGCAGACAGGGAATCGAAATGCCAAATGCATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCa      >  2:2785126/1‑93 (MQ=255)

                                                                                    ttAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaa     >  1:2251152/1‑93 (MQ=255)

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

                                                                                     tAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    <  1:1937958/93‑1 (MQ=255)

                                                                                      aaGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAAt   <  2:353349/93‑1 (MQ=255)

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

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

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

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

GCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616602‑1616781

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: `‑`