BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	827,018	C→T	S396N (AGT→AAT)	ybhF ←	YbhF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	827,018	0	C	T	99.1%	400.5 / ‑5.0	114	S396N (AGT→AAT)	ybhF	YbhF
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (52/61); total (52/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 121 total aligned reads displayed.

AAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTACCC  >  NC_000913/826930‑827107

aaaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttctc                                                                                       >  2:771409/1‑93 (MQ=255)

aaaGCGTAAATTCTGTTCGACCGTCAGGTTACCGAAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttctc                                                                                       >  2:3659565/1‑93 (MQ=255)

 aaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttctct                                                                                      <  2:816811/93‑1 (MQ=255)

 aaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttctct                                                                                      >  1:3984190/1‑93 (MQ=255)

  aGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCtt                                                                                     >  1:1174798/1‑93 (MQ=255)

  aGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCtt                                                                                     >  2:1916703/1‑93 (MQ=255)

   gCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGAGCGCTTTACCGGAATTCTCttt                                                                                    >  2:2747491/1‑93 (MQ=255)

    cGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTc                                                                                   <  1:379761/93‑1 (MQ=255)

    cGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTc                                                                                   <  1:148663/93‑1 (MQ=255)

     gTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCa                                                                                  <  1:4183475/93‑1 (MQ=255)

       aaaTTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGa                                                                                <  2:2797841/93‑1 (MQ=255)

        aaTTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGat                                                                               <  1:2114509/93‑1 (MQ=255)

          ttCTGTTCGACCGTCAGGGTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatcc                                                                             >  2:26929/1‑93 (MQ=255)

           tCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatcca                                                                            >  2:4214185/1‑93 (MQ=255)

           tCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatcca                                                                            <  1:4396305/93‑1 (MQ=255)

           tCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatcca                                                                            <  2:1601434/93‑1 (MQ=255)

            cTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccat                                                                           <  2:4673647/93‑1 (MQ=255)

             tGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccatc                                                                          >  1:3634251/1‑93 (MQ=255)

             tGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccatc                                                                          >  1:4248850/1‑93 (MQ=255)

             tGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccatc                                                                          <  2:1298179/93‑1 (MQ=255)

             tGTTCGACAGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccatc                                                                          >  1:2871723/1‑93 (MQ=255)

              gTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccatcc                                                                         >  1:4200036/1‑93 (MQ=255)

               ttCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATccc                                                                        <  1:356181/93‑1 (MQ=255)

                  gACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAg                                                                     <  1:953856/93‑1 (MQ=255)

                   aCCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGc                                                                    >  1:3113752/1‑93 (MQ=255)

                      gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCAcc                                                                 >  1:4375274/1‑93 (MQ=255)

                         aGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAgc                                                              >  1:1062902/1‑93 (MQ=255)

                         aGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAgc                                                              <  2:1862593/93‑1 (MQ=255)

                         aGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAgc                                                              <  2:2068357/93‑1 (MQ=255)

                          ggTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAgcg                                                             >  2:1789957/1‑93 (MQ=255)

                           gTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAgcgc                                                            >  2:1871207/1‑93 (MQ=255)

                             tACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCt                                                          <  1:1552925/93‑1 (MQ=255)

                             tACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCt                                                          <  2:1125103/93‑1 (MQ=255)

                              tccGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTg                                                         <  1:3733511/92‑1 (MQ=255)

                              aCCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTg                                                         >  2:678622/1‑93 (MQ=255)

                              aCCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTg                                                         >  1:1303319/1‑93 (MQ=255)

                                 gTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGcc                                                      <  2:2598281/93‑1 (MQ=255)

                                  tAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCg                                                     <  1:1318639/93‑1 (MQ=255)

                                  tAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCg                                                     <  2:3549920/93‑1 (MQ=255)

                                  tAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCg                                                     >  2:2884731/1‑93 (MQ=255)

                                    gagCGAAATTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGa                                                   <  2:1764061/93‑1 (MQ=255)

                                    gagCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGa                                                   >  1:1077719/1‑93 (MQ=255)

                                     agCGAAAATTTTTGCGCCATACAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGaa                                                  <  2:3055780/93‑1 (MQ=255)

                                       cGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGt                                                <  2:2816487/93‑1 (MQ=255)

                                        gAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTc                                               <  2:4271431/93‑1 (MQ=255)

                                         aaaaTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCg                                              >  2:2158197/1‑93 (MQ=255)

                                            aTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCa                                           >  2:3716436/1‑93 (MQ=255)

                                            aTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCa                                           <  1:3358561/93‑1 (MQ=255)

                                                  gcgcCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCa                                     >  2:4213515/1‑93 (MQ=255)

                                                  gcgcCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCa                                     >  2:111911/1‑93 (MQ=255)

                                                   cgcCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCaa                                    <  1:4046326/93‑1 (MQ=255)

                                                       atatAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCg                                >  2:29665/1‑93 (MQ=255)

                                                       atatAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCg                                <  1:4405279/93‑1 (MQ=255)

                                                        tataGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGc                               <  1:3673182/93‑1 (MQ=255)

                                                        tataGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGc                               <  2:2903854/93‑1 (MQ=255)

                                                         ataGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGcc                              >  1:4143680/1‑92 (MQ=255)

                                                          taGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGcgc                             >  2:3417008/1‑91 (MQ=255)

                                                          taGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGcac                             >  1:1670589/1‑93 (MQ=255)

                                                          taGCCGAGATGCTGGCGCGCTTTACCAGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGcac                             >  2:837799/1‑93 (MQ=255)

                                                           aGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGcaca                            <  2:4455896/93‑1 (MQ=255)

                                                           aGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGcaca                            <  1:2536361/93‑1 (MQ=255)

                                                             ccGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCCCCAGCGCATGGCCGGAAGTCGGCACCAGCAAACCGCAcatc                          <  1:3365765/93‑1 (MQ=255)

                                                              cGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatca                         >  1:4006693/1‑93 (MQ=255)

                                                              cGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatca                         <  2:815753/93‑1 (MQ=255)

                                                              cGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatca                         >  2:4162942/1‑93 (MQ=255)

                                                               gagaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatcat                        <  1:3000951/93‑1 (MQ=255)

                                                                agaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACatcatc                       >  1:4503095/1‑93 (MQ=255)

                                                                 gaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCt                      >  2:2242539/1‑93 (MQ=255)

                                                                 gaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATACCCAGCACCAGCGCCTGGGCGGAAGTCGGCACCAGCGAACCGCACATCATCt                      >  2:3137021/1‑93 (MQ=255)

                                                                  aTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACCTCATCtt                     <  2:4454779/93‑1 (MQ=255)

                                                                   tGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTa                    >  1:1133924/1‑93 (MQ=255)

                                                                     cTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTaaa                  <  2:2254867/93‑1 (MQ=255)

                                                                      tGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAg                 <  2:320320/93‑1 (MQ=255)

                                                                       ggCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAgg                <  2:1451487/93‑1 (MQ=255)

                                                                         cgcgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggtg              <  2:4071979/93‑1 (MQ=255)

                                                                          gcgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggtgg             >  1:1228768/1‑93 (MQ=255)

                                                                          gcgcTTTACCGGAATTCTCTTTCAGATCCAACCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggtgg             <  1:4622168/93‑1 (MQ=255)

                                                                           cgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggtggt            <  2:457864/93‑1 (MQ=255)

                                                                           cgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggtggt            >  2:121786/1‑93 (MQ=255)

                                                                            gcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTc           >  2:4633063/1‑93 (MQ=255)

                                                                            gcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTc           >  2:1655620/1‑93 (MQ=255)

                                                                            gcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTc           >  2:2419459/1‑93 (MQ=255)

                                                                              tttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGa         <  2:3570532/93‑1 (MQ=255)

                                                                              tttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGa         <  2:1456100/93‑1 (MQ=255)

                                                                              tttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGa         >  2:4647587/1‑93 (MQ=255)

                                                                              tttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGa         <  2:1757774/93‑1 (MQ=255)

                                                                               ttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAt        <  2:645575/93‑1 (MQ=255)

                                                                               ttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAt        <  2:4463975/93‑1 (MQ=255)

                                                                                tACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGTCGGTCGAtt       >  2:3297190/1‑93 (MQ=255)

                                                                                tACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAtt       <  2:4575438/93‑1 (MQ=255)

                                                                                tACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAtt       >  2:191129/1‑93 (MQ=255)

                                                                                tACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAtt       >  2:213271/1‑93 (MQ=255)

                                                                                   cGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAc    <  2:3796756/93‑1 (MQ=255)

                                                                                   cGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAc    <  1:802137/93‑1 (MQ=255)

                                                                                    ggAATTCTCTTTCAGATCCATCCCCAGCGCCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAcc   >  1:2968958/1‑93 (MQ=255)

                                                                                    ggAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAcc   >  2:2782117/1‑93 (MQ=255)

                                                                                     gAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAccc  >  1:2413805/1‑93 (MQ=255)

                                                                                     gAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAccc  <  2:1562284/93‑1 (MQ=255)

                                                                                     gAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAccc  >  2:2932506/1‑93 (MQ=255)

                                                                                     gAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAccc  >  1:2864227/1‑93 (MQ=255)

AAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTACCC  >  NC_000913/826930‑827107

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