BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	671,943	T→C	A688A (GCA→GCG)	leuS ←	leucine tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	671,943	0	T	C	99.4%	1116.4 / ‑4.8	346	A688A (GCA→GCG)	leuS	leucine tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (186/158); total (187/159)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 360 total aligned reads displayed.

ACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTG  >  NC_000913/671853‑672032

  ggCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCCGCGCATCAACGTTCAGCGCCg                                                                                       >  1:8617436/1‑93 (MQ=255)

   gCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATTGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGcgccgc                                                                                      <  2:5959842/93‑1 (MQ=255)

   gCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGcgccgc                                                                                      >  1:6138265/1‑93 (MQ=255)

     cGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGcgccgcca                                                                                    >  1:4843393/1‑91 (MQ=255)

      gATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAAc                                                                                   >  2:8821314/1‑93 (MQ=255)

      gATATCATCGGTCACTCGAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCAGCAACGTTCAGCGTCGCAAc                                                                                   >  2:8952664/1‑93 (MQ=255)

        tatCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACAt                                                                                 >  2:7945246/1‑93 (MQ=255)

        tatCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACAt                                                                                 >  1:5183858/1‑93 (MQ=255)

         atcatcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATc                                                                                >  1:4388538/1‑93 (MQ=255)

           catcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCAc                                                                              >  2:8048886/1‑93 (MQ=255)

            atcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCAcc                                                                             >  2:7742274/1‑93 (MQ=255)

            atcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCAcc                                                                             >  2:8564903/1‑93 (MQ=255)

             tcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCt                                                                            <  2:5585248/93‑1 (MQ=255)

                gTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCtttt                                                                         >  2:2933699/1‑93 (MQ=255)

                 tCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCATCGTTCAGCGCCGCAACATCACCTTTTg                                                                        >  1:3110940/1‑93 (MQ=255)

                 tCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTg                                                                        >  1:3977280/1‑93 (MQ=255)

                  cACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGc                                                                       <  1:4710514/93‑1 (MQ=255)

                    cTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtg                                                                     >  2:7409964/1‑93 (MQ=255)

                    cTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtg                                                                     <  2:6655693/93‑1 (MQ=255)

                     tttAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgt                                                                    <  1:8010048/93‑1 (MQ=255)

                     tttAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgt                                                                    <  2:8288598/93‑1 (MQ=255)

                     tttAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgt                                                                    <  1:1577173/93‑1 (MQ=255)

                     tttAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgt                                                                    >  1:7428609/1‑93 (MQ=255)

                     tttAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgt                                                                    <  1:3534689/93‑1 (MQ=255)

                     tttAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgt                                                                    >  1:8495441/1‑93 (MQ=255)

                       tAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgtgt                                                                  <  2:1968484/93‑1 (MQ=255)

                       tAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgtgt                                                                  <  1:2248554/93‑1 (MQ=255)

                        aGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCtgtgtg                                                                 >  1:5398112/1‑93 (MQ=255)

                         gCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGc                                                                <  2:3184137/93‑1 (MQ=255)

                         gCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGc                                                                >  1:9469262/1‑93 (MQ=255)

                         gCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGc                                                                <  1:3681948/93‑1 (MQ=255)

                         gCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGc                                                                <  2:208599/93‑1 (MQ=255)

                         gCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGc                                                                >  2:9106074/1‑93 (MQ=255)

                          cGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCt                                                               <  1:4106637/93‑1 (MQ=255)

                          cGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCt                                                               <  2:7145539/93‑1 (MQ=255)

                           gATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTc                                                              >  2:8895035/1‑93 (MQ=255)

                            aTCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCg                                                             >  2:9608254/1‑93 (MQ=255)

                            aTCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCg                                                             >  1:920867/1‑93 (MQ=255)

                             tCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGt                                                            >  2:6165311/1‑93 (MQ=255)

                               gTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTaa                                                          <  1:9988291/93‑1 (MQ=255)

                               gTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTaa                                                          <  1:2592667/93‑1 (MQ=255)

                                ttttATGCACATCGCGACGCGGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTaaa                                                         >  2:9685953/1‑93 (MQ=255)

                                ttttATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTaaa                                                         <  1:7219835/93‑1 (MQ=255)

                                  ttATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAAcc                                                       <  2:5436207/93‑1 (MQ=255)

                                  ttATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAAcc                                                       <  1:2829195/93‑1 (MQ=255)

                                   tATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCa                                                      <  1:8834724/93‑1 (MQ=255)

                                    aTGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAg                                                     <  1:919113/93‑1 (MQ=255)

                                    aTGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAg                                                     >  2:916788/1‑93 (MQ=255)

                                     tGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGt                                                    >  1:3784589/1‑93 (MQ=255)

                                     tGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGt                                                    >  2:2504147/1‑93 (MQ=255)

                                      gCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGtt                                                   <  2:8211488/93‑1 (MQ=255)

                                      gCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGtt                                                   >  2:4426366/1‑93 (MQ=255)

                                      gCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGCAAACCAGtt                                                   >  2:8361909/1‑93 (MQ=255)

                                          aTCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCa                                               >  2:4773605/1‑93 (MQ=255)

                                           tCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAg                                              >  1:6914682/1‑93 (MQ=255)

                                               gACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGacac                                          <  2:5214990/93‑1 (MQ=255)

                                               gACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGacac                                          >  2:5051658/1‑93 (MQ=255)

                                               gACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGacac                                          >  1:124031/1‑93 (MQ=255)

                                               gACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGacac                                          >  1:3270835/1‑93 (MQ=255)

                                                aCGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACg                                         >  2:3729697/1‑93 (MQ=255)

                                                aCGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACg                                         >  1:5244217/1‑93 (MQ=255)

                                                aCGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACg                                         <  1:8219980/93‑1 (MQ=255)

                                                   cAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGttt                                      >  1:6435713/1‑93 (MQ=255)

                                                    aGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTc                                     <  2:3309164/93‑1 (MQ=255)

                                                    aGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTc                                     <  2:8702956/93‑1 (MQ=255)

                                                    aGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTc                                     >  2:7491220/1‑93 (MQ=255)

                                                     gcgcTTTCTGATTTTCAGTCAGCTCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCa                                    <  2:6936539/93‑1 (MQ=255)

                                                     gcgcTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCa                                    <  1:2388360/93‑1 (MQ=255)

                                                       gcTTTCTGTTTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACACGTCTCCAGACACGTTTCAgg                                  >  1:2038094/1‑93 (MQ=255)

                                                        cTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGa                                 >  2:4350090/1‑93 (MQ=255)

                                                        cTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTAAGGa                                 >  1:5789215/1‑93 (MQ=255)

                                                          ttCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAg                               <  1:10035299/93‑1 (MQ=255)

                                                          ttCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAg                               >  2:3875222/1‑93 (MQ=255)

                                                          ttCTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAg                               <  2:6794495/93‑1 (MQ=255)

                                                            cTGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCg                             >  1:5901919/1‑93 (MQ=255)

                                                             tGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCgg                            <  1:9659261/93‑1 (MQ=255)

                                                             tGATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCgg                            >  2:3927720/1‑93 (MQ=255)

                                                              gATTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACTTTTCGCTGTGTGCTCGTAAACACGTTTACAGACACGGTTAAGGAAgggct                           >  1:291546/1‑89 (MQ=255)

                                                               aTTTTCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGtt                          <  2:2645137/93‑1 (MQ=255)

                                                                ttttCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTa                         <  2:9831342/93‑1 (MQ=255)

                                                                 tttCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAg                        >  2:5840230/1‑93 (MQ=255)

                                                                  ttCAGTCAGCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGc                       >  2:1672811/1‑93 (MQ=255)

                                                                   tcagtcagCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGcc                      >  1:2692020/1‑93 (MQ=255)

                                                                     agtcagCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGcccc                    >  2:269519/1‑93 (MQ=255)

                                                                     agtcagCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGcccc                    <  1:4035241/93‑1 (MQ=255)

                                                                       tcagCGCATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCtt                  >  1:7558068/1‑93 (MQ=255)

                                                                          gcgcATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCca               <  2:4880079/93‑1 (MQ=255)

                                                                             cATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCacac            >  2:9275355/1‑93 (MQ=255)

                                                                             cATCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCacac            >  2:5800638/1‑93 (MQ=255)

                                                                              aTCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACAcc           >  2:7657474/1‑93 (MQ=255)

                                                                               tCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCg          >  1:7742400/1‑93 (MQ=255)

                                                                               tCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCg          >  1:5109134/1‑93 (MQ=255)

                                                                               tCAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCg          >  2:3585743/1‑93 (MQ=255)

                                                                                cAACGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCgg         <  2:5108572/93‑1 (MQ=255)

                                                                                 aaCGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGa        >  2:1247796/1‑93 (MQ=255)

                                                                                 aaCGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGa        <  1:1793906/93‑1 (MQ=255)

                                                                                 aaCGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGa        <  2:1431293/93‑1 (MQ=255)

                                                                                 aaCGTTCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGa        >  2:1477100/1‑93 (MQ=255)

                                                                                     ttCAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTcc    >  1:9786059/1‑93 (MQ=255)

                                                                                       cAGCGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTg  >  2:4013661/1‑93 (MQ=255)

ACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTG  >  NC_000913/671853‑672032

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