BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,086,377	T→C	E258E (GAA→GAG)	b2016 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,086,377	0	T	C	85.2%	596.7 / 94.3	256	E258E (GAA→GAG)	b2016	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (17/21); new base C (108/110); total (125/131)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

Only 100 of 257 total aligned reads displayed.

CCCCTTGCAGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAGTAAACGGGC  >  NC_000913/2086287‑2086469

ccccTTGCAGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTc                                                                                            <  1:6885200/93‑1 (MQ=255)

ccccTTGCAGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTc                                                                                            <  2:5754599/93‑1 (MQ=255)

  ccTTGCAGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCAt                                                                                          >  2:2298176/1‑93 (MQ=255)

   cTTGCAGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCAtt                                                                                         <  2:2620708/93‑1 (MQ=255)

    ttGCAGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTa                                                                                        >  2:5769989/1‑93 (MQ=255)

     tGCAGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTAc                                                                                       >  2:6494597/1‑93 (MQ=255)

        aGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACaaa                                                                                    <  2:5450024/93‑1 (MQ=255)

        aGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGGATTCAAATCCCAGCTCATTACaaa                                                                                    >  1:1485448/1‑93 (MQ=255)

           gCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATcc                                                                                 <  1:5292210/93‑1 (MQ=255)

             aTGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGa                                                                               <  2:4990502/93‑1 (MQ=255)

              tGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGAc                                                                              >  1:6497136/1‑93 (MQ=255)

              tGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGAc                                                                              <  1:6337881/93‑1 (MQ=255)

              tGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGAc                                                                              >  1:682085/1‑93 (MQ=255)

                cggcgTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTg                                                                            <  1:1976026/93‑1 (MQ=255)

                 ggcgTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGc                                                                           >  1:4278642/1‑93 (MQ=255)

                  gcgTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGcc                                                                          >  2:68595/1‑93 (MQ=255)

                     tCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATc                                                                       >  2:306597/1‑93 (MQ=255)

                     tCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACGGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATc                                                                       >  2:4514770/1‑93 (MQ=255)

                      ccTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCa                                                                      <  1:5090370/93‑1 (MQ=255)

                       ctgctgGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCaa                                                                     >  2:4687951/1‑93 (MQ=255)

                        tgctgGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCaat                                                                    >  1:6567042/1‑93 (MQ=255)

                        tgctgGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACACATCCGACTGCCATCgat                                                                    >  2:3329610/1‑93 (MQ=255)

                        tgctgGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCaat                                                                    >  1:6427339/1‑93 (MQ=255)

                        tgctgGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCaat                                                                    >  2:6266692/1‑93 (MQ=255)

                         gctgGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCaata                                                                   <  2:5299187/93‑1 (MQ=255)

                         gctgGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCaata                                                                   <  1:5869191/93‑1 (MQ=255)

                            ggCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATc                                                                <  1:5637532/93‑1 (MQ=255)

                             gCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCt                                                               <  2:6200941/93‑1 (MQ=255)

                                tGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTAc                                                            >  2:6628775/1‑93 (MQ=255)

                                 gATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTAcc                                                           <  1:3175978/93‑1 (MQ=255)

                                 gATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTAcc                                                           <  1:1066544/93‑1 (MQ=255)

                                 gATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTAcc                                                           >  2:382825/1‑93 (MQ=255)

                                  aTCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCt                                                          <  1:6118720/93‑1 (MQ=255)

                                   tCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCtt                                                         >  1:1424655/1‑93 (MQ=255)

                                   tCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCtt                                                         <  2:4099829/93‑1 (MQ=255)

                                   tCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCtt                                                         <  2:6117144/93‑1 (MQ=255)

                                   tCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCtt                                                         <  2:3253303/93‑1 (MQ=255)

                                    cGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCttt                                                        <  1:830755/93‑1 (MQ=255)

                                        tACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCg                                                    >  1:5811593/1‑93 (MQ=255)

                                         aCTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGc                                                   <  1:553346/93‑1 (MQ=255)

                                         aCTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGc                                                   >  2:5173974/1‑93 (MQ=255)

                                              aGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTcc                                              <  1:7270811/93‑1 (MQ=255)

                                              aGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTcc                                              <  2:1944503/93‑1 (MQ=255)

                                              aGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTcc                                              >  1:2748956/1‑93 (MQ=255)

                                               gCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCa                                             >  2:6025773/1‑93 (MQ=255)

                                               gCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCa                                             <  1:2400900/93‑1 (MQ=255)

                                               gCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCa                                             >  2:7385789/1‑93 (MQ=255)

                                                  gCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGAc                                          >  1:49027/1‑93 (MQ=255)

                                                  gCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGAc                                          >  1:3436379/1‑93 (MQ=255)

                                                  gCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGAc                                          >  2:2474817/1‑93 (MQ=255)

                                                    aTTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTa                                        <  2:7502449/93‑1 (MQ=255)

                                                    aTTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTa                                        <  1:2553403/93‑1 (MQ=255)

                                                     ttACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTAt                                       <  1:4592615/93‑1 (MQ=255)

                                                     ttACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTAt                                       >  2:673509/1‑93 (MQ=255)

                                                      tACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTAtt                                      >  2:16613/1‑93 (MQ=255)

                                                      tACCAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTAtt                                      <  2:3418987/93‑1 (MQ=255)

                                                        ccAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTc                                    <  2:3293343/93‑1 (MQ=255)

                                                        ccAGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTc                                    >  2:5905203/1‑93 (MQ=255)

                                                          aGCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTg                                  <  1:6388129/93‑1 (MQ=255)

                                                           gCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGa                                 <  1:576230/93‑1 (MQ=255)

                                                           gCGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGa                                 <  2:130423/93‑1 (MQ=255)

                                                            cGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGaa                                <  1:2085111/93‑1 (MQ=255)

                                                            cGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGaa                                >  1:1610310/1‑93 (MQ=255)

                                                            cGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGaa                                >  1:4342865/1‑93 (MQ=255)

                                                            cGGATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCAAGACTATTTCTGaa                                <  2:6451332/93‑1 (MQ=255)

                                                             ggATCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAAc                               >  2:2408996/1‑93 (MQ=255)

                                                              gATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACt                              <  2:3899888/93‑1 (MQ=255)

                                                               aTCGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTg                             <  1:5991569/93‑1 (MQ=255)

                                                                 cGGGATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCg                           >  2:3657993/1‑93 (MQ=255)

                                                                  gggATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCgg                          >  2:2186333/1‑93 (MQ=255)

                                                                   ggATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGa                         <  1:3804661/93‑1 (MQ=255)

                                                                   ggATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGa                         >  1:370357/1‑93 (MQ=255)

                                                                   ggATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGa                         <  2:5509671/93‑1 (MQ=255)

                                                                    gATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAg                        >  1:3046537/1‑93 (MQ=255)

                                                                    gATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAg                        <  1:3926702/93‑1 (MQ=255)

                                                                    gATACTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAg                        >  1:2718428/1‑93 (MQ=255)

                                                                     aTACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAgg                       <  2:2071864/93‑1 (MQ=255)

                                                                      tACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGt                      >  2:2975175/1‑93 (MQ=255)

                                                                       aCTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGtt                     >  2:2760496/1‑93 (MQ=255)

                                                                        cTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTc                    <  2:506105/93‑1 (MQ=255)

                                                                        cTGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTc                    >  2:468995/1‑93 (MQ=255)

                                                                         tGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTcc                   >  1:5220047/1‑93 (MQ=255)

                                                                         tGGTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTcc                   <  1:3368813/93‑1 (MQ=255)

                                                                          ggTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCa                  <  1:2665919/93‑1 (MQ=255)

                                                                           gTATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAg                 <  2:7325366/93‑1 (MQ=255)

                                                                            tATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGc                >  2:1336310/1‑93 (MQ=255)

                                                                            tATTCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGc                <  2:1109591/93‑1 (MQ=255)

                                                                              ttCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGccc              >  2:3977984/1‑93 (MQ=255)

                                                                              ttCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGccc              <  2:1038448/93‑1 (MQ=255)

                                                                               tCAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGcccc             <  1:81879/93‑1 (MQ=255)

                                                                                cAAATCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCa            <  2:7103976/93‑1 (MQ=255)

                                                                                 aaaTCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAg           <  2:7162163/93‑1 (MQ=255)

                                                                                 aaaTCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAg           <  2:2896128/93‑1 (MQ=255)

                                                                                   aTCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAGTa         >  2:3594738/1‑93 (MQ=255)

                                                                                   aTCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAGTa         <  1:4636322/93‑1 (MQ=255)

                                                                                   aTCCCAGCTCATTACAAATCCGAATGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGCAGGTTCCAGCCCCAGTa         >  2:4853125/1‑93 (MQ=255)

                                                                                    tCCCAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAGTaa        <  2:3357656/93‑1 (MQ=255)

                                                                                      ccAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAGTAAAc      >  1:4709660/1‑93 (MQ=255)

                                                                                      ccAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAGTAAAc      >  1:4689603/1‑93 (MQ=255)

                                                                                      ccAGCTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAGTAAAc      >  1:4506217/1‑93 (MQ=255)

CCCCTTGCAGTGCATGCGGCGTCCTGCTGGCGTGATCGTTTACTCCAGCGGCATTACCAGCGGATCGGGATACTGGTATTCAAATCCCAGTTCATTACAAATCCGACTGCCATCAATAATCTTACCTTTGCCGCTGTCCAGACTATTTCTGAACTGCGGAGGTTCCAGCCCCAGTAAACGGGC  >  NC_000913/2086287‑2086469

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