BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,108,443	C→A	noncoding (56/76 nt)	pheV →	tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,108,443	0	C	A	99.6%	855.4 / ‑5.0	230	noncoding (56/76 nt)	pheV	tRNA‑Phe
Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (105/124); total (106/124)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 234 total aligned reads displayed.

CGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATT  >  NC_000913/3108353‑3108536

tgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGAt                                                                                            <  1:3869230/93‑1 (MQ=255)

 gagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGAtt                                                                                           >  1:1730536/1‑94 (MQ=255)

 gagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGAta                                                                                           >  1:1284016/1‑93 (MQ=255)

  ggacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTc                                                                                          <  1:6337898/93‑1 (MQ=255)

   gacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTcg                                                                                         >  1:364343/1‑93 (MQ=255)

   gacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTcg                                                                                         >  1:2220638/1‑93 (MQ=255)

    acgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTcca                                                                                        >  1:293978/1‑93 (MQ=255)

      gagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAg                                                                                      >  1:5829828/1‑94 (MQ=255)

       tggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGt                                                                                     <  1:6991090/93‑1 (MQ=255)

       agGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAga                                                                                     >  1:3165527/1‑93 (MQ=255)

       agGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAga                                                                                     >  1:7372745/1‑93 (MQ=255)

        ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGta                                                                                    >  1:3298592/1‑93 (MQ=255)

          ggAATCAGGTTTAATGCGCCCTGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCg                                                                                  <  1:4964364/93‑1 (MQ=255)

          ggAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCg                                                                                  <  1:7401577/93‑1 (MQ=255)

           gAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgg                                                                                 <  1:5584428/94‑1 (MQ=255)

           caaTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgg                                                                                 <  1:4322483/93‑1 (MQ=25)

             aTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGAtt                                                                                           >  2:6304496/1‑82 (MQ=25)

             aTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgggg                                                                               >  1:5851301/1‑93 (MQ=21)

             aTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGc                                                                               >  1:1466752/1‑94 (MQ=21)

              tCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTc                                                                                          >  2:6649332/1‑82 (MQ=18)

              tCAGGTTTAATGCGCCCCGTTGCCCGGATAGATCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTc                                                                                          <  2:6431936/82‑1 (MQ=255)

                    ttAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTa                                                                        >  1:1640512/1‑94 (MQ=18)

                      aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaaa                                                                      >  1:4718482/1‑93 (MQ=21)

                        ggcgcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTc                                                                    <  1:1830236/93‑1 (MQ=255)

                          ggcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCtt                                                                  <  1:4634835/93‑1 (MQ=255)

                            ccccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTat                                                                >  1:5172763/1‑93 (MQ=255)

                              ccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcaccac                                                                          >  2:3942583/1‑82 (MQ=14)

                              ccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAgc                                                              >  1:7081106/1‑93 (MQ=255)

                               cGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACt                                                                         >  2:6130686/1‑82 (MQ=17)

                               cGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACt                                                                         <  2:6711676/82‑1 (MQ=17)

                               cGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGag                                                             >  1:1172571/1‑93 (MQ=255)

                                 ttGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                                       >  2:4761793/1‑82 (MQ=25)

                                  tGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAAt                                                                      >  2:697381/1‑82 (MQ=37)

                                   gCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAAtt                                                                     >  2:2371364/1‑82 (MQ=255)

                                    cccGGATAGCTCAGTCGGTAGAGCAGGGCATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTc                                                                    <  2:1435166/82‑1 (MQ=14)

                                     tcGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGcc                                                       <  1:2453743/93‑1 (MQ=255)

                                     gcGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGcc                                                       <  1:6840324/93‑1 (MQ=255)

                                      cGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCtt                                                                  >  2:5796187/1‑82 (MQ=255)

                                      cGGATAGCGCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCtt                                                                  <  2:5464757/82‑1 (MQ=255)

                                       ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTa                                                                 <  2:2186504/82‑1 (MQ=255)

                                       ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccct                                                     >  1:3457550/1‑93 (MQ=255)

                                        gATAGCTCAGTCGGTAGAGCAGGGGATTGAAGATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcac                                                    <  1:2256016/94‑1 (MQ=255)

                                        gATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                                <  2:1237365/82‑1 (MQ=255)

                                        gATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                                >  2:5576600/1‑82 (MQ=255)

                                         gtAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcaca                                                   <  1:3802979/93‑1 (MQ=255)

                                          taacTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGa                                                              >  2:7335678/4‑82 (MQ=255)

                                          tAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGa                                                              <  2:226058/82‑1 (MQ=255)

                                          tAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                                  <  1:4942531/94‑1 (MQ=255)

                                          gaGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                                  <  1:6469569/93‑1 (MQ=255)

                                           tgCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAg                                                 <  1:6046994/93‑1 (MQ=255)

                                            gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATGCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAgg                                                <  1:6069004/94‑1 (MQ=255)

                                               gaGTCGGTAGAGCTGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCgg                                             <  1:2852194/93‑1 (MQ=255)

                                                 gTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                                           >  1:625051/1‑94 (MQ=255)

                                                  tCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccc                                                      >  2:1957423/1‑82 (MQ=255)

                                                  tCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccc                                                      <  2:6714950/82‑1 (MQ=255)

                                                   tggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttt                                         <  1:5080330/93‑1 (MQ=255)

                                                    tgTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtttt                                        <  1:6780364/93‑1 (MQ=255)

                                                    ggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcac                                                    >  2:54319/1‑82 (MQ=255)

                                                       ggagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGc                                     <  1:4574009/93‑1 (MQ=255)

                                                       agagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAg                                                 >  2:2817808/1‑82 (MQ=255)

                                                         agCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGc                                               >  2:3523343/1‑82 (MQ=255)

                                                           gaGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtttt                                 <  1:298151/93‑1 (MQ=255)

                                                            gggggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg                                <  1:6902701/93‑1 (MQ=255)

                                                            aGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCttttt                                >  1:5660149/1‑93 (MQ=255)

                                                              gggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtt                                          <  2:26911/82‑1 (MQ=255)

                                                                taTTGAAAATCCCCTTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGACCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                            <  1:108600/93‑1 (MQ=37)

                                                                gATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtttt                                        >  2:3794570/1‑82 (MQ=255)

                                                                    gAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCt                                    >  2:4428224/1‑82 (MQ=255)

                                                                     aaaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt                                   <  2:1782849/82‑1 (MQ=255)

                                                                      aaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCttt                                  <  2:7311902/82‑1 (MQ=255)

                                                                      aaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACat                      >  1:2085137/1‑93 (MQ=255)

                                                                       aaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtttt                                 <  2:2461576/82‑1 (MQ=255)

                                                                       aaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtttt                                 <  2:6503791/82‑1 (MQ=255)

                                                                       aaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACaaa                     >  1:5746929/1‑93 (MQ=255)

                                                                        aTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg                                >  2:3830587/1‑82 (MQ=255)

                                                                         tCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTgg                               <  2:3223558/82‑1 (MQ=255)

                                                                         gccccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATaa                   <  1:1998270/93‑1 (MQ=255)

                                                                          ccccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAc                  <  1:7261174/94‑1 (MQ=255)

                                                                           gccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAcc                 <  1:2527405/93‑1 (MQ=255)

                                                                           cccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGAt                             >  2:355947/1‑82 (MQ=255)

                                                                            gcGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCt                <  1:1907026/93‑1 (MQ=255)

                                                                             cGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                           <  2:5687736/82‑1 (MQ=255)

                                                                              ttgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTCAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTc              <  1:888043/93‑1 (MQ=255)

                                                                              gtgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCttg              >  1:590784/1‑93 (MQ=255)

                                                                              gtgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTc              <  1:2694279/94‑1 (MQ=255)

                                                                               tgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGa                         >  2:6210822/1‑82 (MQ=255)

                                                                                gtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAc                        <  2:3568741/82‑1 (MQ=255)

                                                                                gtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCaa            >  1:4120392/1‑93 (MQ=255)

                                                                                gtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCAc            <  1:3987039/94‑1 (MQ=255)

                                                                                  ccTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGACCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACaa                      <  2:3082181/82‑1 (MQ=255)

                                                                                   cTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAAt                     <  2:5024371/82‑1 (MQ=255)

                                                                                   cTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAAt                     <  2:4584560/82‑1 (MQ=255)

                                                                                    ttGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATa                    <  2:4206723/82‑1 (MQ=255)

                                                                                    ttGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATa                    >  2:3714695/1‑82 (MQ=255)

                                                                                     tGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATaa                   <  2:5837158/82‑1 (MQ=255)

                                                                                     gggTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaa       <  1:3056088/93‑1 (MQ=255)

                                                                                      ggTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaa      <  1:3215546/94‑1 (MQ=255)

                                                                                      agTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaa      <  1:5060005/93‑1 (MQ=255)

                                                                                       gTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaag     >  1:367691/1‑93 (MQ=255)

                                                                                          cGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAAtt  <  1:456940/94‑1 (MQ=255)

CGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATT  >  NC_000913/3108353‑3108536

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