BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,185,099	T→C	S40S (AGT→AGC)	b2107 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,185,099	0	T	C	97.8%	582.9 / 0.5	184	S40S (AGT→AGC)	b2107	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/1); new base C (77/103); total (80/104)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 189 total aligned reads displayed.

GGTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCA  >  NC_000913/2185007‑2185188

ggTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagt                                                                                           <  2:1970591/93‑1 (MQ=255)

 gTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagcag                                                                                         <  1:80280/93‑1 (MQ=255)

     cgctttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGttt                                                                                      <  1:1662649/90‑1 (MQ=255)

      gctttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTc                                                                                     <  2:4142023/91‑1 (MQ=255)

       ctttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCa                                                                                    >  1:1708551/2‑93 (MQ=255)

        tttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCAt                                                                                   <  1:989380/93‑1 (MQ=255)

        tttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCAt                                                                                   >  1:2531970/1‑93 (MQ=255)

         ttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCAtt                                                                                  <  1:690874/93‑1 (MQ=255)

         ttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCAtt                                                                                  <  1:646804/93‑1 (MQ=255)

         ttGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCAtt                                                                                  >  1:2086219/1‑93 (MQ=255)

          tGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTg                                                                                 <  1:2862518/93‑1 (MQ=255)

            cTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCt                                                                               >  2:1987626/1‑93 (MQ=255)

             tGGTTCCCAGTGCCGCATGGACCGCACCAGCCACCGCGAGTTCGAACAATACCACGGGAATTTCTAAGTAAGAGTTAAGCAGTTTCATTACTg                                                                              <  2:4360556/93‑1 (MQ=255)

             tGGTTACCAGTTCCGCCTGGGCCGCACCACGCACCGCGGGTTCGGCCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGcgg                                                                              <  1:1603755/93‑3 (MQ=255)

             tGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTg                                                                              >  2:4533035/1‑93 (MQ=255)

             tGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTg                                                                              <  1:1603756/93‑1 (MQ=255)

             tGGTTACCAGTGCCGCCGCGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTg                                                                              <  1:1426595/93‑1 (MQ=255)

              ggTTACCAGTGCCGCCTGGGCCGCCCCAGCCACCGCGGGTTCGACCAATCCCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGa                                                                             <  2:2772847/93‑1 (MQ=255)

              ggTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGa                                                                             <  2:2900571/93‑1 (MQ=255)

               gTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGAc                                                                            >  2:666521/1‑93 (MQ=255)

                 tACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACtt                                                                          >  1:556266/1‑93 (MQ=255)

                 tACCAGTGCCGCCTGGACCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACtt                                                                          <  2:259676/93‑1 (MQ=255)

                  aCCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACttt                                                                         <  1:4230111/93‑1 (MQ=255)

                   ccAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTa                                                                        <  1:3922032/93‑1 (MQ=255)

                    cAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTaa                                                                       >  1:4511775/1‑93 (MQ=255)

                    cAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTaa                                                                       <  1:4414337/93‑1 (MQ=255)

                    cAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTaa                                                                       <  1:984367/93‑1 (MQ=255)

                    cAGTGCCGCCTGGGCCGCACCAGCCACCACGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTaa                                                                       <  1:216149/93‑1 (MQ=255)

                     aGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAg                                                                      <  1:500474/93‑1 (MQ=255)

                     aGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAg                                                                      >  1:1233225/1‑93 (MQ=255)

                        gccgccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCAt                                                                   >  2:4495525/1‑93 (MQ=255)

                         gcgccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCAtt                                                                  <  1:1752573/92‑1 (MQ=255)

                         ccgccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCAtt                                                                  <  1:2331181/93‑1 (MQ=255)

                         ccgccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCAtt                                                                  <  2:191104/93‑1 (MQ=255)

                          cgccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCAttt                                                                 >  1:3631791/1‑93 (MQ=255)

                           gccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTc                                                                <  2:3604763/93‑1 (MQ=255)

                           gccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTc                                                                <  1:4565491/93‑1 (MQ=255)

                           gccTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTc                                                                <  1:4603351/93‑1 (MQ=255)

                             cTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTcca                                                              >  2:454734/1‑91 (MQ=255)

                             cTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCaa                                                              >  2:2501060/1‑93 (MQ=255)

                               gggCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAAc                                                            >  2:3677262/1‑93 (MQ=255)

                               gggCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAAc                                                            >  2:510098/1‑93 (MQ=255)

                                  ccGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAg                                                         <  1:2566283/93‑1 (MQ=255)

                                  ccGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAg                                                         <  2:4292940/93‑1 (MQ=255)

                                   cGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAgg                                                        <  2:2158433/93‑1 (MQ=255)

                                    gCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAggg                                                       <  2:2062941/93‑1 (MQ=255)

                                      cccAGCCACCGAGGGTTCGGCCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCTGACCTgggaa                                                     <  2:4220724/92‑3 (MQ=255)

                                        cAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGaca                                                   >  1:2953652/1‑93 (MQ=255)

                                          gCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACAcc                                                 <  1:3240071/93‑1 (MQ=255)

                                           gcctcgcgGGTTCGACCAATTCCTTGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCg                                                <  1:854517/89‑1 (MQ=255)

                                             aCCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTa                                              >  1:3448466/1‑93 (MQ=255)

                                                gcgGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCa                                           <  2:2532938/93‑1 (MQ=255)

                                                  gggTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGa                                         <  2:842173/93‑1 (MQ=255)

                                                   ggTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGaa                                        <  2:2933573/93‑1 (MQ=255)

                                                       cGACCAATAACTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGATGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGACATgg                                    >  2:3343738/1‑92 (MQ=255)

                                                        gACCAATACCTCGGGAATTTCTAAGTATGAGTTATGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgta                                   <  1:2905923/93‑1 (MQ=255)

                                                        gACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtc                                   >  1:1159321/1‑92 (MQ=255)

                                                        gACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgta                                   <  2:155461/93‑1 (MQ=255)

                                                        gACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgta                                   >  1:1839347/1‑93 (MQ=255)

                                                        gACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgta                                   >  1:1499143/1‑93 (MQ=255)

                                                        gACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgta                                   <  2:4082045/93‑1 (MQ=255)

                                                         aCCAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtac                                  <  1:507225/93‑1 (MQ=255)

                                                         aCCAATACCTCGAGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtac                                  <  2:988234/93‑1 (MQ=255)

                                                           cAATACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccg                                <  2:4604733/93‑1 (MQ=255)

                                                             aTACCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgca                              >  2:350534/1‑91 (MQ=255)

                                                               aCCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgtacc                            >  1:2885296/1‑93 (MQ=255)

                                                               aCCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgtacc                            >  2:647832/1‑93 (MQ=255)

                                                               aCCTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATgtaccgtacc                            >  2:3201370/1‑93 (MQ=255)

                                                                ccTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGtaccgtaccg                           >  1:1819636/1‑93 (MQ=255)

                                                                 cTCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGCCTCCGTACCAGATATGTCCCGTTCCGa                          <  2:4669171/93‑1 (MQ=255)

                                                                  tCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGAt                         <  2:1292456/93‑1 (MQ=255)

                                                                  tCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGAt                         >  2:1167661/1‑93 (MQ=255)

                                                                  tCGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGAt                         >  2:2659811/1‑93 (MQ=255)

                                                                   cGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTGCCAGAAATGTACCGTACCGATg                        <  1:3302106/93‑1 (MQ=255)

                                                                   cGGGAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATg                        <  2:4021271/93‑1 (MQ=255)

                                                                    gggAATTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGa                       >  1:1176373/1‑93 (MQ=255)

                                                                       aaTTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTa                    >  2:4008277/1‑93 (MQ=255)

                                                                       aaTTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTa                    <  1:1075111/93‑1 (MQ=255)

                                                                        aTTTCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTac                   <  2:135682/93‑1 (MQ=255)

                                                                          ttCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaa                 >  1:2445731/1‑93 (MQ=255)

                                                                           tCTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaac                >  1:2224349/1‑93 (MQ=255)

                                                                            cTAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTacaaca               >  1:156871/1‑93 (MQ=255)

                                                                             tAAGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACAt              >  1:2025906/1‑93 (MQ=255)

                                                                              aaGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACAtt             >  2:3993112/1‑93 (MQ=255)

                                                                               aGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTa            >  1:2890446/1‑93 (MQ=255)

                                                                               aGTATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTa            >  2:766629/1‑93 (MQ=255)

                                                                                 tATGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGCAATGTACCGTACCGATGAGTACAACATTAAg          >  2:208114/1‑93 (MQ=255)

                                                                                  aTGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGc         <  1:1477058/93‑1 (MQ=255)

                                                                                   tGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCa        <  1:1047658/93‑1 (MQ=255)

                                                                                   tGAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCa        <  2:3258632/93‑1 (MQ=255)

                                                                                    gAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAg       >  2:387104/1‑93 (MQ=255)

                                                                                    gAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAg       <  1:3511568/93‑1 (MQ=255)

                                                                                    gAGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAg       >  2:2286718/1‑93 (MQ=255)

                                                                                     aGTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGt      >  1:3338783/1‑93 (MQ=255)

                                                                                      gTTAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTg     <  2:2272745/93‑1 (MQ=255)

                                                                                        tAAGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGc   <  1:375935/93‑1 (MQ=255)

                                                                                         aaGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCa  <  1:3959251/93‑1 (MQ=255)

                                                                                         aaGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCa  >  1:4133974/1‑93 (MQ=255)

                                                                                         aaGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCa  <  1:4543788/93‑1 (MQ=255)

                                                                                         aaGCAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCa  <  2:1679882/93‑1 (MQ=255)

GGTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCA  >  NC_000913/2185007‑2185188

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