BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,268,771	G→A	A214V (GCC→GTC)	b3124 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,268,771	0	G	A	99.6%	900.0 / ‑5.3	246	A214V (GCC→GTC)	b3124	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (129/116); total (130/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 253 total aligned reads displayed.

ACATCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGGCTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAATGACGCAATCT  >  NC_000913/3268681‑3268860

aCATCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACt                                                                                         <  2:4130613/93‑1 (MQ=255)

aCATCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACt                                                                                         <  2:2106403/93‑1 (MQ=255)

 cATCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTc                                                                                        >  2:3837933/1‑93 (MQ=255)

 cATCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTc                                                                                        <  2:1492332/93‑1 (MQ=255)

  aTCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTcc                                                                                       >  2:1589225/1‑93 (MQ=255)

   tCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTccc                                                                                      <  2:831966/93‑1 (MQ=255)

   tCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTccc                                                                                      <  2:4349944/93‑1 (MQ=255)

   tCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTccc                                                                                      >  2:2786088/1‑93 (MQ=255)

    cTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCt                                                                                     <  2:1414993/93‑1 (MQ=255)

    cTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCt                                                                                     <  2:3729369/93‑1 (MQ=255)

     tttCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCtt                                                                                    <  2:3205087/93‑1 (MQ=255)

      ttCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCttt                                                                                   <  1:4940585/93‑1 (MQ=255)

      ttCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCttt                                                                                   <  2:4850584/93‑1 (MQ=255)

      ttCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCttt                                                                                   <  1:1597200/93‑1 (MQ=255)

        cacaTCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTtg                                                                                 >  2:2247829/1‑93 (MQ=255)

        cacaTCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTtg                                                                                 >  2:5816493/1‑93 (MQ=255)

         acaTCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTtgt                                                                                <  1:1893810/93‑1 (MQ=255)

          caTCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTtgtg                                                                               <  1:3257014/93‑1 (MQ=255)

           aTCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTgg                                                                              <  2:5623654/93‑1 (MQ=255)

           aTCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTgg                                                                              >  2:6709919/1‑93 (MQ=255)

            tCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTggg                                                                             >  1:2543099/1‑93 (MQ=255)

              aaCATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGcc                                                                           >  1:2110084/1‑93 (MQ=255)

              aaCATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGCGcc                                                                           >  1:2648604/1‑93 (MQ=255)

               aCATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCa                                                                          >  1:1958737/1‑93 (MQ=255)

               aCATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCa                                                                          >  2:287420/1‑93 (MQ=255)

                 aTGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCaaa                                                                        <  1:1387285/93‑1 (MQ=255)

                 aTGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCaaa                                                                        <  1:2318985/93‑1 (MQ=255)

                   gCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAga                                                                      >  1:2721958/1‑93 (MQ=255)

                   gCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAga                                                                      >  1:6189061/1‑93 (MQ=255)

                   gCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAga                                                                      >  1:5691741/1‑93 (MQ=255)

                    cAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAgat                                                                     <  1:438021/93‑1 (MQ=255)

                     aGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAgatg                                                                    >  1:4594155/1‑93 (MQ=255)

                     aGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAgatg                                                                    >  1:5147397/1‑93 (MQ=255)

                      gcgcTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAgatgc                                                                   >  2:6513006/1‑93 (MQ=255)

                      gcgcTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAgatgc                                                                   >  2:2373489/1‑93 (MQ=255)

                           tttttAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAgatgcgcgat                                                              <  2:6082169/93‑1 (MQ=255)

                             tttAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAgatgcgcgatgc                                                            <  1:2793236/93‑1 (MQ=255)

                               tAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAgatgcgcgatgcgc                                                          <  1:1718065/93‑1 (MQ=255)

                               tAATGAACTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAgatgcgcgatgcgc                                                          <  1:547422/93‑1 (MQ=255)

                                  tGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGt                                                       <  2:2500114/93‑1 (MQ=255)

                                  tGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGt                                                       >  2:6610563/1‑93 (MQ=255)

                                  tGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGt                                                       >  1:3442243/1‑93 (MQ=255)

                                  tGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGt                                                       <  1:6666995/93‑1 (MQ=255)

                                  tGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTCACTCCCTTTAGTGGGCCAAAGATGCGCGATGCGCCGt                                                       <  2:3415028/93‑1 (MQ=255)

                                   gACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGtt                                                      >  2:632591/1‑93 (MQ=255)

                                   gACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGtt                                                      <  2:2520488/93‑1 (MQ=255)

                                   gACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGtt                                                      <  2:543997/93‑1 (MQ=255)

                                     ccTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTAt                                                    >  1:969253/1‑93 (MQ=255)

                                       tCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCCAAGATGCGCGATGCGCCGTTATcc                                                  >  1:6588568/1‑92 (MQ=255)

                                         ggCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGcc                                                >  2:526849/1‑93 (MQ=255)

                                         ggCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGcc                                                >  1:1223235/1‑93 (MQ=255)

                                         ggCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGcc                                                <  1:4878191/93‑1 (MQ=255)

                                          gCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGccc                                               >  1:3155175/1‑93 (MQ=255)

                                          gCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGccc                                               <  2:2897935/93‑1 (MQ=255)

                                          gCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGccc                                               <  1:4004269/93‑1 (MQ=255)

                                            aTAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCccac                                             >  2:394540/1‑93 (MQ=255)

                                            aTAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCccac                                             >  1:705744/1‑93 (MQ=255)

                                             tAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCccacc                                            <  2:726552/93‑1 (MQ=255)

                                             tAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCccacc                                            >  2:3121705/1‑93 (MQ=255)

                                             tAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCccacc                                            >  2:3492043/1‑93 (MQ=255)

                                              aGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGGccacca                                           >  1:3859668/1‑93 (MQ=255)

                                               gTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAg                                          >  2:2427149/1‑93 (MQ=255)

                                               gTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAg                                          >  2:6409381/1‑93 (MQ=255)

                                               gTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAg                                          >  2:4687898/1‑93 (MQ=255)

                                                 gagagaGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCg                                        <  1:3267302/93‑1 (MQ=255)

                                                  agagagGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCgg                                       <  2:1047000/93‑1 (MQ=255)

                                                   gagagGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGa                                      <  1:5952955/93‑1 (MQ=255)

                                                    agagGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGAt                                     <  2:1278841/93‑1 (MQ=255)

                                                    agagGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGAt                                     <  2:2720389/93‑1 (MQ=255)

                                                    agagGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGAt                                     <  1:4082306/93‑1 (MQ=255)

                                                    agagGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGAt                                     <  1:3009796/93‑1 (MQ=255)

                                                     gagGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGAtt                                    <  2:4734169/93‑1 (MQ=255)

                                                     gagGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGAtt                                    >  2:4397961/1‑93 (MQ=255)

                                                     gagGTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGAtt                                    >  2:5950494/1‑93 (MQ=255)

                                                       ggTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTgg                                  >  1:1111411/1‑93 (MQ=255)

                                                       ggTTATTGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTgg                                  <  1:4163331/93‑1 (MQ=255)

                                                            ttGTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACa                             <  2:3909001/93‑1 (MQ=255)

                                                              gTCCAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATc                           >  2:6637133/1‑93 (MQ=255)

                                                                ccAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATcac                         <  2:2269377/93‑1 (MQ=255)

                                                                 cAGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATcaca                        >  1:5239874/1‑93 (MQ=255)

                                                                  aGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACaa                       <  1:5503462/93‑1 (MQ=255)

                                                                  aGCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACaa                       <  1:3402941/93‑1 (MQ=255)

                                                                   gCTCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAg                      >  1:6213470/1‑93 (MQ=255)

                                                                     tCAACAATCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCg                    >  2:1982145/1‑93 (MQ=255)

                                                                       aacaaTCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGAc                  <  1:4261531/93‑1 (MQ=255)

                                                                        acaaTCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGAcg                 >  2:5144161/1‑93 (MQ=255)

                                                                        acaaTCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGAcg                 <  2:4885654/93‑1 (MQ=255)

                                                                         caaTCATCGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGAcgc                <  1:3849416/93‑1 (MQ=255)

                                                                           atcatcGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGa              >  1:5210164/1‑93 (MQ=255)

                                                                            tcatcGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGaa             >  1:213441/1‑93 (MQ=255)

                                                                            tcatcGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGaa             <  2:2998721/93‑1 (MQ=255)

                                                                             catcGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAAt            <  1:2875685/93‑1 (MQ=255)

                                                                              atcGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATGACAAGCGACGCGAATg           >  2:4157377/1‑93 (MQ=255)

                                                                                cGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAATGAc         >  1:234028/1‑93 (MQ=255)

                                                                                cGCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAATGAc         <  2:1145403/93‑1 (MQ=255)

                                                                                 gCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAATGACg        >  1:906688/1‑93 (MQ=255)

                                                                                 gCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAATGACg        >  2:5320746/1‑93 (MQ=255)

                                                                                 gCTTCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAATGACg        >  1:1019959/1‑93 (MQ=255)

                                                                                    tCACTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAATGACGCaa     >  2:4204027/1‑93 (MQ=255)

                                                                                      aCTGACTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAATGACGCAATc   <  2:2751546/93‑1 (MQ=255)

ACATCTTTCACATCAACATGCAGCGCTTTTTTAATGACCTCGGCATAGTGAGAGAGGTTATTGTCCAGCTCAACAATCATCGCTTCACTGGCTCCCTTTTGTGGGCCAAAGATGCGCGATGCGCCGTTATCGCCCACCAGCGGATTGGTGACATCACAAGCGACGCGAATGACGCAATCT  >  NC_000913/3268681‑3268860

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