BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,956,046	A→C	R287R (CGT→CGG)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,956,046	0	A	C	84.9%	758.7 / 113.3	258	R287R (CGT→CGG)	ptr	protease III
Reads supporting (aligned to +/- strand): ref base A (21/18); new base C (117/102); total (138/120)
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.71e-01

Only 100 of 262 total aligned reads displayed.

CGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTT  >  NC_000913/2955954‑2956138

        aaTCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAg                                                                                      >  1:3643789/1‑93 (MQ=255)

        aaTCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAg                                                                                      >  2:266004/1‑93 (MQ=255)

        aaTCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAg                                                                                      >  2:5190294/1‑93 (MQ=255)

        aaTCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAg                                                                                      >  2:7494212/1‑93 (MQ=255)

        aaTCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAg                                                                                      >  1:19599/1‑93 (MQ=255)

        aaTCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAg                                                                                      <  2:2468982/93‑1 (MQ=255)

          tCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAgcg                                                                                    >  2:781659/1‑93 (MQ=255)

          tCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAgcg                                                                                    >  1:2070294/1‑93 (MQ=255)

            aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCa                                                                                  <  2:5621982/93‑1 (MQ=255)

             gATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAg                                                                                 >  1:6606463/1‑93 (MQ=255)

             gATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAg                                                                                 >  1:7410436/1‑93 (MQ=255)

              aTAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAgg                                                                                <  2:7010516/93‑1 (MQ=255)

               tAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAggg                                                                               >  2:6624132/1‑93 (MQ=255)

                   tAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGt                                                                           >  2:3208095/1‑93 (MQ=255)

                   tAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGt                                                                           <  2:1670686/93‑1 (MQ=255)

                    aaTCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTa                                                                          >  2:3040524/1‑93 (MQ=255)

                    aaTCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTa                                                                          >  2:1291390/1‑93 (MQ=255)

                     aTCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaa                                                                         >  2:288412/1‑93 (MQ=255)

                     aTCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaa                                                                         <  1:1122468/93‑1 (MQ=255)

                      tCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaat                                                                        <  1:7311591/93‑1 (MQ=255)

                       cAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaatg                                                                       >  2:1142633/1‑93 (MQ=255)

                        aaTTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaatga                                                                      <  1:2645188/93‑1 (MQ=255)

                         aTTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaatgaa                                                                     >  2:5663756/1‑93 (MQ=255)

                         aTTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaatgaa                                                                     <  2:5817501/93‑1 (MQ=255)

                          ttCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaatgaat                                                                    <  2:6976847/93‑1 (MQ=255)

                          ttCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaatgaat                                                                    >  2:1129619/1‑93 (MQ=255)

                           tCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTaatgaatg                                                                   >  2:1026480/1‑93 (MQ=255)

                            cATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAatgaatga                                                                  <  2:3790079/93‑1 (MQ=255)

                             aTCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGAt                                                                 >  1:4717623/1‑93 (MQ=255)

                              tCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATa                                                                <  1:7105538/93‑1 (MQ=255)

                              tCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATa                                                                >  1:2962056/1‑93 (MQ=255)

                               cGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATaa                                                               <  1:5770066/93‑1 (MQ=255)

                               cGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATaa                                                               <  2:6084429/93‑1 (MQ=255)

                                   tttACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGcc                                                           >  2:4089035/1‑93 (MQ=255)

                                   tttACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGcc                                                           >  2:1873967/1‑93 (MQ=255)

                                    ttACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGccc                                                          >  2:4531844/1‑93 (MQ=255)

                                    ttACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGccc                                                          >  1:1527733/1‑93 (MQ=255)

                                      actacGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCtt                                                        <  2:3396820/93‑1 (MQ=255)

                                      actacGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCtt                                                        <  2:7447459/93‑1 (MQ=255)

                                       ctacGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCttt                                                       <  1:5151683/93‑1 (MQ=255)

                                       ctacGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCttt                                                       <  2:1338832/93‑1 (MQ=255)

                                        tacGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCtttt                                                      >  2:4125969/1‑93 (MQ=255)

                                        tacGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCtttt                                                      <  1:4629122/93‑1 (MQ=255)

                                        tacGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCtttt                                                      >  2:5956188/1‑93 (MQ=255)

                                        tacGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCtttt                                                      >  1:3189254/1‑93 (MQ=255)

                                        tacGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCtttt                                                      <  1:2105927/93‑1 (MQ=255)

                                         acGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTg                                                     <  2:6546819/93‑1 (MQ=255)

                                          cGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTgc                                                    >  2:4733568/1‑93 (MQ=255)

                                          cGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTgc                                                    <  1:6645862/93‑1 (MQ=255)

                                          cGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTgc                                                    >  1:4935618/1‑93 (MQ=255)

                                           ggAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTgcg                                                   >  1:5439506/1‑93 (MQ=255)

                                           ggAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTgcg                                                   <  2:887020/93‑1 (MQ=255)

                                           ggAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTgcg                                                   >  2:6304835/1‑93 (MQ=255)

                                             aaCTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGcgcg                                                 <  1:1039069/93‑1 (MQ=255)

                                                 tCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCgg                                             <  2:3519058/93‑1 (MQ=255)

                                                  cGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGt                                            <  2:7329020/93‑1 (MQ=255)

                                                     tGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGAc                                         >  1:2715411/1‑93 (MQ=255)

                                                     tGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGAc                                         <  1:3862557/93‑1 (MQ=255)

                                                       aGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTa                                       <  1:421728/93‑1 (MQ=255)

                                                        gttgttATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTAc                                      <  2:2772617/93‑1 (MQ=255)

                                                         ttgttATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTAcc                                     >  2:648127/1‑93 (MQ=255)

                                                          tgttATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCg                                    <  1:5999457/93‑1 (MQ=255)

                                                           gttATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCgg                                   >  1:5741441/1‑93 (MQ=255)

                                                             tATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCa                                 >  1:6202562/1‑93 (MQ=255)

                                                             tATCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATACCCTTTTGCGCGTCGGTGACTACCGGCa                                 >  2:143617/1‑93 (MQ=255)

                                                               tCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGCCGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACg                               >  1:4301678/1‑93 (MQ=255)

                                                               tCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACg                               >  1:651064/1‑93 (MQ=255)

                                                               tCGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACg                               <  1:2476016/93‑1 (MQ=255)

                                                               tCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACg                               >  1:1251268/1‑93 (MQ=255)

                                                               gcGATGCGAATCTAAGCG‑GCAACCCTTTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACg                               <  1:2371668/92‑1 (MQ=255)

                                                                ctaTGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACgg                              >  2:2382877/3‑93 (MQ=255)

                                                                cGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACgg                              <  1:5904252/93‑1 (MQ=255)

                                                                cGATGCGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACgg                              <  1:1634124/93‑1 (MQ=255)

                                                                 gATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGt                             >  2:7083183/1‑93 (MQ=255)

                                                                     cGAAACTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGAtt                         <  2:2086977/93‑1 (MQ=255)

                                                                        aaCTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTcc                      <  2:3327121/93‑1 (MQ=255)

                                                                         aCTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCg                     >  2:5387260/1‑93 (MQ=255)

                                                                         aCTCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCg                     <  2:5737575/93‑1 (MQ=255)

                                                                         aCTCAACGCGCAACAC‑TTTCCGCGGAAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCg                     >  2:6426207/1‑93 (MQ=255)

                                                                         aCTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCg                     <  1:5219613/93‑1 (MQ=255)

                                                                          cTCAACGCGCAACAC‑TTTCCGCGGCAGCGCTGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCgg                    <  1:4615598/93‑1 (MQ=255)

                                                                          cTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCgg                    >  2:3001905/1‑93 (MQ=255)

                                                                           tCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGt                   >  1:3955167/1‑93 (MQ=255)

                                                                           tCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGt                   >  2:2042535/1‑93 (MQ=255)

                                                                           tCAACGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGt                   <  1:5052203/93‑1 (MQ=255)

                                                                           tCAACGCACAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGt                   >  1:2224416/1‑93 (MQ=255)

                                                                             aaCGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGttt                 <  1:5528228/93‑1 (MQ=255)

                                                                              aCGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGtttt                <  2:5128893/93‑1 (MQ=255)

                                                                              aCGCGCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGtttt                >  1:6108302/1‑93 (MQ=255)

                                                                               cgcgCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGttttt               >  1:3349919/1‑93 (MQ=255)

                                                                               cgcgCAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGttttt               <  2:6048635/93‑1 (MQ=255)

                                                                                gcgcAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGtttttt              <  1:2870327/93‑1 (MQ=255)

                                                                                   cAACAC‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGct           >  1:290399/1‑93 (MQ=255)

                                                                                      cac‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGctctc        <  1:6328401/93‑1 (MQ=255)

                                                                                      cac‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGctctc        >  2:4953843/1‑93 (MQ=255)

                                                                                       ac‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGctctct       <  1:197485/93‑1 (MQ=255)

                                                                                       ac‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGctctct       >  1:4435221/1‑93 (MQ=255)

                                                                                       ac‑TTTCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGctctct       <  2:5529791/93‑1 (MQ=255)

                                                                                          tttCCGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCttt     >  1:5336569/1‑93 (MQ=255)

                                                                                          tttACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCttt     <  1:3705216/93‑1 (MQ=255)

CGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACAC‑TTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTT  >  NC_000913/2955954‑2956138

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