BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,576,170	A→G	R246R (CGT→CGC)	mcrB ←	McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,576,170	0	A	G	99.4%	592.2 / ‑4.9	167	R246R (CGT→CGC)	mcrB	McrB
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (117/49); total (118/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 168 total aligned reads displayed.

AAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCA  >  NC_000913/4576080‑4576258

aaaaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCg                                                                                        >  2:6034886/1‑93 (MQ=255)

 aaaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGa                                                                                       >  1:2396485/1‑93 (MQ=255)

  aaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGAt                                                                                      >  2:316157/1‑93 (MQ=255)

   aaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATa                                                                                     >  1:4459793/1‑93 (MQ=255)

   aaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATa                                                                                     >  2:6135508/1‑93 (MQ=255)

    atatACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAg                                                                                    >  1:156271/1‑93 (MQ=255)

    atatACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAg                                                                                    >  2:3074805/1‑93 (MQ=255)

     tataCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGc                                                                                   >  1:2926424/1‑93 (MQ=255)

     tataCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGc                                                                                   >  1:937302/1‑93 (MQ=255)

     tataCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGc                                                                                   >  2:1258026/1‑93 (MQ=255)

       taCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGccc                                                                                 >  1:2465989/1‑93 (MQ=255)

       taCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGccc                                                                                 <  2:3678709/93‑1 (MQ=255)

       taCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGccc                                                                                 >  2:4791972/1‑93 (MQ=255)

        aCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCt                                                                                <  1:5588499/93‑1 (MQ=255)

            tttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAATCGGGCGATAGCCCTGtat                                                                            >  2:4752224/1‑93 (MQ=255)

            tgtctcTGGCTGCTCTTTAGCGTGCTGACAAAAATTCTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGtat                                                                            <  1:2857616/91‑1 (MQ=255)

             ttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGtata                                                                           >  1:2171624/1‑93 (MQ=255)

             ttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGtata                                                                           >  1:6085957/1‑93 (MQ=255)

              tctcTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaa                                                                          >  2:4527895/1‑93 (MQ=255)

              tctcTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaa                                                                          >  2:5838199/1‑93 (MQ=255)

              tctcTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaa                                                                          >  1:1138048/1‑93 (MQ=255)

               ctctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaaa                                                                         >  1:3984815/1‑93 (MQ=255)

               ctctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaaa                                                                         >  1:5298252/1‑93 (MQ=255)

                tctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaaaa                                                                        <  1:4417754/93‑1 (MQ=255)

                 ctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAAt                                                                       <  1:1036991/93‑1 (MQ=255)

                  tGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATc                                                                      >  2:715728/1‑93 (MQ=255)

                  tGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATc                                                                      >  1:898939/1‑93 (MQ=255)

                   ggCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATcc                                                                     >  2:1810428/1‑93 (MQ=255)

                   ggCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATcc                                                                     <  2:5343937/93‑1 (MQ=255)

                   ggCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGAAGCCACTCGGGCGATAGCCCTGTATAAAATcc                                                                     >  1:3497554/1‑93 (MQ=255)

                    gctgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCt                                                                    >  2:1763365/1‑93 (MQ=255)

                    gctgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCt                                                                    >  2:558145/1‑93 (MQ=255)

                     ctgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTc                                                                   >  1:4155922/1‑93 (MQ=255)

                     ctgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTc                                                                   <  1:5227596/93‑1 (MQ=255)

                      tgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATCGCCCTGTATAAAATCCTCa                                                                  >  2:152139/1‑93 (MQ=255)

                      tgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCa                                                                  >  2:500280/1‑93 (MQ=255)

                      tgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCa                                                                  >  2:381523/1‑93 (MQ=255)

                      tgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCa                                                                  >  1:3162352/1‑93 (MQ=255)

                       gctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGGTAGCCCTGTATAAAATCCTCAt                                                                 >  1:1283865/1‑93 (MQ=255)

                       gctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCAt                                                                 >  1:461654/1‑93 (MQ=255)

                       gctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCAt                                                                 <  2:2995205/93‑1 (MQ=255)

                       gctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCAt                                                                 >  1:3816926/1‑93 (MQ=255)

                       gctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCAt                                                                 <  1:3532987/93‑1 (MQ=255)

                        ctctTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATa                                                                <  2:2457934/93‑1 (MQ=255)

                         tctTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAg                                                               <  2:3702290/93‑1 (MQ=255)

                             tAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCtat                                                           >  1:4776907/1‑93 (MQ=255)

                               gCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATaa                                                         >  2:5079326/1‑93 (MQ=255)

                               gCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATaa                                                         <  2:1813529/93‑1 (MQ=255)

                               gCTTGCTGACAAAAATTGTAAAATATGCCGCCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATaa                                                         >  2:4411539/1‑93 (MQ=255)

                                 ttGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGa                                                       >  2:4977902/1‑93 (MQ=255)

                                 ttGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGa                                                       <  2:4287699/93‑1 (MQ=255)

                                   gCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGAtt                                                     >  1:4370264/1‑93 (MQ=255)

                                   gCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGAtt                                                     >  2:61698/1‑93 (MQ=255)

                                    cTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCtata                                                          >  1:3694572/1‑87 (MQ=255)

                                    cTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTg                                                    >  2:2609125/1‑93 (MQ=255)

                                      gACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt                                                  >  1:2821652/1‑93 (MQ=255)

                                      gACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt                                                  >  1:5012869/1‑93 (MQ=255)

                                      gACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt                                                  >  1:5665661/1‑93 (MQ=255)

                                        cAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATgg                                                >  1:5789625/1‑93 (MQ=255)

                                          aaaaTTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGaa                                              >  2:762504/1‑93 (MQ=255)

                                            aaTTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACt                                            <  1:2615243/93‑1 (MQ=255)

                                                gTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAAc                                        >  2:2966726/1‑93 (MQ=255)

                                                  aaaaTATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCa                                      >  1:2435597/1‑93 (MQ=255)

                                                   aaaTATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCat                                     >  1:5054874/1‑93 (MQ=255)

                                                     atatGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCatat                                   <  1:1237755/93‑1 (MQ=255)

                                                      tatGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATAtt                                  <  2:4303917/93‑1 (MQ=255)

                                                        tGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGa                                <  2:3815103/93‑1 (MQ=255)

                                                        tGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGa                                >  2:3573256/1‑93 (MQ=255)

                                                        tGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGa                                >  2:3513654/1‑93 (MQ=255)

                                                          ccGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcg                              <  1:3358362/93‑1 (MQ=255)

                                                           cGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgc                             >  1:5099821/1‑93 (MQ=255)

                                                           cGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgc                             >  2:3186502/1‑93 (MQ=255)

                                                            gTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgcg                            >  1:5874458/1‑93 (MQ=255)

                                                             tCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGt                           >  1:198588/1‑93 (MQ=255)

                                                             tCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGt                           <  2:4838089/93‑1 (MQ=255)

                                                              cTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGtt                          >  1:2483148/1‑93 (MQ=255)

                                                              cTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGtt                          >  1:2298092/1‑93 (MQ=255)

                                                              cTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGtt                          <  2:3026970/93‑1 (MQ=255)

                                                                ttACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGc                        <  1:1978903/93‑1 (MQ=255)

                                                                 tACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCg                       >  2:4122655/1‑93 (MQ=255)

                                                                  aCGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCgg                      >  1:3052915/1‑93 (MQ=255)

                                                                    gTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAg                    >  2:3917688/1‑93 (MQ=255)

                                                                       ggAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCt                 <  2:5851305/93‑1 (MQ=255)

                                                                       ggAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCt                 >  2:458662/1‑93 (MQ=255)

                                                                        gAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtt                >  2:1624246/1‑93 (MQ=255)

                                                                          aGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtttt              >  1:3225996/1‑93 (MQ=255)

                                                                          aGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtttt              >  2:3237522/1‑93 (MQ=255)

                                                                          aGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCAAATTGCCGCGTGGGGGAGCCtttt              >  1:548900/1‑93 (MQ=255)

                                                                           gCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTtc             >  1:3275296/1‑93 (MQ=255)

                                                                            ccGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTtct            >  1:4195099/1‑93 (MQ=255)

                                                                              gACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTcc          >  2:556354/1‑93 (MQ=255)

                                                                              gACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTcc          <  1:2584838/93‑1 (MQ=255)

                                                                              gACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTcc          <  2:2365484/93‑1 (MQ=255)

                                                                               aCGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCt         >  2:4571811/1‑93 (MQ=255)

                                                                                cGCCATTCGGGCGATGGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTg        >  1:2278248/1‑93 (MQ=255)

                                                                                cGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTg        >  1:1027829/1‑93 (MQ=255)

                                                                                 gCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGt       >  1:4798128/1‑93 (MQ=255)

                                                                                  ccATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTc      >  1:284816/1‑93 (MQ=255)

                                                                                    aTTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAg    >  2:3278864/1‑93 (MQ=255)

                                                                                     ttCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGc   >  2:1907514/1‑93 (MQ=255)

AAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCA  >  NC_000913/4576080‑4576258

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