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.3%	1523.6 / ‑5.2	427	G82E (GGA→GAA)	lyxK	L‑xylulose kinase
Reads supporting (aligned to +/- strand): ref base G (3/0); new base A (190/233); total (194/233)
Fisher's exact test for biased strand distribution p-value = 9.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 450 total aligned reads displayed.

GGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGGAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAATCGTTCGTC  >  NC_000913/3744868‑3745047

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

   aaTGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAgg                                                                                      <  1:1624093/93‑1 (MQ=255)

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

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

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

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

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

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

      gctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCt                                                                                   >  1:5546345/1‑93 (MQ=255)

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

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

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

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

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

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

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

            tGGCTGTCATTCGCGCCTTTCGTACTCATCCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGtttt                                                                             <  2:2540456/93‑1 (MQ=255)

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

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

               cTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGc                                                                          >  2:5031764/1‑93 (MQ=255)

                tGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCt                                                                         <  2:5303409/93‑1 (MQ=255)

                tGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCt                                                                         >  1:1787761/1‑93 (MQ=255)

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

                  tCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTgg                                                                       >  1:7121798/1‑93 (MQ=255)

                     ttCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATa                                                                    >  1:7427872/1‑93 (MQ=255)

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

                        gcgcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATaaaa                                                                 >  2:2188463/1‑93 (MQ=255)

                         cgcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATaaaaa                                                                <  1:6292862/93‑1 (MQ=255)

                         cgcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATaaaaa                                                                <  2:5709589/93‑1 (MQ=255)

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

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

                             cTGCTTACTCATTCTGGTGTTAGCGGGGGACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGAc                                                            >  2:1530233/1‑93 (MQ=255)

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

                             cTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCGCCGCACAGGAAAAGGGCTTGTTTTTGCTGGAGAAAAACGAc                                                            >  2:7030037/1‑93 (MQ=255)

                               gCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACaa                                                          <  1:2082991/93‑1 (MQ=255)

                                cTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACaaa                                                         >  2:3489886/1‑93 (MQ=255)

                                cTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACaaa                                                         >  1:2876856/1‑93 (MQ=255)

                                 ttACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAAc                                                        >  2:6033823/1‑93 (MQ=255)

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

                                 ttACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAAc                                                        >  1:3922884/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                                         tCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCggg                                                <  1:3659713/93‑1 (MQ=255)

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

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

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

                                           tGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAGAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGaa                                              <  1:355885/93‑1 (MQ=255)

                                           tGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGaa                                              >  2:6103291/1‑93 (MQ=255)

                                            ggTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAAt                                             >  1:236336/1‑93 (MQ=255)

                                                ttAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTa                                         <  2:1437644/93‑1 (MQ=255)

                                                ttAGCGGGGAACAAATTGTCGGTATCGACATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTa                                         <  1:4164302/93‑1 (MQ=255)

                                                 tAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAt                                        >  1:2025815/1‑93 (MQ=255)

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

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

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

                                                    cGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAtttt                                     <  2:4706724/93‑1 (MQ=255)

                                                    cGGGCAACAAATTGTCGGTATCGCCATCTCCGCTCAGGAAAAGGACTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAtttt                                     <  2:5451227/93‑1 (MQ=255)

                                                     ggggAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTg                                    <  2:1952649/93‑1 (MQ=255)

                                                       ggAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTc                                  >  1:3602729/1‑93 (MQ=255)

                                                        gAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGGTATTTTGTcc                                 >  2:879768/1‑93 (MQ=255)

                                                        gAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTcc                                 <  2:3893712/93‑1 (MQ=255)

                                                        gAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTcc                                 <  1:7860439/93‑1 (MQ=255)

                                                        gAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTcc                                 >  1:6607532/1‑93 (MQ=255)

                                                          aCAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTc                               >  1:2741776/1‑93 (MQ=255)

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

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

                                                           cAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCg                              <  1:1109774/93‑1 (MQ=255)

                                                           cAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCg                              <  2:2682679/93‑1 (MQ=255)

                                                            aaaTTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATATAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTctg                             >  1:7044244/1‑91 (MQ=255)

                                                              aTTGTCGGTATCGGCAGCTCCGCTCAGGAAAATGGCGTGCTTTTCCTGGATAAAAACTACAAAACCCTCGGGAATGCTATTTTGTCCTCGGAc                           <  2:2496902/93‑1 (MQ=255)

                                                                tGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAAGGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGAccg                         >  2:5549591/1‑93 (MQ=255)

                                                                   cGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGAccgccg                      <  1:2526210/93‑1 (MQ=255)

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

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

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

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

                                                                           gCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATgg              >  1:4013726/1‑93 (MQ=255)

                                                                           gCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATgg              >  2:6724572/1‑93 (MQ=255)

                                                                             aTCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACGAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGaa            <  1:7563466/93‑1 (MQ=255)

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

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

                                                                                tcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtc         >  1:2949398/1‑93 (MQ=255)

                                                                                  cGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgt       <  2:7180277/93‑1 (MQ=255)

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

                                                                                   gCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgtt      >  2:3128393/1‑93 (MQ=255)

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

                                                                                   gCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgtt      >  1:4630718/1‑93 (MQ=255)

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

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

                                                                                    cacaGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgttc     <  2:5951181/93‑1 (MQ=255)

                                                                                      caGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgttcgt   <  2:4423538/93‑1 (MQ=255)

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

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

                                                                                       aGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAATCGTtcgtc  <  2:1150625/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` `< 28 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`