BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,721,090	C→T	R316* (CGA→TGA)	b1645 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,721,090	0	C	T	99.6%	830.2 / ‑5.7	236	R316* (CGA→TGA)	b1645	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (102/132); total (103/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 240 total aligned reads displayed.

TCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCATACCAGACTCCCA  >  NC_000913/1721000‑1721179

tCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGa                                                                                         <  2:6652707/93‑1 (MQ=255)

 agCAAGCAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGAc                                                                                        <  1:5288970/92‑1 (MQ=255)

   cAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTa                                                                                      <  2:1744231/93‑1 (MQ=255)

   cAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTa                                                                                      >  1:3248304/1‑93 (MQ=255)

    aaaCAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTAc                                                                                     <  1:7033326/93‑1 (MQ=255)

     aaCAGATGTTTACAACGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACCATCGGCACGTCGCCTTCTGGCAGTGACTACg                                                                                    >  2:381555/1‑93 (MQ=255)

      aCAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGt                                                                                   <  1:1807194/93‑1 (MQ=255)

      aCAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGt                                                                                   <  1:6136457/93‑1 (MQ=255)

       cAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGtt                                                                                  >  2:7927097/1‑93 (MQ=255)

       cAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGtt                                                                                  >  2:3052561/1‑93 (MQ=255)

        aGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTa                                                                                 <  1:5370835/93‑1 (MQ=255)

        aGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTa                                                                                 >  1:1519465/1‑93 (MQ=255)

         gATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTAt                                                                                >  1:6724082/1‑93 (MQ=255)

          aTGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTAtt                                                                               >  1:2736158/1‑93 (MQ=255)

           tGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTAttt                                                                              >  2:4452685/1‑93 (MQ=255)

           tGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTAttt                                                                              >  2:4861171/1‑93 (MQ=255)

           tGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTAttt                                                                              >  1:6171870/1‑93 (MQ=255)

           tGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTAttt                                                                              >  1:4893330/1‑93 (MQ=255)

           tGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGCCGCCTTCTGGCAGTGACTACGTTAttt                                                                              >  2:1593846/1‑93 (MQ=255)

             tttACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTAttttt                                                                            <  2:3145736/93‑1 (MQ=255)

              ttACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTg                                                                           <  1:114953/93‑1 (MQ=255)

              ttACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTg                                                                           >  1:2563258/1‑93 (MQ=255)

                 cacCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTgccg                                                                        >  2:474769/1‑93 (MQ=255)

                  acCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTgccgc                                                                       >  1:726102/1‑93 (MQ=255)

                  acCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTgccgc                                                                       >  2:3629916/1‑93 (MQ=255)

                   ccGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTgccgcc                                                                      <  1:2639815/93‑1 (MQ=255)

                     gTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCtt                                                                    <  1:4029640/93‑1 (MQ=255)

                       cGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTa                                                                  <  2:2419920/93‑1 (MQ=255)

                        gCACGTATTATCTCCCCGCTACGCCTGTCCATCGCCGCCGGCTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTAt                                                                 <  2:7412991/93‑1 (MQ=255)

                        gCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTAt                                                                 <  2:5687295/93‑1 (MQ=255)

                        gCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTAt                                                                 >  1:1788671/1‑93 (MQ=255)

                         cACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATc                                                                >  1:5078050/1‑93 (MQ=255)

                         cACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATc                                                                >  1:800849/1‑93 (MQ=255)

                         cACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATc                                                                >  1:5519594/1‑93 (MQ=255)

                          ccGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCt                                                               <  1:7340751/92‑1 (MQ=255)

                          aCGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCt                                                               <  1:3505330/93‑1 (MQ=255)

                          aCGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCt                                                               >  1:5736337/1‑93 (MQ=255)

                           cGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTg                                                              >  1:2320538/1‑93 (MQ=255)

                           cGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTAAGTTATTTTTGCCGCCTTTATCTg                                                              >  2:4303015/1‑93 (MQ=255)

                            gTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGc                                                             >  2:3460301/1‑93 (MQ=255)

                               ttatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCaaa                                                          >  2:5142037/1‑93 (MQ=255)

                               ttatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCaaa                                                          >  2:6975587/1‑93 (MQ=255)

                               ttatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCaaa                                                          >  2:3605894/1‑93 (MQ=255)

                                tatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAag                                                         >  2:2193950/1‑93 (MQ=255)

                                tatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAag                                                         <  1:6512274/93‑1 (MQ=255)

                                tatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAag                                                         >  2:5806214/1‑93 (MQ=255)

                                tatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATCTGCAAag                                                         <  1:3481657/93‑1 (MQ=255)

                                 atCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAagt                                                        <  2:7504093/93‑1 (MQ=255)

                                  tCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAagta                                                       <  1:4772766/93‑1 (MQ=255)

                                   cGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAagtag                                                      <  1:131896/93‑1 (MQ=255)

                                   cGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAagtag                                                      <  2:6238207/93‑1 (MQ=255)

                                            aCGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAAtt                                             <  2:6693252/93‑1 (MQ=255)

                                             cGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTa                                            <  2:2585782/93‑1 (MQ=255)

                                              gCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTAc                                           >  2:299335/1‑93 (MQ=255)

                                               cccGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACa                                          <  1:282070/93‑1 (MQ=255)

                                               cccGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACa                                          <  1:3643401/93‑1 (MQ=255)

                                               cccGACCAACGTCCCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACa                                          <  2:6743551/93‑1 (MQ=255)

                                                 cgcccAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATc                                        <  2:5039749/90‑1 (MQ=255)

                                                 cGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATc                                        <  1:3068242/93‑1 (MQ=255)

                                                 cGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATc                                        >  1:6204287/1‑93 (MQ=255)

                                                    ccAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTc                                     <  2:5241432/93‑1 (MQ=255)

                                                    ccAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTc                                     <  1:1044427/93‑1 (MQ=255)

                                                    ccAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTc                                     <  2:89107/93‑1 (MQ=255)

                                                     cAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCt                                    <  2:5343163/93‑1 (MQ=255)

                                                      caCGTCGCCGGCTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTg                                   <  1:2077946/92‑1 (MQ=255)

                                                      aaCGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTg                                   >  2:2153837/1‑93 (MQ=255)

                                                       aCGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGc                                  <  1:7760501/93‑1 (MQ=255)

                                                         gTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCaa                                >  2:6263036/1‑93 (MQ=255)

                                                         gTCGCCGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCaa                                <  1:6866945/93‑1 (MQ=255)

                                                            cccGACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGc                             <  2:2724615/92‑1 (MQ=255)

                                                               gACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAATGCggg                          >  2:6505687/1‑92 (MQ=255)

                                                               gACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGgt                          >  1:3879551/1‑93 (MQ=255)

                                                               gACTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGgt                          <  2:194185/93‑1 (MQ=255)

                                                                aCTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGgtg                         <  1:5855664/93‑1 (MQ=255)

                                                                 cTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGgtgt                        <  1:7487241/93‑1 (MQ=255)

                                                                 cTATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGgtgt                        >  2:1155026/1‑93 (MQ=255)

                                                                  tATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGTCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTa                       <  1:5913726/93‑1 (MQ=255)

                                                                  tATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTa                       >  1:4690822/1‑93 (MQ=255)

                                                                  tATCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTa                       >  1:7161154/1‑93 (MQ=255)

                                                                    tCGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAga                     <  2:2675409/93‑1 (MQ=255)

                                                                     cGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAgat                    <  2:5956011/93‑1 (MQ=255)

                                                                     cGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAgat                    <  2:6372068/93‑1 (MQ=255)

                                                                     cGGCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAgat                    <  1:417776/93‑1 (MQ=255)

                                                                      ggCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTACAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAgatg                   <  2:605025/93‑1 (MQ=255)

                                                                       gcgcGTCGCCTTCTGGCCGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAgatga                  <  1:2643717/90‑1 (MQ=255)

                                                                       gCACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAgatga                  <  1:5593333/93‑1 (MQ=255)

                                                                        cACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAgatgat                 <  2:2276865/93‑1 (MQ=255)

                                                                        cACGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAgatgat                 <  2:4770824/93‑1 (MQ=255)

                                                                         aCGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATc                <  1:4819371/93‑1 (MQ=255)

                                                                          cGTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCa               >  2:2741692/1‑93 (MQ=255)

                                                                           gTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCCAAGCGGTGTAgctgacaaa              >  2:4938365/1‑89 (MQ=255)

                                                                           gTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCAt              >  1:3002474/1‑93 (MQ=255)

                                                                           gTCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCAt              >  1:5982353/1‑93 (MQ=255)

                                                                            tCGCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCATa             >  2:728323/1‑93 (MQ=255)

                                                                              gCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCATAcc           <  2:1355789/93‑1 (MQ=255)

                                                                              gCCTTCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCATAcc           <  1:7807989/93‑1 (MQ=255)

                                                                                 ttCTGGCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCATACCAGa        <  1:1159186/93‑1 (MQ=255)

                                                                                     ggCAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCATACCAGACTcc    <  1:6437449/93‑1 (MQ=255)

                                                                                       cAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCATACCAGACTCCCa  <  1:4249681/93‑1 (MQ=255)

                                                                                       cAGTGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCATACCAGACTCCCa  <  1:1649891/93‑1 (MQ=255)

TCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGTCAGGAATTACATCGTCTGCAAAGCGGTGTAGATGATCATACCAGACTCCCA  >  NC_000913/1721000‑1721179

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