BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,927,241	(C)_6→7	coding (1599/2061 nt)	ptrB ←	oligopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,927,235	1	.	C	99.6%	982.5 / ‑3.2	233	coding (1605/2061 nt)	ptrB	oligopeptidase B
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (107/125); total (108/125)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 277 total aligned reads displayed.

ATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑CCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGT  >  NC_000913/1927135‑1927330

 ttCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑c                                                                                                >  2:3525917/1‑101 (MQ=255)

 ttCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑c                                                                                                >  1:3525916/1‑101 (MQ=255)

  tcatcaAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cc                                                                                               <  1:3525099/101‑1 (MQ=255)

   catcaAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑ccc                                                                                              <  1:3525342/101‑1 (MQ=255)

   catcaAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑ccc                                                                                              >  2:3525921/1‑101 (MQ=255)

    atcaAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccc                                                                                             >  1:3525923/1‑101 (MQ=255)

      caAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccc                                                                                           >  1:3525927/1‑101 (MQ=255)

       aaGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑ccccccc                                                                                          >  2:3525930/1‑100 (MQ=255)

         gCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccccgc                                                                                        >  2:3525933/1‑98 (MQ=255)

          cATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccccgca                                                                                       >  1:3525934/1‑97 (MQ=255)

          cATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccccgca                                                                                       <  2:3525468/101‑5 (MQ=255)

          cATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccccgca                                                                                       <  1:3525400/101‑5 (MQ=255)

           aTCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccccgcac                                                                                      >  1:3525940/1‑96 (MQ=255)

           aTCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccccgcac                                                                                      >  2:3525939/1‑96 (MQ=255)

               ttgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTccc                                                                                   <  2:3525512/101‑1 (MQ=255)

                tgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGTAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCg                                                                                  >  1:3525943/1‑101 (MQ=255)

                tgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCg                                                                                  <  2:3525326/101‑1 (MQ=255)

                tgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCg                                                                                  >  2:3525944/1‑101 (MQ=255)

                 gttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGc                                                                                 <  1:3525483/101‑1 (MQ=255)

                  ttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGcc                                                                                <  2:3525526/101‑1 (MQ=255)

                  ttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGcc                                                                                <  1:3525585/101‑1 (MQ=255)

                   tACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGccc                                                                               <  2:3525508/101‑1 (MQ=255)

                    acaacaTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCa                                                                              >  2:3525949/1‑101 (MQ=255)

                    acaacaTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCa                                                                              >  1:3525950/1‑101 (MQ=255)

                     caacaTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCAt                                                                             <  1:3525460/101‑1 (MQ=255)

                      aacaTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATc                                                                            >  2:3525953/1‑101 (MQ=255)

                      aacaTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATc                                                                            >  1:3525957/1‑101 (MQ=255)

                       acaTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCg                                                                           >  2:3525958/1‑101 (MQ=255)

                         aTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCa                                                                         >  2:3525962/1‑101 (MQ=255)

                            aacaaACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATaa                                                                      >  1:3525965/1‑101 (MQ=255)

                                aaCGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACaaa                                                                  <  2:3525809/101‑1 (MQ=255)

                                 aCGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAg                                                                 >  1:3525970/1‑101 (MQ=255)

                                 aCGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAg                                                                 >  2:3525969/1‑101 (MQ=255)

                                  cGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGc                                                                <  1:3525636/101‑1 (MQ=255)

                                  cGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGc                                                                <  2:3525256/101‑1 (MQ=255)

                                   ggTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCg                                                               <  2:3525725/101‑1 (MQ=255)

                                   ggTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCg                                                               <  2:3525584/101‑1 (MQ=255)

                                   ggTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCg                                                               >  2:3525973/1‑101 (MQ=255)

                                    gTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGc                                                                >  1:3525978/1‑99 (MQ=255)

                                        cTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCg                                                               <  2:3525633/96‑1 (MQ=255)

                                        cTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGga                                                          <  1:3525670/101‑1 (MQ=255)

                                        cTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAACTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGga                                                          >  2:3525980/1‑101 (MQ=255)

                                          gggCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGgaga                                                        >  1:3525983/1‑101 (MQ=255)

                                            gCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGc                                                      <  2:3525600/101‑1 (MQ=255)

                                             cGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGcc                                                     <  2:3525451/101‑1 (MQ=255)

                                              gATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCa                                                    <  2:3525424/101‑1 (MQ=255)

                                              gATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCa                                                    >  2:3525988/1‑101 (MQ=255)

                                              gATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCa                                                    >  1:3525987/1‑101 (MQ=255)

                                              gATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccc                                                                                           >  1:3525990/1‑61 (MQ=255)

                                               aTAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCAt                                                   <  1:3525156/101‑1 (MQ=255)

                                                tAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGc                                                                >  1:3525993/1‑87 (MQ=255)

                                                 aaCGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAg                                                 <  1:3525417/101‑1 (MQ=255)

                                                 aaCGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAg                                                 >  2:3525997/1‑101 (MQ=255)

                                                   cGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGgaga                                                        >  2:3526000/1‑92 (MQ=255)

                                                   cGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGcc                                               <  1:3525629/101‑1 (MQ=255)

                                                     ccGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCa                                             >  1:3526003/1‑101 (MQ=255)

                                                     ccGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCa                                             <  2:3525492/101‑1 (MQ=255)

                                                       gTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGt                                           >  2:3526007/1‑101 (MQ=255)

                                                         ggAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGttt                                         <  2:3525539/101‑1 (MQ=255)

                                                         ggAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGttt                                         >  1:3526009/1‑101 (MQ=255)

                                                         ggAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGttt                                         <  1:3525709/101‑1 (MQ=255)

                                                          gAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGtttt                                        <  2:3525594/101‑1 (MQ=255)

                                                              aatTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACa                                    <  2:3525746/101‑1 (MQ=255)

                                                               atTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACaa                                   >  2:3526013/1‑101 (MQ=255)

                                                               atTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACaa                                   <  1:3525607/101‑1 (MQ=255)

                                                                ttCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAAt                                  >  2:3526016/1‑101 (MQ=255)

                                                                 tCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATg                                 >  1:3526017/1‑101 (MQ=255)

                                                                 tCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATg                                 <  2:3525215/101‑1 (MQ=255)

                                                                  ccGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGc                                <  1:3525502/101‑1 (MQ=255)

                                                                   cGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCa                               <  2:3525477/101‑1 (MQ=255)

                                                                    gggCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCAt                              <  1:3525579/101‑1 (MQ=255)

                                                                     ggCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATc                             <  2:3525661/101‑1 (MQ=255)

                                                                      gCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCg                            <  2:3525571/101‑1 (MQ=255)

                                                                      gCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCg                            >  2:3526021/1‑101 (MQ=255)

                                                                       cGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGc                           >  2:3526023/1‑101 (MQ=255)

                                                                        gTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCa                          <  2:3525674/101‑1 (MQ=255)

                                                                        gTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCa                          >  1:3526024/1‑101 (MQ=255)

                                                                        gTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCa                          <  2:3525718/101‑1 (MQ=255)

                                                                        gTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCa                          <  1:3525522/101‑1 (MQ=255)

                                                                          tGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAgg                        >  2:3526028/1‑101 (MQ=255)

                                                                            aTTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCa                      <  1:3525673/101‑1 (MQ=255)

                                                                             ttAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCAt                     <  1:3525438/101‑1 (MQ=255)

                                                                             ttAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCAt                     >  2:3526032/1‑101 (MQ=255)

                                                                              tAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATc                    >  1:3526033/1‑101 (MQ=255)

                                                                               aaTTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCa                   >  1:3526037/1‑101 (MQ=255)

                                                                                aTTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCaa                  <  1:3525557/101‑1 (MQ=255)

                                                                                aTTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCaa                  <  2:3525550/101‑1 (MQ=255)

                                                                                aTTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCaa                  <  1:3525889/101‑1 (MQ=255)

                                                                                  tGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCa                   >  1:3526040/1‑98 (MQ=255)

                                                                                  tGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGa                <  2:3525370/101‑1 (MQ=255)

                                                                                    cAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATa              <  2:3525641/101‑1 (MQ=255)

                                                                                     gaCGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATaa             >  1:3526042/2‑101 (MQ=255)

                                                                                     aaCGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATaa             <  1:3525648/101‑1 (MQ=255)

                                                                                       cGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATc           >  2:3526045/1‑101 (MQ=255)

                                                                                        gCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCa          <  2:3525517/101‑1 (MQ=255)

                                                                                         cccATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCAt         >  1:3526047/1‑101 (MQ=255)

                                                                                            aTCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTaa      <  1:3525808/101‑1 (MQ=255)

                                                                                             tCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTaaa     <  2:3525615/101‑1 (MQ=255)

                                                                                              cAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAAc    <  2:3525536/101‑1 (MQ=255)

                                                                                                aCATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGt  >  1:3526052/1‑101 (MQ=255)

ATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑CCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGT  >  NC_000913/1927135‑1927330

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