BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,018,487	C→T	R436C (CGT→TGT)	b2879 →	putative proteoglycan

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,018,487	0	C	T	99.7%	1265.2 / ‑5.9	356	R436C (CGT→TGT)	b2879	putative proteoglycan
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (151/203); total (152/204)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 367 total aligned reads displayed.

AATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTT  >  NC_000913/3018397‑3018576

aaTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGgtgt                                                                                         >  1:5270585/1‑93 (MQ=255)

 aTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTc                                                                                        <  2:3226492/93‑1 (MQ=255)

  tGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCg                                                                                       >  2:3040867/1‑93 (MQ=255)

   ggTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCgg                                                                                      >  2:7574745/1‑93 (MQ=255)

    gtgtGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCgga                                                                                     >  1:1508619/1‑93 (MQ=255)

      gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCggatg                                                                                   <  2:7750713/93‑1 (MQ=255)

      gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCggatg                                                                                   <  1:661966/93‑1 (MQ=255)

      gtgATGGTCTATGAAGACCCTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCggatg                                                                                   <  1:1182775/93‑1 (MQ=255)

        gATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCggatgga                                                                                 <  1:4039198/93‑1 (MQ=255)

        gATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCggatgga                                                                                 >  1:5135845/1‑93 (MQ=255)

        gATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCggatgga                                                                                 <  2:3417609/93‑1 (MQ=255)

         aTGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCAAAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCggatggat                                                                                <  2:5711677/93‑1 (MQ=255)

          tGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGgatggatg                                                                               >  2:3578633/1‑93 (MQ=255)

              cTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCt                                                                           <  1:5464635/93‑1 (MQ=255)

                 tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGc                                                                        >  2:6370658/1‑93 (MQ=255)

                    aGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATa                                                                     >  2:5118018/1‑93 (MQ=255)

                     gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATaa                                                                    >  2:256556/1‑93 (MQ=255)

                     gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATaa                                                                    >  2:875049/1‑93 (MQ=255)

                      aCCGTCAGTTTAACTTCGATTGCGATTCCATTTAGGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATaaa                                                                   <  1:1620403/93‑1 (MQ=255)

                        cGTCAGTTTAACTTAGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATg                                                                 <  1:7920625/93‑1 (MQ=255)

                         gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGa                                                                >  2:5420167/1‑93 (MQ=255)

                         gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGa                                                                >  1:2107162/1‑93 (MQ=255)

                          tCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGAc                                                               <  2:5949697/93‑1 (MQ=255)

                          tCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGAc                                                               <  1:2883623/93‑1 (MQ=255)

                            aGTTTAACTTCGATGGCGATCCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAg                                                             <  2:2237592/93‑1 (MQ=255)

                             gTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGa                                                            <  2:7863419/93‑1 (MQ=255)

                             gTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGa                                                            <  1:7307981/93‑1 (MQ=255)

                             gTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGa                                                            <  2:4833781/93‑1 (MQ=255)

                             gTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGa                                                            <  2:5003427/93‑1 (MQ=255)

                             gTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGa                                                            <  2:5331773/93‑1 (MQ=255)

                             gTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGa                                                            >  1:1722009/1‑93 (MQ=255)

                             gTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCAAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGa                                                            <  2:2064970/93‑1 (MQ=255)

                              tttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGAt                                                           <  1:6174708/93‑1 (MQ=255)

                                tAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGc                                                         <  2:134158/93‑1 (MQ=255)

                                  aCTTCGATTGCGTTTCCTTATATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGACTAAAAGAAAGAtaacc                                                       <  2:7031191/93‑5 (MQ=255)

                                  aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGAGCTGGCATAAATGACAGATGccc                                                       >  2:3840702/1‑93 (MQ=255)

                                   cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCct                                                      <  1:2501643/93‑1 (MQ=255)

                                   cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCct                                                      >  1:326009/1‑93 (MQ=255)

                                    ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCctc                                                     <  1:4737200/93‑1 (MQ=255)

                                     tCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCctct                                                    <  1:4396530/93‑1 (MQ=255)

                                       gATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTc                                                  <  1:1020814/93‑1 (MQ=255)

                                        aTTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTcc                                                 >  1:3991809/1‑93 (MQ=255)

                                        aTTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTcc                                                 >  1:156574/1‑93 (MQ=255)

                                         ttGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTccc                                                <  2:6038653/93‑1 (MQ=255)

                                         ttGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTccc                                                >  2:1396719/1‑93 (MQ=255)

                                          tGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCg                                               >  1:7814784/1‑93 (MQ=255)

                                             gATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAg                                            <  2:1284723/93‑1 (MQ=255)

                                                tCCAGTTATGCAGAAGACAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGa                                         <  1:3316763/93‑1 (MQ=255)

                                                  cATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAag                                       <  1:3800277/93‑1 (MQ=255)

                                                  cATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAag                                       <  1:872685/93‑1 (MQ=255)

                                                   aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAagg                                      >  2:2149855/1‑92 (MQ=255)

                                                   aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAaga                                      <  2:6142721/93‑1 (MQ=255)

                                                   aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAaga                                      >  2:2612905/1‑93 (MQ=255)

                                                    tttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAagag                                     <  1:3104640/93‑1 (MQ=255)

                                                      tATGCACAAGGCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAggg                                   <  1:4206168/93‑1 (MQ=255)

                                                      tATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAggg                                   >  2:7702243/1‑93 (MQ=255)

                                                      tATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAggg                                   >  2:6594053/1‑93 (MQ=255)

                                                       aTGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGt                                  <  2:1698586/93‑1 (MQ=255)

                                                       aTGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGt                                  <  2:2403921/93‑1 (MQ=255)

                                                        tGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGtt                                 <  1:4494201/93‑1 (MQ=255)

                                                         gCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTa                                <  1:3480442/93‑1 (MQ=255)

                                                          cacaAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTaa                               <  2:1519386/93‑1 (MQ=255)

                                                           acacGCCAGAAAAGCCGCTGCGAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCATCTTCCCGCAGGGAATAGGGTTaaa                              <  2:6670270/93‑1 (MQ=255)

                                                           acaAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTaaa                              <  1:3040868/93‑1 (MQ=255)

                                                            caAGGCAGAAAAGCCCGTGCCAGTGTGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAt                             <  1:6203570/93‑1 (MQ=255)

                                                             aaGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATa                            <  2:7650786/93‑1 (MQ=255)

                                                             aaGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATa                            >  1:183014/1‑93 (MQ=255)

                                                             aaGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATa                            <  1:7008025/93‑1 (MQ=255)

                                                              aGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATaa                           <  2:1975405/93‑1 (MQ=255)

                                                              aGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATaa                           <  1:490054/93‑1 (MQ=255)

                                                              aGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATaa                           <  1:5756779/93‑1 (MQ=255)

                                                                 cAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGa                        <  2:7334054/93‑1 (MQ=255)

                                                                  aGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGaa                       <  2:42392/93‑1 (MQ=255)

                                                                   gAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGaaa                      >  1:1288805/1‑93 (MQ=255)

                                                                    aaaaGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaag                     >  2:4079343/1‑93 (MQ=255)

                                                                     aaaGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaaga                    <  2:4784224/93‑1 (MQ=255)

                                                                      aaGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaa                   <  2:1003275/93‑1 (MQ=255)

                                                                      aaGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaa                   <  1:2981757/93‑1 (MQ=255)

                                                                      aaGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaa                   <  1:2800997/93‑1 (MQ=255)

                                                                      aaGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaa                   >  2:7115649/1‑93 (MQ=255)

                                                                         ccGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                >  1:2537874/1‑93 (MQ=255)

                                                                         ccGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                <  1:5070523/93‑1 (MQ=255)

                                                                          cGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatg               <  2:1833166/93‑1 (MQ=255)

                                                                           gCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga              <  2:5232721/93‑1 (MQ=255)

                                                                           gCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga              <  1:2189371/93‑1 (MQ=255)

                                                                             tGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt            <  2:4352759/93‑1 (MQ=255)

                                                                             tGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt            <  1:1902552/93‑1 (MQ=255)

                                                                             tGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt            >  1:2243718/1‑93 (MQ=255)

                                                                             tGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt            >  2:6793685/1‑93 (MQ=255)

                                                                             tGCCAGTATGTGGTGTCGGAT‑GATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGGAGAGGGGTAAATAAGGAAAGAAGATGATTg           >  2:3137173/1‑93 (MQ=255)

                                                                               ccAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGa          >  1:3372454/1‑93 (MQ=255)

                                                                                cAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa         >  1:5693067/1‑93 (MQ=255)

                                                                                 aGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAAc        <  2:3702895/93‑1 (MQ=255)

                                                                                  gTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACa       >  1:6323901/1‑93 (MQ=255)

                                                                                  gTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACa       >  1:2691311/1‑93 (MQ=255)

                                                                                  gTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACa       >  1:4313300/1‑93 (MQ=255)

                                                                                     tgtgGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAtt    >  2:554192/1‑93 (MQ=255)

                                                                                     tgtgGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAtt    <  2:4372638/93‑1 (MQ=255)

                                                                                     tgtgGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAtt    <  1:7516183/93‑1 (MQ=255)

                                                                                      gtggtgTCGGATGGATGCGCTGACATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttt   <  2:4664788/93‑1 (MQ=255)

AATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTT  >  NC_000913/3018397‑3018576

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