BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,197,191	T→C	Q31R (CAG→CGG)	glnE ←	GlnE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,197,191	0	T	C	98.0%	856.7 / 2.7	247	Q31R (CAG→CGG)	glnE	GlnE
Reads supporting (aligned to +/- strand): ref base T (2/3); new base C (113/129); total (115/132)
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.97e-01

Only 100 of 254 total aligned reads displayed.

GGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCGGCTTC  >  NC_000913/3197101‑3197280

ggcggTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCcgt                                                                                         >  1:2010549/1‑93 (MQ=255)

 gcggTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCcgtg                                                                                        >  2:3655027/1‑93 (MQ=255)

  cggTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTgc                                                                                       <  1:4351174/93‑1 (MQ=255)

    gTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTgcgc                                                                                     <  2:1190099/93‑1 (MQ=255)

    gTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTgcgc                                                                                     >  2:5208054/1‑93 (MQ=255)

     ttGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCt                                                                                    >  1:3140081/1‑93 (MQ=255)

     ttGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCt                                                                                    >  2:3986116/1‑93 (MQ=255)

     ttGCCTTTCCAGTTCCCTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCt                                                                                    <  1:3376062/93‑1 (MQ=255)

      tGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTa                                                                                   >  1:2956799/1‑93 (MQ=255)

       ggCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTaa                                                                                  >  2:5567247/1‑93 (MQ=255)

       ggCTTCCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTaa                                                                                  <  2:2556528/93‑1 (MQ=255)

        gCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAg                                                                                 >  2:3139285/1‑93 (MQ=255)

        gCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAg                                                                                 <  1:316371/93‑1 (MQ=255)

         cTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGt                                                                                >  1:5500883/1‑93 (MQ=255)

          tttCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTg                                                                               >  1:1965462/1‑93 (MQ=255)

           ttCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGa                                                                              <  2:3696305/93‑1 (MQ=255)

           ttCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGa                                                                              >  1:4408591/1‑93 (MQ=255)

           ttCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGa                                                                              <  2:5469095/93‑1 (MQ=255)

             ccAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGAtt                                                                            <  1:1168122/93‑1 (MQ=255)

             ccAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGAtt                                                                            <  2:6119172/93‑1 (MQ=255)

             ccAGTCCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGAtt                                                                            <  1:3069551/93‑1 (MQ=255)

              cAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTc                                                                           <  2:854968/93‑1 (MQ=255)

              cAGTCCCGCCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTc                                                                           <  2:2232241/93‑1 (MQ=255)

                gtcccGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCt                                                                         <  2:4368301/90‑1 (MQ=255)

                gTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCt                                                                         >  2:3935087/1‑93 (MQ=255)

                gTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCt                                                                         >  1:379472/1‑93 (MQ=255)

                gTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCt                                                                         <  1:4553038/93‑1 (MQ=255)

                gTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCt                                                                         >  2:2214779/1‑93 (MQ=255)

                   ccGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCgg                                                                      <  1:943229/93‑1 (MQ=255)

                        aGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTaaa                                                                 <  1:5113022/93‑1 (MQ=255)

                         gCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAg                                                                >  2:5579717/1‑93 (MQ=255)

                           cACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGc                                                              >  1:706635/1‑93 (MQ=255)

                            cctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGct                                                             <  1:5206533/92‑1 (MQ=255)

                             ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                                            <  1:3255383/93‑1 (MQ=255)

                             ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                                            <  2:3244213/93‑1 (MQ=255)

                             ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                                            >  1:3601351/1‑93 (MQ=255)

                              tctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGctct                                                           >  2:5262795/1‑93 (MQ=255)

                               ctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTg                                                          <  1:1090038/93‑1 (MQ=255)

                               ctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTg                                                          <  1:1281061/93‑1 (MQ=255)

                                     gcgcAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGcc                                                    <  1:1007345/93‑1 (MQ=255)

                                        cAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCt                                                 >  2:1273166/1‑93 (MQ=255)

                                        cAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCt                                                 <  1:201183/93‑1 (MQ=255)

                                           tCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcac                                              >  2:1256839/1‑92 (MQ=255)

                                            cACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaac                                             <  1:2269451/93‑1 (MQ=255)

                                             aCGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaaca                                            >  2:2407675/1‑93 (MQ=255)

                                              cGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaacaa                                           >  1:6031642/1‑93 (MQ=255)

                                              cGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaacaa                                           <  2:3537360/93‑1 (MQ=255)

                                               gCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCaacaac                                          >  1:2688879/1‑93 (MQ=255)

                                                 tGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACgg                                        <  2:2036252/93‑1 (MQ=255)

                                                 tGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACgg                                        >  2:1167995/1‑93 (MQ=255)

                                                 tGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACgg                                        <  1:2382830/93‑1 (MQ=255)

                                                  gTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGt                                       >  1:5349582/1‑93 (MQ=255)

                                                   tCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTc                                      <  2:2741635/93‑1 (MQ=255)

                                                   tCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTc                                      <  2:1378764/93‑1 (MQ=255)

                                                    ccTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCt                                     >  2:3462146/1‑93 (MQ=255)

                                                    ccTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGACt                                     >  2:1088838/1‑93 (MQ=255)

                                                    ccTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCt                                     >  2:3038609/1‑93 (MQ=255)

                                                      tGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGc                                   >  1:1352709/1‑93 (MQ=255)

                                                         acaAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAg                                <  1:413264/93‑1 (MQ=255)

                                                         acaAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAg                                >  1:1856504/1‑93 (MQ=255)

                                                          caAAATCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGt                               >  2:2334821/1‑93 (MQ=255)

                                                           aaaaTCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTa                              >  2:5003056/1‑93 (MQ=255)

                                                            aaaTCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAc                             <  2:1658750/93‑1 (MQ=255)

                                                            aaaTCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAc                             <  1:3286071/93‑1 (MQ=255)

                                                             aaTCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAct                            >  1:2637668/1‑93 (MQ=255)

                                                              aTCACTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActg                           >  1:1114347/1‑93 (MQ=255)

                                                                 aCTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgctg                        <  2:3830164/93‑1 (MQ=255)

                                                                  cTAAAAGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGt                       >  2:2001815/1‑93 (MQ=255)

                                                                     aaaGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAAc                    >  2:448596/1‑93 (MQ=255)

                                                                     aaaGTAAGTACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAAc                    <  1:3394247/93‑1 (MQ=255)

                                                                      aagtaagtACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACg                   >  2:4735054/1‑93 (MQ=255)

                                                                      aagtaagtACTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACg                   <  1:4778330/93‑1 (MQ=255)

                                                                       agtaagtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACgg                  >  1:2503907/1‑93 (MQ=255)

                                                                        gtaagtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGt                 >  2:5896508/1‑93 (MQ=255)

                                                                        gtaagtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGt                 <  1:123307/93‑1 (MQ=255)

                                                                        ataagtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGt                 <  2:4589459/92‑1 (MQ=255)

                                                                         taagtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTg                >  1:5264533/1‑93 (MQ=255)

                                                                          aagtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGa               <  1:4821860/93‑1 (MQ=255)

                                                                          aagtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGa               <  1:5544083/93‑1 (MQ=255)

                                                                          aagtaCTGACTTCGCCCGTGCGCTAAGTGACTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGa               <  1:871107/93‑1 (MQ=255)

                                                                           agtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGaa              <  1:4589862/93‑1 (MQ=255)

                                                                           agtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGaa              <  1:559077/93‑1 (MQ=255)

                                                                            gtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAag             <  1:5838568/93‑1 (MQ=255)

                                                                            gtaCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAag             >  2:6073730/1‑93 (MQ=255)

                                                                             taCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGATCTGGCAGCCGCTCAACACCGGTCTGCCAGTACTGCTGTCACGGTGAagg            >  2:5820041/1‑92 (MQ=255)

                                                                              aCTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAagag           >  2:4569026/1‑93 (MQ=255)

                                                                               cTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAagaga          >  2:284223/1‑93 (MQ=255)

                                                                               cTGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAagaga          >  2:2442372/1‑93 (MQ=255)

                                                                                tGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAagagag         <  2:2575103/93‑1 (MQ=255)

                                                                                tGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAagagag         <  1:608510/93‑1 (MQ=255)

                                                                                tGACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAagagag         <  2:3087555/93‑1 (MQ=255)

                                                                                 gACTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGc        <  2:1559083/93‑1 (MQ=255)

                                                                                  aCTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCg       >  2:5446799/1‑93 (MQ=255)

                                                                                  aCTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCg       >  1:5300142/1‑93 (MQ=255)

                                                                                   cTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCgg      >  2:4865816/1‑93 (MQ=255)

                                                                                   cTTCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCgg      >  2:3430815/1‑93 (MQ=255)

                                                                                    ttCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCGGc     >  2:4268056/1‑93 (MQ=255)

                                                                                    ttCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCGGc     <  1:5003927/93‑1 (MQ=255)

                                                                                    ttCGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCGGc     >  1:4867526/1‑93 (MQ=255)

                                                                                      cGCCCGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCGGCtt   >  1:1574802/1‑93 (MQ=255)

GGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGAAGAGAGCGGCTTC  >  NC_000913/3197101‑3197280

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