BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,555	(G)_5→6	coding (55/618 nt)	b3646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,550	1	.	G	94.9%	664.5 / 26.8	214	coding (50/618 nt)	b3646	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (5/6); new base G (74/129); total (79/135)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-02

Only 100 of 247 total aligned reads displayed.

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

cacaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAc                                                                                          <  2:4570330/93‑1 (MQ=255)

cacaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAc                                                                                          <  1:4529771/93‑1 (MQ=255)

 acaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑g                                                                                        >  2:1874785/1‑93 (MQ=255)

   aaaTCGCTCCCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg                                                                                      <  2:4209098/93‑1 (MQ=255)

   aaaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg                                                                                      <  1:4476197/93‑1 (MQ=255)

    aaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑gggg                                                                                     <  2:4711387/93‑1 (MQ=255)

     aTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggg                                                                                    >  2:2024356/1‑93 (MQ=255)

      tCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑gggggg                                                                                   >  1:3915258/1‑92 (MQ=255)

        gCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggggcg                                                                                 <  2:4342033/93‑4 (MQ=255)

         cTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggggcgc                                                                                >  1:1526356/1‑89 (MQ=255)

           aCCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTg                                                                               >  2:3874426/1‑93 (MQ=255)

            cccTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTgg                                                                              <  1:3643249/93‑1 (MQ=255)

            cccTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTggc                                                                             <  1:3730490/93‑1 (MQ=255)

               tGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTggcgg                                                                           <  2:2094748/93‑1 (MQ=255)

                gCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTggcggc                                                                          <  1:2640521/93‑1 (MQ=255)

                 ccAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGcc                                                                         >  2:2403216/1‑93 (MQ=255)

                     acagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgac                                                                     >  2:1073184/1‑93 (MQ=255)

                     acagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgac                                                                     <  2:2712409/93‑1 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacgac                                                                  >  1:4606839/1‑93 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacgac                                                                  >  1:1052165/1‑93 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacgac                                                                  >  1:4317221/1‑93 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacgac                                                                  <  2:3467725/93‑1 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacgac                                                                  >  2:822274/1‑93 (MQ=255)

                         aTTTTTAGGGAGGGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGacgacg                                                                 <  1:2495735/93‑1 (MQ=255)

                         aTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGacgacg                                                                 <  1:4268263/93‑1 (MQ=255)

                         aTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGacgacg                                                                 <  2:566117/93‑1 (MQ=255)

                         aTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGacgacg                                                                 <  2:20065/93‑1 (MQ=255)

                          tttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCAGGACGACgg                                                                <  1:553001/93‑1 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgc                                                               <  2:2840986/93‑1 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgc                                                               <  1:3961998/93‑1 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgc                                                               <  2:2428346/93‑1 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgc                                                               <  1:3370009/93‑1 (MQ=255)

                            tttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgcg                                                              >  1:1924266/1‑93 (MQ=255)

                             ttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgcgc                                                             >  1:2464496/1‑93 (MQ=255)

                             ttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCAGGACGGGGGGCGCTGGCGGCCGGACGACGgcgc                                                             <  1:1184917/93‑1 (MQ=255)

                              tAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCa                                                            <  2:362497/93‑1 (MQ=255)

                              tAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGCACGACGGCGCa                                                            <  1:4667830/93‑1 (MQ=255)

                              tAGGGAGAGAACCATGCGGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTTGCGGCCGGACGACGGCGCa                                                            >  2:4212421/1‑93 (MQ=255)

                               aGGTAGAGAACCATGCTGTTACACATTTTGTATTGGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCAt                                                           <  2:1131787/93‑1 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCAt                                                           >  2:4355614/1‑93 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCAt                                                           <  1:1436045/93‑1 (MQ=255)

                                gggAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATg                                                          <  1:4013110/93‑1 (MQ=255)

                                gggAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATg                                                          <  1:2836852/93‑1 (MQ=255)

                                 ggAGCGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATgg                                                         >  2:3281688/1‑93 (MQ=255)

                                 ggAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATgg                                                         >  2:2920839/1‑93 (MQ=255)

                                   agagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGAt                                                       <  2:1168329/93‑1 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    >  2:4516097/1‑93 (MQ=255)

                                       aaCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACAt                                                   >  2:65017/1‑93 (MQ=255)

                                       aaCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACAt                                                   <  2:1349400/93‑1 (MQ=255)

                                        aCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACAtt                                                  <  1:316278/93‑1 (MQ=255)

                                            tGCTGTTGCACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGc                                              <  1:2134170/93‑1 (MQ=255)

                                            tGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGc                                              >  1:77349/1‑93 (MQ=255)

                                               tGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTa                                           <  1:3320771/93‑1 (MQ=255)

                                                 ttACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAat                                         <  1:1323248/93‑1 (MQ=255)

                                                 ttACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAat                                         <  1:124442/93‑1 (MQ=255)

                                                   acacATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatta                                       <  1:4725495/93‑1 (MQ=255)

                                                     acaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatt                                     >  2:4261447/1‑93 (MQ=255)

                                                     acaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatt                                     >  1:4267837/1‑93 (MQ=255)

                                                     acaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatta                                    <  1:3925639/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   <  1:4454755/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   <  2:3246619/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   >  1:1927458/1‑93 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  2:126864/93‑1 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  2:2151508/93‑1 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  1:2740571/93‑1 (MQ=255)

                                                        ttttGTATTTGGTGGGCATTACGGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  1:2115081/93‑1 (MQ=255)

                                                         tttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTg                                 <  1:56050/93‑1 (MQ=255)

                                                         tttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTg                                 <  1:1490731/93‑1 (MQ=255)

                                                         tttGTATTGGGTTGGCATTACTGCCGAAGCCATGCCGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTg                                 <  1:2990788/93‑1 (MQ=255)

                                                          ttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGc                                <  1:3830698/93‑1 (MQ=255)

                                                          ttGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCg                               >  1:744185/1‑93 (MQ=255)

                                                           tGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCg                               >  1:2998393/1‑93 (MQ=255)

                                                            gTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACCTTTGGCGTAATTATTATTGCGa                              >  2:1769579/1‑93 (MQ=255)

                                                            gTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGa                              <  1:326895/93‑1 (MQ=255)

                                                            gTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGa                              <  2:3569627/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             >  1:4235883/1‑93 (MQ=255)

                                                              ttttGATTGGCATTACTGCAGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAcc                            <  2:3151396/92‑1 (MQ=255)

                                                              aTTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAcc                            >  2:3644182/1‑93 (MQ=255)

                                                               tttGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGc                          >  1:4414619/1‑93 (MQ=255)

                                                                ttggttgcCATTACTGCCGAAGCCATGACGGGGGGCTCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGc                          <  2:756084/93‑1 (MQ=255)

                                                                 tggttggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGcc                         <  1:1038677/93‑1 (MQ=255)

                                                                 tggttggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGcc                         <  2:2536102/93‑1 (MQ=255)

                                                                   gttggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCAc                       <  2:4458112/93‑1 (MQ=255)

                                                                      ggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGc                    <  2:2230811/93‑1 (MQ=255)

                                                                        cATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  1:3791678/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  2:3164331/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  1:1840872/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  1:2488692/93‑1 (MQ=255)

                                                                          ttACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAtt                <  2:4643488/93‑1 (MQ=255)

                                                                          ttACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAtt                >  1:4121855/1‑93 (MQ=255)

                                                                           tACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTg               >  2:2630681/1‑93 (MQ=255)

                                                                             cTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGc             <  1:2586291/93‑1 (MQ=255)

                                                                               gCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCgg           >  1:3531219/1‑93 (MQ=255)

                                                                               gCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACAACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCgg           <  1:1129688/93‑1 (MQ=255)

                                                                                ccGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGa          <  2:1623721/93‑1 (MQ=255)

                                                                                  gAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAgg        >  2:1409380/1‑93 (MQ=255)

                                                                                   aaGACATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAggg       <  2:1960361/93‑1 (MQ=255)

                                                                                     gCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTc     >  1:1791420/1‑93 (MQ=255)

                                                                                      ccATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCa    <  1:532062/93‑1 (MQ=255)

                                                                                        aTGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGt  >  2:2555842/1‑93 (MQ=255)

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

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