BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,745,316	C→T	I201I (ATC→ATT)	lyxK →	L‑xylulose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,745,316	0	C	T	99.6%	810.0 / ‑5.6	230	I201I (ATC→ATT)	lyxK	L‑xylulose kinase
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (135/93); total (135/95)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

Only 100 of 241 total aligned reads displayed.

AGGCTGTGAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATCGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCGGTCT  >  NC_000913/3745226‑3745405

aGGCTGTGAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGc                                                                                         <  2:3182016/93‑1 (MQ=255)

 ggCTGTGAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCt                                                                                        >  1:3111/1‑93 (MQ=255)

 ggCTGTGAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCt                                                                                        <  1:2780950/93‑1 (MQ=255)

 ggCTGTGAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCt                                                                                        >  2:2756030/1‑93 (MQ=255)

   cTGTGAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGa                                                                                      >  2:695791/1‑93 (MQ=255)

     gtgAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGaaa                                                                                    <  2:1229528/93‑1 (MQ=255)

       gAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAaatc                                                                                  >  1:4457130/1‑93 (MQ=255)

       gAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAaatc                                                                                  >  1:1814540/1‑93 (MQ=255)

       gAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAaatc                                                                                  >  2:648780/1‑93 (MQ=255)

       gAAGAGAGCAATATTTCAGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATAGCTGAaatc                                                                                  <  1:3545978/93‑1 (MQ=255)

        aaGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAaatca                                                                                 <  2:4637002/93‑1 (MQ=255)

         agagagCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAaatcca                                                                                >  1:1768447/1‑91 (MQ=255)

          gagagCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGGaatcaat                                                                               >  2:3397254/1‑93 (MQ=255)

          gagagCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAaatcaat                                                                               <  2:2965039/93‑1 (MQ=255)

             agCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCAc                                                                            <  1:3216628/93‑1 (MQ=255)

              gCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACg                                                                           <  1:893549/93‑1 (MQ=255)

              gCAATATTTCCGAGTCCAACCTCTACAACATGAGCCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACg                                                                           >  1:4078840/1‑93 (MQ=255)

               cAATATTTCCGAGTCCACCCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGc                                                                          <  2:1628886/93‑1 (MQ=255)

               cAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGc                                                                          <  1:1712561/93‑1 (MQ=255)

               cAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGc                                                                          <  2:2336459/93‑1 (MQ=255)

                 atatTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGccc                                                                        <  2:13293/93‑1 (MQ=255)

                 atatTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGccc                                                                        >  2:4067509/1‑93 (MQ=255)

                 atatTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGccc                                                                        <  2:982181/93‑1 (MQ=255)

                  tatTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCt                                                                       >  1:597836/1‑93 (MQ=255)

                  tatTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCt                                                                       <  2:538082/93‑1 (MQ=255)

                    tttCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTgc                                                                     <  2:3658862/93‑1 (MQ=255)

                    tttCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTgc                                                                     >  2:3755744/1‑93 (MQ=255)

                       ccGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTgccgc                                                                  >  2:1138487/1‑93 (MQ=255)

                        cGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTgccgcc                                                                 >  2:3613124/1‑93 (MQ=255)

                           gTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTcccccccct                                                              >  1:603731/1‑90 (MQ=255)

                           gTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCtgt                                                              <  2:2501842/93‑1 (MQ=255)

                           gTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCtgt                                                              <  1:30014/93‑1 (MQ=255)

                           gTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCcgt                                                              >  2:3996669/1‑93 (MQ=255)

                              cAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCtgttgt                                                           >  1:1741373/1‑93 (MQ=255)

                              cAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCtgttgt                                                           <  2:2228359/93‑1 (MQ=255)

                                aCCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCg                                                         >  2:1762520/1‑93 (MQ=255)

                                aCCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCg                                                         <  1:1989879/93‑1 (MQ=255)

                                 ccTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCgg                                                        >  2:2694034/1‑93 (MQ=255)

                                   tctACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGAt                                                      >  1:1153414/1‑93 (MQ=255)

                                   tctACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGAt                                                      >  1:3151316/1‑93 (MQ=255)

                                     tACAACATGAGTCTTGGGGGATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATcg                                                    >  1:1854781/1‑92 (MQ=255)

                                     tACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCt                                                    >  2:3779450/1‑93 (MQ=255)

                                     tACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCt                                                    >  1:3451077/1‑93 (MQ=255)

                                      acaacaTGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTg                                                   >  1:4045130/1‑93 (MQ=255)

                                       caacaTGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGc                                                  <  2:2639888/93‑1 (MQ=255)

                                         acaTGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCg                                                >  2:1897796/1‑93 (MQ=255)

                                         acaTGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCg                                                >  2:140604/1‑93 (MQ=255)

                                          caTGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGa                                               >  1:4330624/1‑93 (MQ=255)

                                           aTGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGaa                                              >  2:329104/1‑93 (MQ=255)

                                           aTGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGaa                                              >  1:3353915/1‑93 (MQ=255)

                                           aTGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGaa                                              <  1:3816606/93‑1 (MQ=255)

                                            tGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGaaa                                             >  1:4167614/1‑93 (MQ=255)

                                             gAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAAt                                            <  2:3347617/93‑1 (MQ=255)

                                              aGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATc                                           <  2:3560481/93‑1 (MQ=255)

                                               gTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCt                                          >  1:4677711/1‑93 (MQ=255)

                                                  ttGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCg                                       >  2:2562740/1‑93 (MQ=255)

                                                  ttGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCg                                       <  1:2396563/93‑1 (MQ=255)

                                                  ttGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAAACACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCg                                       >  1:2136754/1‑93 (MQ=255)

                                                   tGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCgg                                      <  1:692302/93‑1 (MQ=255)

                                                    ggtgAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAACCAATCCCGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCggg                                     <  2:2076489/90‑1 (MQ=255)

                                                    ggggAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCggg                                     >  1:2672260/1‑93 (MQ=255)

                                                    ggggAAGATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCggg                                     <  2:38164/93‑1 (MQ=255)

                                                     gggAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCgggg                                    <  1:3760868/93‑1 (MQ=255)

                                                     gggAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCgggg                                    >  1:3955449/1‑93 (MQ=255)

                                                      ggAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGga                                   >  1:1020606/1‑93 (MQ=255)

                                                      ggAATATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGga                                   <  1:2280311/93‑1 (MQ=255)

                                                       gAATATGACCCGTGCCTCACCGACTGGCTGGGGATCGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGgag                                  <  1:4490338/93‑1 (MQ=255)

                                                        aaTATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGgaga                                 <  2:1653943/93‑1 (MQ=255)

                                                        aaTATGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGgaga                                 <  1:966411/93‑1 (MQ=255)

                                                        aaTATGACCCGTGCCTCACCGACTGGCTCGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGgaga                                 <  2:2575497/93‑1 (MQ=255)

                                                         atatGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGAt                                >  2:2428390/1‑93 (MQ=255)

                                                          tatGACCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATc                               >  2:4565571/1‑93 (MQ=255)

                                                              aCCCGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCg                           >  1:1340532/1‑93 (MQ=255)

                                                               cccGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGc                          <  1:4140596/93‑1 (MQ=255)

                                                                ccGTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCt                         >  2:1174803/1‑93 (MQ=255)

                                                                 cGTGCCTCTCCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTCTTGTCGGATCTGCCTAAATCTGCGGGGAGATCACCGCTc                        <  2:1559387/93‑1 (MQ=255)

                                                                 cGTGCCTCACCGACTGGCTGGTGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTc                        <  1:702736/93‑1 (MQ=255)

                                                                  gTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCa                       >  2:2615521/1‑93 (MQ=255)

                                                                  gTGCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCa                       >  2:2878824/1‑93 (MQ=255)

                                                                    gCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGa                     >  1:4636679/1‑93 (MQ=255)

                                                                    gCCTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGa                     <  1:1569869/93‑1 (MQ=255)

                                                                      cTCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACa                   <  1:2481176/93‑1 (MQ=255)

                                                                       tCACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGATAg                  >  2:434376/1‑93 (MQ=255)

                                                                        cgccGACTGTCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAgc                 <  1:2525081/91‑1 (MQ=255)

                                                                        cACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAgc                 <  2:2383088/93‑1 (MQ=255)

                                                                        cACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAgc                 >  2:3932442/1‑93 (MQ=255)

                                                                        cACCGACTGGCTGGGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAgc                 <  1:2801300/93‑1 (MQ=255)

                                                                               tggctggGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCCGACAGCCGCCCTg          >  2:199748/1‑93 (MQ=255)

                                                                               tggctggGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTg          <  2:3886286/93‑1 (MQ=255)

                                                                                ggctggGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGa         >  2:4151500/1‑93 (MQ=255)

                                                                                ggctggGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGa         >  1:1258816/1‑93 (MQ=255)

                                                                                 gatggGGATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGAc        <  1:959700/91‑1 (MQ=255)

                                                                                    ggggATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCgg     <  1:1629017/93‑1 (MQ=255)

                                                                                    ggggATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCgg     <  1:3502565/93‑1 (MQ=255)

                                                                                    ggggATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCgg     >  2:2800259/1‑93 (MQ=255)

                                                                                    ggggATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCgg     >  2:772483/1‑93 (MQ=255)

                                                                                      ggATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCGGTc   <  2:2269476/93‑1 (MQ=255)

                                                                                      ggATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCGGTc   >  1:1662927/1‑93 (MQ=255)

                                                                                       gATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCGGTCt  >  1:1441235/1‑93 (MQ=255)

                                                                                       gATTGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCGGTCt  >  2:2912727/1‑93 (MQ=255)

AGGCTGTGAAGAGAGCAATATTTCCGAGTCCAACCTCTACAACATGAGTCTTGGGGAATATGACCCGTGCCTCACCGACTGGCTGGGGATCGCTGAAATCAATCACGCCCTGCCGCCTGTTGTCGGATCTGCCGAAATCTGCGGGGAGATCACCGCTCAGACAGCCGCCCTGACCGGTCT  >  NC_000913/3745226‑3745405

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