BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,954,264	A→T	D881E (GAT→GAA)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,954,264	0	A	T	99.5%	674.4 / NA	189	D881E (GAT→GAA)	ptr	protease III
Reads supporting (aligned to +/- strand): ref base A (1/0); new base T (83/105); total (84/105)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

CTGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGATCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTGCGC  >  NC_000913/2954174‑2954353

 tGATGTATGAAATCAGCAAGTTTTTGCTGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCg                                                                                        <  2:2045846/93‑1 (MQ=255)

   aTGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGTGTTCGaa                                                                                      <  2:3650053/93‑1 (MQ=255)

    tGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGTGCATATTGCCGCGTTCGaaa                                                                                     <  2:3831977/93‑1 (MQ=255)

    tGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGaaa                                                                                     <  1:4636252/93‑1 (MQ=255)

     ggAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGTTTTTATCGCGCGAATCGAAGCGCATATTGCCGCGTTCGAAAt                                                                                    <  1:2507648/93‑1 (MQ=255)

     ggAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAAt                                                                                    >  2:1653179/1‑93 (MQ=255)

        agaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCtt                                                                                 >  1:1475617/1‑93 (MQ=255)

        agaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCtt                                                                                 <  2:3168351/93‑1 (MQ=255)

        agaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCtt                                                                                 <  1:2252572/93‑1 (MQ=255)

         gaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCttt                                                                                >  1:4466680/1‑93 (MQ=255)

          aaaTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTa                                                                               <  1:464892/93‑1 (MQ=255)

           aaTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTAc                                                                              <  1:2988988/93‑1 (MQ=255)

           aaTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTAc                                                                              <  1:2305124/93‑1 (MQ=255)

            aTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACt                                                                             <  2:2955535/93‑1 (MQ=255)

             tCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACtt                                                                            >  1:3215392/1‑93 (MQ=255)

              cAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTa                                                                           >  1:335152/1‑93 (MQ=255)

                gCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAAc                                                                         >  2:4455373/1‑93 (MQ=255)

                gCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAAc                                                                         <  1:1699391/93‑1 (MQ=255)

                gCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAAc                                                                         <  2:4031913/93‑1 (MQ=255)

                   aGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTc                                                                      >  2:2382855/1‑93 (MQ=255)

                    gTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCg                                                                     <  2:3588048/93‑1 (MQ=255)

                     tttttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGa                                                                    <  1:3814552/93‑1 (MQ=255)

                     tttttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGa                                                                    <  1:3524322/93‑1 (MQ=255)

                      ttttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGAt                                                                   <  2:2021196/93‑1 (MQ=255)

                      ttttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGAt                                                                   <  1:1619651/93‑1 (MQ=255)

                         tGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGct                                                                <  2:534612/93‑1 (MQ=255)

                           cggcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttc                                                              <  1:1647714/93‑1 (MQ=255)

                            ggcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCACATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttct                                                             >  1:2687150/1‑93 (MQ=255)

                            ggccTCAGCAGTTTTATCTGGTCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCTCGTTCGAAATCTTTACTTAACTTCGATGcttct                                                             <  2:1037413/93‑1 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttctt                                                            >  2:387908/1‑93 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttctt                                                            <  1:588313/93‑1 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttctt                                                            <  2:2744834/93‑1 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttctt                                                            <  2:2314123/93‑1 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttctt                                                            <  1:3082197/93‑1 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCCATGcttctt                                                            >  2:2367680/1‑93 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGACACGATTTTATGACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttctt                                                            >  2:3064069/1‑93 (MQ=255)

                                tCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGc                                                         >  1:388272/1‑93 (MQ=255)

                                 cagcagTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGcc                                                        <  1:3912813/93‑1 (MQ=255)

                                 cagcagTTTTATCGGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGcc                                                        <  2:2099676/93‑1 (MQ=255)

                                   gcagTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGa                                                      >  1:4582392/1‑93 (MQ=255)

                                     agTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGc                                                    <  1:1581963/93‑1 (MQ=255)

                                      ggtttATCTGGGCCACGATGTTATCATGCGAAGCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACATCGATGCTGCTTCGCCGAGCg                                                   >  1:4687020/3‑93 (MQ=255)

                                      gTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCg                                                   >  2:3322781/1‑93 (MQ=255)

                                      gTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCg                                                   >  2:3085145/1‑93 (MQ=255)

                                      gTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCCAGCg                                                   >  2:3852422/1‑93 (MQ=255)

                                       ttttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTNCGATGCTTCTTCGCCGAGCGt                                                  >  1:2230162/1‑93 (MQ=255)

                                        tttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGtt                                                 <  1:3809167/93‑1 (MQ=255)

                                         ttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGttt                                                <  2:963664/93‑1 (MQ=255)

                                            tCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCg                                             >  2:4558178/1‑93 (MQ=255)

                                            tCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCg                                             >  2:571830/1‑93 (MQ=255)

                                            tCTGGGCCACGATTTTATCACGCGAATCGAAGAGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGACGAGCGTTTGCg                                             >  2:4347830/1‑93 (MQ=255)

                                             cTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCgg                                            <  2:4152773/93‑1 (MQ=255)

                                             cTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCgg                                            <  2:507090/93‑1 (MQ=255)

                                             cTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCgg                                            <  2:4534270/93‑1 (MQ=255)

                                                 gCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGcc                                        >  2:1577389/1‑93 (MQ=255)

                                                 gCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGAGGTGcc                                        >  1:754369/1‑93 (MQ=255)

                                                  ccACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCt                                       >  1:106408/1‑93 (MQ=255)

                                                   cACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTg                                      >  1:3483773/1‑93 (MQ=255)

                                                   cACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTg                                      <  2:3384404/93‑1 (MQ=255)

                                                    aCGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgc                                     <  1:1203685/93‑1 (MQ=255)

                                                    aCGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgc                                     >  1:4327357/1‑93 (MQ=255)

                                                    aCGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgc                                     >  2:2648740/1‑93 (MQ=255)

                                                     cGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgca                                    <  2:3491001/93‑1 (MQ=255)

                                                     cGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgca                                    <  1:499260/93‑1 (MQ=255)

                                                      gATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcag                                   >  2:63646/1‑93 (MQ=255)

                                                      gATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcag                                   >  1:3425968/1‑93 (MQ=255)

                                                       aTTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcagc                                  >  2:237154/1‑93 (MQ=255)

                                                          ttATCACGCGAATCGAAGCGCATATTGCCGCGTTTGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATc                               <  2:367119/93‑1 (MQ=255)

                                                           tATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCt                              >  1:701174/1‑93 (MQ=255)

                                                             tCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTgg                            <  2:171475/93‑1 (MQ=255)

                                                             tCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTgg                            <  1:280007/93‑1 (MQ=255)

                                                             tCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTgg                            >  2:1589305/1‑93 (MQ=255)

                                                              cACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTggg                           >  1:1979688/1‑93 (MQ=255)

                                                               aCGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGt                          >  1:2042922/1‑93 (MQ=255)

                                                               aCGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGt                          <  1:2045823/93‑1 (MQ=255)

                                                               aCGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGt                          >  2:2662722/1‑93 (MQ=255)

                                                               aCGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGt                          >  1:1969458/1‑93 (MQ=255)

                                                                 gcgAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTaa                        <  1:363782/93‑1 (MQ=255)

                                                                  cgAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAAt                       <  2:2556881/93‑1 (MQ=255)

                                                                  cgAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAAt                       <  2:2481068/93‑1 (MQ=255)

                                                                  cgAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAAt                       <  1:3086592/93‑1 (MQ=255)

                                                                   gAATCGAAGCGAATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAAtt                      >  1:2908962/1‑93 (MQ=255)

                                                                    aaTCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTa                     >  2:2830965/1‑93 (MQ=255)

                                                                    aaTCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTa                     <  2:76596/93‑1 (MQ=255)

                                                                     aTCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTAc                    <  2:3951124/93‑1 (MQ=255)

                                                                      tCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTAcc                   >  1:1109553/1‑93 (MQ=255)

                                                                      tCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTAcc                   >  2:4507201/1‑93 (MQ=255)

                                                                        gAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGc                 >  2:1103690/1‑93 (MQ=255)

                                                                         aaGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGcc                <  2:1489464/93‑1 (MQ=255)

                                                                            cgcAGATTGCCGCGTTCGAAATCTTTACTTAGCTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctgc             <  1:3186221/93‑1 (MQ=255)

                                                                             gcATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGgctgct            >  1:3562555/1‑93 (MQ=255)

                                                                             gcATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctgct            >  2:2178959/1‑93 (MQ=255)

                                                                              cATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctgctg           <  2:1136884/93‑1 (MQ=255)

                                                                                   tGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCACCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGAttt      <  2:4096170/93‑1 (MQ=255)

                                                                                    gCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTg     >  2:1853944/1‑93 (MQ=255)

                                                                                    gCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTg     >  2:1300421/1‑93 (MQ=255)

                                                                                     ccGCGTTCGAAATCTTTACTTAGCTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTgc    <  2:4491269/93‑1 (MQ=255)

                                                                                     ccGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTgc    <  1:973654/93‑1 (MQ=255)

                                                                                       gcgTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTgcgc  >  2:4235337/1‑93 (MQ=255)

                                                                                       gcgATCGAAATCTCTACTTAACTCCGATGCTTCTCCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGTGTAATTACCGACCGCTGGATTcgcga  >  2:3115115/1‑92 (MQ=255)

CTGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGATCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTGCGC  >  NC_000913/2954174‑2954353

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