BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,643,412	T→C	E82E (GAA→GAG)	relE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,643,412	0	T	C	99.6%	855.7 / NA	233	E82E (GAA→GAG)	relE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (104/128); minor base A (1/0); total (105/128)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

AAGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGTTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAAT  >  NC_000913/1643322‑1643501

 aGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGctct                                                                                        <  1:3462737/93‑1 (MQ=255)

   cAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGctctct                                                                                      >  1:1458007/1‑93 (MQ=255)

    aGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCtt                                                                                     >  1:5873739/1‑93 (MQ=255)

     gAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCttt                                                                                    <  1:1153971/93‑1 (MQ=255)

     gAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCttt                                                                                    >  2:6306430/1‑93 (MQ=255)

       aGTGTCACCTTCGGTGCGACACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTc                                                                                  >  1:6482630/1‑93 (MQ=255)

           tCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCaa                                                                              <  2:5213529/93‑1 (MQ=255)

            cACCTTCGGTGCGAAACAGAGATGGCATGCTTTGGTTCAGAGAATGCGTTTGACAGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAAc                                                                             >  1:4296584/1‑93 (MQ=255)

              ccTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAg                                                                           <  2:5296461/93‑1 (MQ=255)

              ccTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAg                                                                           <  1:3807754/93‑1 (MQ=255)

               cTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGa                                                                          <  1:17712/93‑1 (MQ=255)

                ttCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGaa                                                                         <  1:1560791/93‑1 (MQ=255)

                   ggTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATc                                                                      <  2:710601/93‑1 (MQ=255)

                   ggTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATc                                                                      >  1:4551131/1‑93 (MQ=255)

                        gAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGaa                                                                 >  2:3622170/1‑93 (MQ=255)

                         aaaCAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGaaa                                                                >  1:4058356/1‑93 (MQ=255)

                          aaCAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGaaaa                                                               <  2:93566/93‑1 (MQ=255)

                          aaCAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGaaaa                                                               <  2:4716121/93‑1 (MQ=255)

                            cAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACg                                                             <  1:6629840/93‑1 (MQ=255)

                             agagATGTCATGCTTTGGTTCAGAGCATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGa                                                            >  1:2036387/1‑93 (MQ=255)

                             agagATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGa                                                            <  2:963599/93‑1 (MQ=255)

                              gagaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAACACGAc                                                           >  1:1249762/1‑93 (MQ=255)

                               agaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACa                                                          >  2:1306816/1‑93 (MQ=255)

                               agaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACa                                                          <  1:238963/93‑1 (MQ=255)

                                gaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACaa                                                         <  1:4881462/93‑1 (MQ=255)

                                  tGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACt                                                       >  2:3652087/1‑93 (MQ=255)

                                  tGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACt                                                       >  1:2570140/1‑93 (MQ=255)

                                   gTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACtt                                                      >  2:4350149/1‑93 (MQ=255)

                                      aTGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTtct                                                   <  2:906905/93‑1 (MQ=255)

                                       tGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTtctc                                                  <  2:3196714/93‑1 (MQ=255)

                                         cTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGt                                                >  1:1374889/1‑93 (MQ=255)

                                          tttGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTc                                               >  1:3005578/1‑93 (MQ=255)

                                          tttGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTc                                               <  2:4296984/93‑1 (MQ=255)

                                          tttGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTc                                               >  1:1879619/1‑93 (MQ=255)

                                           ttGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATGCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCt                                              <  2:1882156/93‑1 (MQ=255)

                                           ttGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCt                                              <  2:2595667/93‑1 (MQ=255)

                                           ttGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCt                                              <  2:4214804/93‑1 (MQ=255)

                                            tGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCta                                             >  1:2982508/1‑93 (MQ=255)

                                            tGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCta                                             <  2:6114698/93‑1 (MQ=255)

                                             ggTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCtat                                            <  2:821912/93‑1 (MQ=255)

                                             ggTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCtat                                            <  1:874924/93‑1 (MQ=255)

                                                 cAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAAcc                                        >  2:4045774/1‑93 (MQ=255)

                                                  agagAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCt                                       <  1:3346232/93‑1 (MQ=255)

                                                  agagAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCt                                       <  1:5762852/93‑1 (MQ=255)

                                                   gagaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTg                                      <  2:2466088/93‑1 (MQ=255)

                                                   gagaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTg                                      >  1:2127331/1‑93 (MQ=255)

                                                    agaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTgg                                     <  1:6590328/93‑1 (MQ=255)

                                                    agaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTgg                                     <  2:4063374/93‑1 (MQ=255)

                                                    agaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTgg                                     <  1:557945/93‑1 (MQ=255)

                                                    agaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTgg                                     <  1:402754/93‑1 (MQ=255)

                                                     gaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGt                                    >  2:6337201/1‑93 (MQ=255)

                                                      aaTGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGta                                   <  1:4217258/93‑1 (MQ=255)

                                                       aTGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGtat                                  >  1:5230213/1‑93 (MQ=255)

                                                        tGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGtata                                 >  1:3744222/1‑93 (MQ=255)

                                                        tGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGtata                                 >  2:2106991/1‑93 (MQ=255)

                                                        tGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTGCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGtata                                 <  2:2369713/93‑1 (MQ=255)

                                                         gCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATAc                                <  2:3367715/93‑1 (MQ=255)

                                                           gTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACaa                              <  2:1161873/93‑1 (MQ=255)

                                                           gTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACaa                              <  2:1076370/93‑1 (MQ=255)

                                                           gTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACaa                              <  1:506829/93‑1 (MQ=255)

                                                            tttGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAg                             <  2:6504292/93‑1 (MQ=255)

                                                            tttGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAg                             >  2:5618827/1‑93 (MQ=255)

                                                            tttGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAg                             <  1:5797900/93‑1 (MQ=255)

                                                             ttGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAgg                            <  1:367040/93‑1 (MQ=255)

                                                             ttGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAgg                            <  2:2196498/93‑1 (MQ=255)

                                                              tGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGc                           <  1:1086767/93‑1 (MQ=255)

                                                              tGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGc                           >  2:6091130/1‑93 (MQ=255)

                                                              tGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGc                           >  1:4635248/1‑93 (MQ=255)

                                                              tGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGc                           <  2:4245448/93‑1 (MQ=255)

                                                              tGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGc                           <  2:1674046/93‑1 (MQ=255)

                                                                aCCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGa                         <  2:1052029/93‑1 (MQ=255)

                                                                aCCGCCTCGCTATATACTTCCGAGCGCTCACTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGa                         <  2:281940/93‑1 (MQ=255)

                                                                 ccGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGAt                        <  1:1401913/93‑1 (MQ=255)

                                                                  cGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATa                       <  2:4669365/93‑1 (MQ=255)

                                                                    ccTCTCTATATACTTCCGACCGCTCCCTTCTCCCAACAGAACTCACGAAAACGACACCTTTCTCGTCTATAACCTGGTATACAAGGCGATAGc                     <  2:3465873/93‑1 (MQ=255)

                                                                     cTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGcc                    >  2:5188381/1‑93 (MQ=255)

                                                                      tCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCt                   >  2:2806670/1‑93 (MQ=255)

                                                                      tCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCt                   >  2:5991775/1‑93 (MQ=255)

                                                                      tCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCt                   <  2:4953455/93‑1 (MQ=255)

                                                                      tCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCt                   >  2:1602226/1‑93 (MQ=255)

                                                                         cTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGaa                >  1:3967944/1‑93 (MQ=255)

                                                                         atatataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGaa                <  2:2881598/92‑1 (MQ=255)

                                                                          tatataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAg               >  2:6275061/1‑93 (MQ=255)

                                                                          tatataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGCTATACAAGTCGATAGCCTGAAg               >  1:2520023/1‑93 (MQ=255)

                                                                           atataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGa              >  1:5890685/1‑93 (MQ=255)

                                                                           atataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGa              >  2:3569627/1‑93 (MQ=255)

                                                                             ataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGAcc            <  2:1699160/93‑1 (MQ=255)

                                                                              taCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCg           >  2:4454421/1‑93 (MQ=255)

                                                                              taCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCg           <  1:2923622/93‑1 (MQ=255)

                                                                              taCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCg           <  1:2144166/93‑1 (MQ=255)

                                                                               aCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAATCGACACCTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCgg          <  1:3614005/93‑1 (MQ=255)

                                                                                cTTCCGAGCGCTCTCTTTCCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGa         <  2:5805428/93‑1 (MQ=255)

                                                                                 ttCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAg        >  1:5414318/1‑93 (MQ=255)

                                                                                 ttCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAg        <  2:225400/93‑1 (MQ=255)

                                                                                 ttCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAg        >  1:1211892/1‑93 (MQ=255)

                                                                                   ccGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCt      <  1:1825461/93‑1 (MQ=255)

                                                                                   ccGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCt      <  1:3552508/93‑1 (MQ=255)

                                                                                   ccGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCt      <  1:74140/93‑1 (MQ=255)

                                                                                       gcgcTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAAt  <  1:2717554/93‑1 (MQ=255)

                                                                                       gcgcTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAAt  >  2:3270603/1‑93 (MQ=255)

AAGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGTTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAAT  >  NC_000913/1643322‑1643501

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