BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,359	0	C	T	51.3%	123.5 / inf	407	I179I (ATC→ATT)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (103/95); new base T (103/106); total (206/201)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 417 total aligned reads displayed.

ACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT  >  NC_000913/1949267‑1949451

aCCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATc                                                                                              >  2:5986189/1‑93 (MQ=255)

 ccgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCa                                                                                             <  2:6285048/93‑1 (MQ=255)

  cgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCaa                                                                                            <  1:3622138/93‑1 (MQ=255)

  cgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTaa                                                                                            >  2:4493613/1‑93 (MQ=255)

   gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAAc                                                                                           >  2:3819316/1‑93 (MQ=255)

    cccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTaatc                                                                                          <  2:1978601/93‑3 (MQ=255)

     ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCa                                                                                         <  1:4501142/93‑1 (MQ=255)

     ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCa                                                                                         >  1:5997900/1‑93 (MQ=255)

     ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCa                                                                                         >  1:5598301/1‑93 (MQ=255)

      cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCat                                                                                        >  1:5559000/1‑93 (MQ=255)

      cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCat                                                                                        >  2:1752979/1‑93 (MQ=255)

       gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCata                                                                                       <  1:3285504/93‑1 (MQ=255)

       gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCata                                                                                       >  2:4414491/1‑93 (MQ=255)

        cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCATCGCCGAGCAACACAATAACAGCATTAATCatat                                                                                      <  1:6142919/93‑1 (MQ=255)

        cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCatat                                                                                      <  2:1092166/93‑1 (MQ=255)

        cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCatat                                                                                      >  2:1534855/1‑93 (MQ=255)

        cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCatat                                                                                      <  1:2685733/93‑1 (MQ=255)

        cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCatat                                                                                      <  2:5698227/93‑1 (MQ=255)

          aTGCCTGGGAACTTGGCTTTAATCTGGTGAGTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACCGTATCAACCGTATCt                                                                                    >  1:2289789/1‑93 (MQ=255)

           tGCCTGGTAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                                   <  2:5119765/93‑1 (MQ=255)

           tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACTCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                                   <  1:4127379/93‑1 (MQ=255)

           tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                                   <  1:5033409/93‑1 (MQ=255)

           tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTa                                                                                   >  1:5465278/1‑93 (MQ=255)

            gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAc                                                                                  <  1:2195741/93‑1 (MQ=255)

            gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAc                                                                                  >  1:2527950/1‑93 (MQ=255)

             ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAcc                                                                                 <  2:5793710/93‑1 (MQ=255)

             ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAcc                                                                                 <  1:2002623/93‑1 (MQ=255)

              cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                                >  1:3544928/1‑93 (MQ=255)

                gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgc                                                                              >  2:6169259/1‑93 (MQ=255)

                gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcgc                                                                              >  2:5929166/1‑93 (MQ=255)

                gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcgc                                                                              >  2:2912689/1‑93 (MQ=255)

                  gAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACACgcgc                                                                            >  1:2878886/1‑93 (MQ=255)

                   aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCa                                                                           <  2:8461/93‑1 (MQ=255)

                   aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCa                                                                           >  1:4216178/1‑93 (MQ=255)

                    aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCAt                                                                          >  2:5185421/1‑93 (MQ=255)

                    aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCAt                                                                          >  1:3466350/1‑93 (MQ=255)

                      ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCg                                                                        >  1:5326552/1‑93 (MQ=255)

                       tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGc                                                                       >  2:3102931/1‑93 (MQ=255)

                        ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                      >  2:3349003/1‑93 (MQ=255)

                        ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                      <  1:1536782/93‑1 (MQ=255)

                        ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                      <  1:5229274/93‑1 (MQ=255)

                         gCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGccc                                                                     <  1:4347084/93‑1 (MQ=255)

                          cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCg                                                                    >  1:2051711/1‑93 (MQ=255)

                            ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtg                                                                  <  2:5269891/93‑1 (MQ=255)

                            ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtg                                                                  >  1:4808502/1‑93 (MQ=255)

                            ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgtg                                                                  <  1:3654317/93‑1 (MQ=255)

                            ttAATCTGGTGATTGCCGAAGACGCCTGTAGTCCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtg                                                                  <  2:2547682/93‑1 (MQ=255)

                             tAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtgt                                                                 >  2:6080507/1‑93 (MQ=255)

                              aaTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGtgtg                                                                >  2:3483120/1‑93 (MQ=255)

                                tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCg                                                              >  1:1421602/1‑93 (MQ=255)

                                  tGGTGATTGCCGAAGACGACTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTa                                                            >  1:5139785/1‑93 (MQ=255)

                                     tGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCg                                                         >  2:4005633/1‑93 (MQ=255)

                                      gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                                        >  1:4660509/1‑93 (MQ=255)

                                       aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                                       >  2:3631899/1‑93 (MQ=255)

                                        ttGCCGAAGACGCCTGTAGTGCCGCCAGCGCGGAGCAACACAATAACAGCATTAATCATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTgg                                                      <  2:5336076/93‑1 (MQ=255)

                                        ttGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTgg                                                      >  2:2724622/1‑93 (MQ=255)

                                           ccGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                   <  2:4179193/93‑1 (MQ=255)

                                           ccGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                   <  1:4845841/93‑1 (MQ=255)

                                           ccGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                   <  2:4893884/93‑1 (MQ=255)

                                             gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagag                                                 >  1:724012/1‑93 (MQ=255)

                                             gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagaa                                                 >  2:2354944/1‑92 (MQ=255)

                                              aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                >  1:2967166/1‑93 (MQ=255)

                                               aGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGAt                                               >  2:407741/1‑93 (MQ=255)

                                                 aCGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                             <  1:6270562/93‑1 (MQ=255)

                                                 aCGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                             <  1:3894439/93‑1 (MQ=255)

                                                  cGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCt                                            <  1:1810297/93‑1 (MQ=255)

                                                     cTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCaa                                         >  2:546574/1‑93 (MQ=255)

                                                     cTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCaa                                         <  1:948818/93‑1 (MQ=255)

                                                      tGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCACGTGTGCGTAGCGTGGAAGAGATCCTCAAc                                        <  2:463281/93‑1 (MQ=255)

                                                          gTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGt                                    >  1:3092318/1‑93 (MQ=255)

                                                          gTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGt                                    >  2:1396228/1‑93 (MQ=255)

                                                              cgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATg                                <  2:2214703/93‑1 (MQ=255)

                                                               gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGa                               <  2:3208166/93‑1 (MQ=255)

                                                                ccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAt                              >  1:2506423/1‑93 (MQ=255)

                                                                ccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt                              <  1:2745874/93‑1 (MQ=255)

                                                                 cAGCGCCGAGCAACACAATAACAGCATTAATCATATCTGCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAtt                             <  1:4599872/93‑1 (MQ=255)

                                                                  aGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAttt                            >  2:126795/1‑93 (MQ=255)

                                                                   gcgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa                           >  1:3720882/1‑93 (MQ=255)

                                                                      ccGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACAt                        <  1:1020908/93‑1 (MQ=255)

                                                                        gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                      <  2:5576435/93‑1 (MQ=255)

                                                                        gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                      >  2:4585138/1‑93 (MQ=255)

                                                                         agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCgg                     >  1:5723462/1‑93 (MQ=255)

                                                                         agcaacACAATAACCGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCgg                     <  1:569398/93‑1 (MQ=255)

                                                                         agcaacACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGAGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCgg                     >  1:6628/1‑93 (MQ=255)

                                                                          gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt                    <  1:1233045/93‑1 (MQ=255)

                                                                          gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt                    >  1:6073736/1‑93 (MQ=255)

                                                                             acacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtg                 >  2:6252083/2‑92 (MQ=255)

                                                                               acaATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAcc               >  1:5934928/1‑93 (MQ=255)

                                                                                caATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCg              >  2:2428453/1‑93 (MQ=255)

                                                                                 aaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGc             >  2:3677576/1‑93 (MQ=255)

                                                                                  aTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCa            >  2:4358335/1‑93 (MQ=255)

                                                                                  aTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCa            >  2:497085/1‑93 (MQ=255)

                                                                                  aTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa            >  1:4938841/1‑93 (MQ=255)

                                                                                   tAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCaa           >  2:2629727/1‑93 (MQ=255)

                                                                                   tAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCaa           >  1:1569291/1‑93 (MQ=255)

                                                                                    aaCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAAt          <  1:312197/93‑1 (MQ=255)

                                                                                    aaCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAAt          <  2:3936914/93‑1 (MQ=255)

                                                                                     aCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATGGGTCTACCGCAATg         >  2:4709685/1‑93 (MQ=255)

                                                                                     aCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATg         >  1:4011459/1‑93 (MQ=255)

                                                                                     aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATg         <  2:4207854/93‑1 (MQ=255)

ACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT  >  NC_000913/1949267‑1949451

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