BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	158,711	C→T	R128H (CGC→CAC)	pcnB ←	poly(A) polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	158,711	0	C	T	98.8%	570.1 / ‑4.0	166	R128H (CGC→CAC)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand): ref base C (0/2); new base T (88/75); total (89/77)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

Only 100 of 177 total aligned reads displayed.

TACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGA  >  NC_000913/158621‑158800

 aCAGGCTGTTGATAGTGAAATGGCGGCGCTGGGCGGCTTATTCGATGGAGCCGAAAATGTTGTCGCGCAGCAGGATGCCGTTTTGCCCGTGtt                                                                                        >  2:668506/1‑93 (MQ=255)

   aGGCTGTTGATAGTGAAATTGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTgg                                                                                      >  2:4561627/1‑93 (MQ=255)

    ggCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTggg                                                                                     >  1:976040/1‑93 (MQ=255)

     gCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGa                                                                                    <  2:340784/93‑1 (MQ=255)

     gCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGa                                                                                    <  1:3049660/93‑1 (MQ=255)

      cTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAg                                                                                   <  1:1105865/93‑1 (MQ=255)

         ttGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTc                                                                                <  2:2383826/93‑1 (MQ=255)

          tGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCg                                                                               >  1:637895/1‑93 (MQ=255)

          tGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCg                                                                               >  1:1841822/1‑93 (MQ=255)

               gTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCgg                                                                          >  1:2950556/1‑93 (MQ=255)

               gTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCgg                                                                          >  2:3992644/1‑93 (MQ=255)

                tGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGt                                                                         >  1:69178/1‑93 (MQ=255)

                 gAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTc                                                                        <  1:2917596/93‑1 (MQ=255)

                 gAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTc                                                                        >  1:1784565/1‑93 (MQ=255)

                  aaaTCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCg                                                                       >  1:563523/1‑93 (MQ=255)

                   aaTCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTAGc                                                                      >  2:1685375/1‑93 (MQ=255)

                   aaTCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGc                                                                      >  1:3509390/1‑93 (MQ=255)

                     tCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTg                                                                    >  2:4173110/1‑93 (MQ=255)

                     tCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTg                                                                    >  1:4213756/1‑93 (MQ=255)

                        cggcgCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACg                                                                 >  2:3479860/1‑93 (MQ=255)

                         ggcgCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGt                                                                <  1:2571364/93‑1 (MQ=255)

                            gcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTAc                                                             <  1:540734/93‑1 (MQ=255)

                            gcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTAc                                                             <  1:4093549/93‑1 (MQ=255)

                                  cGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgt                                                       <  1:2140007/93‑1 (MQ=255)

                                   gTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtg                                                      >  2:3328938/1‑93 (MQ=255)

                                   gTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtg                                                      <  1:565252/93‑1 (MQ=255)

                                   gTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtg                                                      >  1:4076304/1‑93 (MQ=255)

                                   gTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtg                                                      >  2:3631530/1‑93 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtgg                                                     >  2:4061271/1‑93 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtgg                                                     >  1:1321576/1‑93 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAAAATGCCGTTTTGCCCGTGTTGGGGGGGCGGGCGGTCGCTGACGTTACTTTCttgg                                                     >  2:496968/1‑93 (MQ=255)

                                     cttcttCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtggt                                                    <  2:2954546/93‑1 (MQ=255)

                                     cttctgCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGCGCGGTCGCTGACGTTACCTTCgtggt                                                    <  1:2284094/93‑1 (MQ=255)

                                      ttcttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtggtg                                                   <  1:3054337/93‑1 (MQ=255)

                                        cttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTc                                                 >  1:4182120/1‑93 (MQ=255)

                                        cttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTAGTGGTGTc                                                 >  2:4675220/1‑93 (MQ=255)

                                         ttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTcc                                                <  1:260567/93‑1 (MQ=255)

                                            gATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACg                                             >  1:1649784/1‑93 (MQ=255)

                                            gATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACg                                             <  2:3002296/93‑1 (MQ=255)

                                             aTGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACgg                                            >  2:3112565/1‑93 (MQ=255)

                                              tGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGa                                           >  1:4650705/1‑93 (MQ=255)

                                              tGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGCCGCTGACGTTACCTTCGTGGTGTCCACGGa                                           >  1:3865240/1‑93 (MQ=255)

                                               ggAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGaa                                          <  1:3695073/93‑1 (MQ=255)

                                               ggAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGaa                                          <  1:831990/93‑1 (MQ=255)

                                               ggAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGaa                                          <  2:1885890/93‑1 (MQ=255)

                                                gAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAg                                         >  1:2127784/1‑93 (MQ=255)

                                                 aGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAgg                                        <  2:570065/93‑1 (MQ=255)

                                                 aGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAgg                                        <  1:278722/93‑1 (MQ=255)

                                                  gCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGt                                       <  1:217814/93‑1 (MQ=255)

                                                  gCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGt                                       <  1:4425350/93‑1 (MQ=255)

                                                   ccGAAAATGTTGTCGCGCAGCAACGAGCCGTTTTGCCAGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTc                                      <  1:3123404/93‑1 (MQ=255)

                                                   ccGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTc                                      >  1:1779908/1‑93 (MQ=255)

                                                   ccGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTc                                      <  1:1923621/93‑1 (MQ=255)

                                                    cGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCg                                     >  2:2256036/1‑93 (MQ=255)

                                                    cGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCg                                     >  1:2614484/1‑93 (MQ=255)

                                                    cGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCg                                     >  1:4377026/1‑93 (MQ=255)

                                                     gAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGc                                    >  1:1203604/1‑93 (MQ=255)

                                                      aaaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCa                                   >  2:2106220/1‑93 (MQ=255)

                                                      aaaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCa                                   >  2:2392065/1‑93 (MQ=255)

                                                      aaaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCa                                   >  1:2111138/1‑93 (MQ=255)

                                                       aaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCaa                                  >  2:4532708/1‑93 (MQ=255)

                                                       aaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCaa                                  <  1:3882686/93‑1 (MQ=255)

                                                        aaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAAc                                 <  2:13619/93‑1 (MQ=255)

                                                        aaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAAc                                 <  1:3569638/93‑1 (MQ=255)

                                                        aaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAAc                                 >  1:3816140/1‑93 (MQ=255)

                                                         aTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACt                                >  2:3674431/1‑93 (MQ=255)

                                                         aTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACt                                >  1:2526432/1‑93 (MQ=255)

                                                          tgttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACtt                               <  2:294160/93‑1 (MQ=255)

                                                          tgttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACtt                               >  2:3705323/1‑93 (MQ=255)

                                                          tgttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACtt                               >  2:4368056/1‑93 (MQ=255)

                                                           gttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTc                              <  2:288171/93‑1 (MQ=255)

                                                             tgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGa                            <  1:352925/93‑1 (MQ=255)

                                                              gtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGAt                           <  1:3562485/93‑1 (MQ=255)

                                                                 gcgcAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAt                        <  2:627495/93‑1 (MQ=255)

                                                                  cgcAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAtc                       <  1:3427390/93‑1 (MQ=255)

                                                                   gcagcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAtct                      >  2:1160537/1‑93 (MQ=255)

                                                                    cagcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAtctc                     >  1:1503896/1‑93 (MQ=255)

                                                                     agcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTcc                    <  2:1658178/93‑1 (MQ=255)

                                                                      gcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCg                   <  2:1243153/93‑1 (MQ=255)

                                                                      gcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCg                   <  1:3792120/93‑1 (MQ=255)

                                                                      gcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCg                   >  2:4409212/1‑93 (MQ=255)

                                                                      gcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGCCGCAACTTCGATAATCTCCg                   >  2:1609208/1‑93 (MQ=255)

                                                                       caACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCgg                  <  1:2211424/93‑1 (MQ=255)

                                                                        aaCATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCggg                 >  1:4285955/1‑93 (MQ=255)

                                                                        aaCATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCggg                 >  1:1221445/1‑93 (MQ=255)

                                                                        aaCATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCggg                 <  2:4262585/93‑1 (MQ=255)

                                                                         aCATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGc                >  2:2531121/1‑93 (MQ=255)

                                                                             ggcGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCaaa            >  1:4530588/3‑93 (MQ=255)

                                                                              ccGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAAc           >  1:1151407/1‑93 (MQ=255)

                                                                              ccGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAAc           >  2:3432078/1‑93 (MQ=255)

                                                                               cGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACa          <  2:4155648/93‑1 (MQ=255)

                                                                                 ttttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAtt        <  1:4290397/93‑1 (MQ=255)

                                                                                 ttttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAtt        >  2:954353/1‑93 (MQ=255)

                                                                                  tttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTa       <  2:1806226/93‑1 (MQ=255)

                                                                                  tttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTa       <  1:3724875/93‑1 (MQ=255)

                                                                                    tGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACa     <  2:883424/93‑1 (MQ=255)

                                                                                     gCCCGTGTTGGGAGGGCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGGCGCAACTTCGATAATCTCCGGGCCAAACATTACAt    >  1:1978267/1‑93 (MQ=255)

                                                                                      cccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATg   <  1:78548/93‑1 (MQ=255)

                                                                                      cccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCATCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATg   <  1:1722429/93‑1 (MQ=255)

                                                                                       ccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGa  <  1:3635053/93‑1 (MQ=255)

TACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGA  >  NC_000913/158621‑158800

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