BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,744,958	G→A	G82E (GGA→GAA)	lyxK →	L‑xylulose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,744,958	0	G	A	99.2%	901.6 / ‑5.4	251	G82E (GGA→GAA)	lyxK	L‑xylulose kinase
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (110/139); total (111/140)
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.33e-01

Only 100 of 267 total aligned reads displayed.

GGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGGAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAATCGTTCGTC  >  NC_000913/3744868‑3745047

ggCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGaaa                                                                                         <  2:762671/93‑1 (MQ=255)

 gCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCTGCATCTCCGCACAGGaaaa                                                                                        <  2:1768442/93‑1 (MQ=255)

 gCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGaaaa                                                                                        <  2:2903338/93‑1 (MQ=255)

 gCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGaaaa                                                                                        <  1:2237011/93‑1 (MQ=255)

  cAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGGAAAg                                                                                       >  2:4219348/1‑93 (MQ=255)

  cAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAg                                                                                       <  2:4242770/93‑1 (MQ=255)

   aaTGCTGCATGGCTGTCATTCGCTCCTTGCTTACGCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAgg                                                                                      <  2:818464/93‑1 (MQ=255)

    aTGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAggg                                                                                     >  1:1200505/1‑93 (MQ=255)

    aTGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAggg                                                                                     <  2:799791/93‑1 (MQ=255)

    aTGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAggg                                                                                     >  2:2759297/1‑93 (MQ=255)

    aTGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAggg                                                                                     <  1:3037186/93‑1 (MQ=255)

    aTGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAggg                                                                                     <  1:1723955/93‑1 (MQ=255)

     tgctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGc                                                                                    >  2:717537/1‑93 (MQ=255)

      gctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCt                                                                                   <  2:3208739/93‑1 (MQ=255)

       ctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCtt                                                                                  <  2:3425980/93‑1 (MQ=255)

        tgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTg                                                                                 >  1:18623/1‑93 (MQ=255)

        tgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTg                                                                                 <  2:26582/93‑1 (MQ=255)

         gcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGt                                                                                >  2:4126874/1‑93 (MQ=255)

          cATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGtt                                                                               <  2:4179429/93‑1 (MQ=255)

          cATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGtt                                                                               <  1:1570071/93‑1 (MQ=255)

           aTGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttt                                                                              >  2:4271162/1‑93 (MQ=255)

           aTGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttt                                                                              <  2:64022/93‑1 (MQ=255)

           aTGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttt                                                                              >  2:551663/1‑93 (MQ=255)

            tGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGtttt                                                                             <  2:4558709/93‑1 (MQ=255)

             ggCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttttt                                                                            <  1:4269643/93‑1 (MQ=255)

             ggCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttttt                                                                            <  2:2313567/93‑1 (MQ=255)

             ggCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttttt                                                                            <  1:3616199/93‑1 (MQ=255)

              gCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGGACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTg                                                                           >  2:1147077/1‑93 (MQ=255)

                 gTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTg                                                                        >  2:3059672/1‑93 (MQ=255)

                 gTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTg                                                                        >  1:1925016/1‑93 (MQ=255)

                   cATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGa                                                                      <  1:2344608/93‑1 (MQ=255)

                       cgcgCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATaaa                                                                  <  1:1348989/93‑1 (MQ=255)

                        gcgcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATaaaa                                                                 <  2:3275091/93‑1 (MQ=255)

                        gcgcCCTGCTTACTCAGTCTGGTGTTAGCGGGGGACAAATTGTCGGTATCGGCATCTCCGCACAGGGAAAGGGCTTGTGTTTGCTGGCTaaaa                                                                 >  2:3613811/1‑93 (MQ=255)

                          gcgcTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGATGGATAAAAAc                                                               <  1:1493684/90‑1 (MQ=255)

                          gcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAAc                                                               <  1:3292643/93‑1 (MQ=255)

                          gcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAAc                                                               >  1:4420346/1‑93 (MQ=255)

                          gcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAAc                                                               >  2:895320/1‑93 (MQ=255)

                          gcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAAc                                                               <  1:658734/93‑1 (MQ=255)

                          ccccTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAAc                                                               <  1:2457992/92‑1 (MQ=255)

                          ccccTGCTTACTCATTCTGGTGTTAGCGGCGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAAc                                                               <  2:2989040/92‑1 (MQ=255)

                            ccTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGa                                                             >  1:404881/1‑93 (MQ=255)

                            ccTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGa                                                             >  1:4138538/1‑93 (MQ=255)

                             cTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGAc                                                            >  1:3148992/1‑93 (MQ=255)

                              tGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGGAAAGGGCTTGTTTTTGCTGGATAATAACGACa                                                           >  1:927908/1‑93 (MQ=255)

                                cTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACaaa                                                         <  1:3333912/93‑1 (MQ=255)

                                 ttACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAAc                                                        <  1:1009020/93‑1 (MQ=255)

                                    cTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGc                                                     <  1:1547588/93‑1 (MQ=255)

                                    cTCATTCTGGTGTTAGCGGAGAACAAATTGTCGGTATCGGCATCTCCGCAAAGGAAACGGTCTTGTTTTTGCTGGATAAAACCGCCAAACCGc                                                     <  2:1574650/93‑1 (MQ=255)

                                      cATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTc                                                   <  2:427417/93‑1 (MQ=255)

                                      cATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTc                                                   >  1:2161366/1‑93 (MQ=255)

                                       aTTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCg                                                  <  2:1345383/93‑1 (MQ=255)

                                        ttCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAACACGACAAACCGCTCgg                                                 >  2:4442504/1‑93 (MQ=255)

                                        ttCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCgt                                                 <  2:814314/93‑2 (MQ=255)

                                        ttCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCgg                                                 >  1:224841/1‑93 (MQ=255)

                                        ttCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCgg                                                 <  1:4625096/93‑1 (MQ=255)

                                        ttCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCgg                                                 >  2:362827/1‑93 (MQ=255)

                                          cTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGa                                               <  2:679227/93‑1 (MQ=255)

                                          cTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGa                                               >  1:473611/1‑93 (MQ=255)

                                           tGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGaa                                              <  2:1013106/93‑1 (MQ=255)

                                           tGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGaa                                              <  1:2153124/93‑1 (MQ=255)

                                            tgtgtTAGCGGGTAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAATCGACAAACCGGTCGGGAAt                                             <  2:4545710/92‑1 (MQ=255)

                                            ggTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAAt                                             <  2:3926289/93‑1 (MQ=255)

                                            ggTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAAt                                             <  1:4571314/93‑1 (MQ=255)

                                            ggTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAAt                                             <  2:3136149/93‑1 (MQ=255)

                                             gtgtTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATg                                            <  1:371011/93‑1 (MQ=255)

                                               gtTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCt                                          <  1:290784/93‑1 (MQ=255)

                                                  aGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAtt                                       <  1:1442770/93‑1 (MQ=255)

                                                  aGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAtt                                       <  1:4322643/93‑1 (MQ=255)

                                                  aGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAtt                                       <  2:192305/93‑1 (MQ=255)

                                                   gCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAttt                                      >  1:2076545/1‑93 (MQ=255)

                                                   gCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAttt                                      <  2:4627456/93‑1 (MQ=255)

                                                     ggggAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTg                                    >  2:2252998/1‑93 (MQ=255)

                                                      gggAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGt                                   >  1:2128318/1‑93 (MQ=255)

                                                      gggAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGt                                   <  2:280439/93‑1 (MQ=255)

                                                      gggAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCACTCGGGAATGCTATTTTGt                                   >  2:379060/1‑93 (MQ=255)

                                                         aaCAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCt                                <  2:2691228/93‑1 (MQ=255)

                                                          aCAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTc                               <  1:177777/93‑1 (MQ=255)

                                                             aaTTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGa                            <  1:4172879/93‑1 (MQ=255)

                                                             aaTTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGa                            >  1:2639888/1‑93 (MQ=255)

                                                                 gTCGGTATCAGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGAccgc                        <  2:1672203/93‑1 (MQ=255)

                                                                    ggTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCgg                     <  1:2782560/93‑1 (MQ=255)

                                                                     gTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCggg                    >  1:80535/1‑93 (MQ=255)

                                                                       aTCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCg                  >  2:4201714/1‑93 (MQ=255)

                                                                       aTCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCg                  <  1:3946294/93‑1 (MQ=255)

                                                                        tCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGa                 >  1:2597950/1‑93 (MQ=255)

                                                                          ggCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATg               <  2:1502481/93‑1 (MQ=255)

                                                                             aTCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGaa            >  2:4005558/1‑93 (MQ=255)

                                                                             aTCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGaa            <  2:2019321/93‑1 (MQ=255)

                                                                             aTCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGaa            >  2:1894337/1‑93 (MQ=255)

                                                                              tctcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGaaa           >  2:3321906/1‑93 (MQ=255)

                                                                              tctcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGaaa           <  1:251780/93‑1 (MQ=255)

                                                                               ctcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAt          <  2:1413239/93‑1 (MQ=255)

                                                                               ctcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAt          <  1:3503483/93‑1 (MQ=255)

                                                                               ctcCGCACAGGAAAAGGGCTTGTTTTTGCGGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAt          <  2:3799309/93‑1 (MQ=255)

                                                                                 ccGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcg        >  1:3766587/1‑93 (MQ=255)

                                                                                   gCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgtt      <  2:2189484/93‑1 (MQ=255)

                                                                                   gCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgtt      <  1:984976/93‑1 (MQ=255)

                                                                                     acaGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgttcg    >  2:951978/1‑93 (MQ=255)

                                                                                       aGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAATCGTtcgtc  <  2:3699878/93‑1 (MQ=255)

GGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGGAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAATCGTTCGTC  >  NC_000913/3744868‑3745047

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