BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,585,577	G→A	intergenic (+126/‑38)	narQ → / → acrD	NarQ/AcrD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,585,577	0	G	A	97.1%	581.2 / 1.9	174	intergenic (+126/‑38)	narQ/acrD	NarQ/AcrD
Reads supporting (aligned to +/- strand): ref base G (0/5); new base A (67/102); total (67/107)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

Only 100 of 179 total aligned reads displayed.

GTTCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAGCCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATCC  >  NC_000913/2585487‑2585666

 ttCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCt                                                                                        <  2:3861100/93‑1 (MQ=255)

 ttCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCt                                                                                        >  1:4542503/1‑93 (MQ=255)

  tCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAGCCTa                                                                                       <  2:4107094/93‑1 (MQ=255)

   cAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTAc                                                                                      >  1:1101504/1‑93 (MQ=255)

    aGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACa                                                                                     >  2:2049139/1‑93 (MQ=255)

    aGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACa                                                                                     >  1:1049489/1‑93 (MQ=255)

      tataAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAAc                                                                                   <  2:1939472/93‑1 (MQ=255)

      tataAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAAc                                                                                   <  2:1723918/93‑1 (MQ=255)

        taAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGa                                                                                 >  2:3345596/1‑93 (MQ=255)

         aaaaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGAAAGTCAAACCTACAACGAt                                                                                <  1:781389/93‑1 (MQ=255)

          aaaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATa                                                                               >  1:4165704/1‑93 (MQ=255)

          aaaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATa                                                                               >  2:913583/1‑93 (MQ=255)

           aaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATAc                                                                              <  1:1145423/93‑1 (MQ=255)

           aaGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATAc                                                                              <  2:276412/93‑1 (MQ=255)

           aaGGGCAGGATAATTTACATTACCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTTCGGTGAAGCAAGTCAAACCTACAACGATAc                                                                              <  1:1832508/93‑1 (MQ=255)

            aGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACg                                                                             >  1:722179/1‑93 (MQ=255)

             gggCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGc                                                                            <  2:1923710/93‑1 (MQ=255)

              ggCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCTCTACAGTGAAGCAAGTCAAACCTACAACGATACGCa                                                                           <  1:1540147/93‑1 (MQ=255)

              ggCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCa                                                                           <  2:1016996/93‑1 (MQ=255)

               gCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAg                                                                          >  1:2023454/1‑93 (MQ=255)

               gCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACATTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAg                                                                          <  2:4058872/93‑1 (MQ=255)

                cATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGa                                                                         <  1:1882260/93‑1 (MQ=255)

                 aTGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGaa                                                                        >  1:164590/1‑93 (MQ=255)

                 aTGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGaa                                                                        >  1:4234876/1‑93 (MQ=255)

                   gATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAac                                                                      >  2:3947630/1‑93 (MQ=255)

                     tAATTTATATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAacac                                                                    <  1:4320733/93‑1 (MQ=255)

                     tAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAacac                                                                    <  1:3357331/93‑1 (MQ=255)

                     tAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAacac                                                                    >  2:7762/1‑93 (MQ=255)

                     tAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAacac                                                                    <  1:2578008/93‑1 (MQ=255)

                     tAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAacac                                                                    >  2:823593/1‑93 (MQ=255)

                     tAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAacac                                                                    >  1:2739898/1‑93 (MQ=255)

                      aaTTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACg                                                                   >  2:1048534/1‑93 (MQ=255)

                       aTTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGa                                                                  <  1:3514901/93‑1 (MQ=255)

                       aTTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGa                                                                  <  2:580339/93‑1 (MQ=255)

                        tttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAg                                                                 <  1:2038509/93‑1 (MQ=255)

                        tttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAg                                                                 <  1:630362/93‑1 (MQ=255)

                        tttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAg                                                                 >  1:516146/1‑93 (MQ=255)

                        tttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGCACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAg                                                                 <  1:3772813/93‑1 (MQ=255)

                         ttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAgg                                                                <  2:518574/93‑1 (MQ=255)

                         ttACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAgg                                                                <  1:3631031/93‑1 (MQ=255)

                          tACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGt                                                               <  2:228989/93‑1 (MQ=255)

                             aTTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCct                                                            <  2:3138521/93‑1 (MQ=255)

                                aaCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCtt                                                         >  2:2986336/1‑93 (MQ=255)

                                aaCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCtt                                                         <  2:2217516/93‑1 (MQ=255)

                                 aCTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCttt                                                        >  2:2737542/1‑93 (MQ=255)

                                   tCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTa                                                      >  1:2604441/1‑93 (MQ=255)

                                   tCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTa                                                      >  2:2363393/1‑93 (MQ=255)

                                        tttttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGc                                                 >  2:626073/1‑93 (MQ=255)

                                          tttCTCCGCGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAGCCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGa                                               <  1:4504927/93‑1 (MQ=255)

                                          tttCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGa                                               >  1:1282112/1‑93 (MQ=255)

                                            tctcCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAAt                                             <  1:650963/93‑1 (MQ=255)

                                              tcCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAAttt                                           <  1:2677066/93‑1 (MQ=255)

                                               ccACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGGAGAAACACGAGGTCCTCTTTTAATTGAGACtttt                                          >  2:188356/1‑92 (MQ=255)

                                               ccACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTc                                          >  2:188355/1‑93 (MQ=255)

                                                cACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCt                                         <  2:3697571/93‑1 (MQ=255)

                                                cACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCt                                         >  1:2563681/1‑93 (MQ=255)

                                                cACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCt                                         >  1:81311/1‑93 (MQ=255)

                                                 ccGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCtt                                        <  2:1304161/92‑1 (MQ=255)

                                                 aCGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCtt                                        <  1:2378469/93‑1 (MQ=255)

                                                 aCGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCtt                                        <  2:667749/93‑1 (MQ=255)

                                                 aCGATGGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACTCGAGGTCCTCTTTTAATGGCGAATTTCtt                                        <  1:675995/93‑1 (MQ=255)

                                                    aTTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTAt                                     <  2:4452147/93‑1 (MQ=255)

                                                    aTTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTAt                                     <  2:1386388/93‑1 (MQ=255)

                                                     ttGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTAtt                                    <  1:327161/93‑1 (MQ=255)

                                                     ttGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTAtt                                    <  1:1787821/93‑1 (MQ=255)

                                                      tGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTg                                   >  1:456738/1‑93 (MQ=255)

                                                       ggCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTgt                                  >  2:1544622/1‑92 (MQ=255)

                                                       ggCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGa                                  >  1:1039895/1‑93 (MQ=255)

                                                       ggCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGGCCTCTTTTAATGGCGAATTTCTTTATTGa                                  >  1:325942/1‑93 (MQ=255)

                                                        gCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGAt                                 <  2:3511368/93‑1 (MQ=255)

                                                        gCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGAt                                 <  1:3364951/93‑1 (MQ=255)

                                                        gCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGAt                                 <  1:1425002/93‑1 (MQ=255)

                                                        gCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGAt                                 <  1:998257/93‑1 (MQ=255)

                                                         cTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATc                                <  2:1474521/93‑1 (MQ=255)

                                                             tACCTTGGCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGccc                            <  1:2493578/93‑1 (MQ=255)

                                                             tACCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGccc                            <  2:953618/93‑1 (MQ=255)

                                                              aCCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGcccc                           <  1:368400/93‑1 (MQ=255)

                                                              aCCTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGGCCTCTTTTAATGGCGAATTTCTTTATTGATCGcccc                           >  2:516922/1‑93 (MQ=255)

                                                               ccTTGCCGCTACAGTGAAGCAAGTCAAACCTACAGCGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCa                          <  1:1239088/93‑1 (MQ=255)

                                                               ccTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCa                          <  1:1708434/93‑1 (MQ=255)

                                                                cTTGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCAt                         >  1:873947/1‑93 (MQ=255)

                                                                 ttGCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCAtt                        <  2:3217030/93‑1 (MQ=255)

                                                                   gCCGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCAtttt                      <  1:442126/93‑1 (MQ=255)

                                                                     cGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCAtttttc                    >  1:3766354/1‑92 (MQ=255)

                                                                     cGCTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTg                    <  1:4110816/93‑1 (MQ=255)

                                                                       cTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGcc                  <  2:2931488/93‑1 (MQ=255)

                                                                       cTACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGcc                  <  2:4479867/93‑1 (MQ=255)

                                                                        tACAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCt                 <  2:3676611/93‑1 (MQ=255)

                                                                         aCAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTg                <  1:1044571/93‑1 (MQ=255)

                                                                         aCAGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTg                <  2:3835291/93‑1 (MQ=255)

                                                                           aGTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTggg              <  2:318684/93‑1 (MQ=255)

                                                                            gTGAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGt             >  2:2415448/1‑93 (MQ=255)

                                                                              gAAGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGc           <  2:3407213/93‑1 (MQ=255)

                                                                                aGCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTg         <  1:2882405/93‑1 (MQ=255)

                                                                                 gCAAGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTgg        <  1:1214448/93‑1 (MQ=255)

                                                                                   aaGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGACCCATTTTTGCCTGGGTGCTGGCa      <  1:1715774/93‑1 (MQ=255)

                                                                                    aGTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCaa     <  1:4257983/93‑1 (MQ=255)

                                                                                     gTCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAAt    >  2:1446633/1‑93 (MQ=255)

                                                                                      tCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATc   <  2:454390/93‑1 (MQ=255)

                                                                                      tCAAACCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATc   <  2:657584/93‑1 (MQ=255)

GTTCAGTATAAAAGGGCATGATAATTTACATTAACTCCTTTTTTTCTCCACGATTGGCTCGTACCTTGCCGCTACAGTGAAGCAAGTCAAGCCTACAACGATACGCAGAAACACGAGGTCCTCTTTTAATGGCGAATTTCTTTATTGATCGCCCCATTTTTGCCTGGGTGCTGGCAATCC  >  NC_000913/2585487‑2585666

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