BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,369,637	T→C	R345R (CGT→CGC)	mopA →	GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,369,637	0	T	C	99.6%	796.7 / NA	249	R345R (CGT→CGC)	mopA	GroEL
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (121/127); minor base A (0/1); total (121/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

AGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGG  >  NC_000913/4369547‑4369726

aGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCgcgt                                                                                         <  1:2289843/93‑1 (MQ=255)

 gACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGtt                                                                                        >  1:5200998/1‑93 (MQ=255)

   ccTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGc                                                                                      >  1:826674/1‑93 (MQ=255)

    cTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCt                                                                                     >  2:4961629/1‑93 (MQ=255)

    cTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCt                                                                                     >  2:1333400/1‑93 (MQ=255)

     tGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTc                                                                                    >  1:5691366/1‑93 (MQ=255)

     tGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTc                                                                                    >  1:550546/1‑93 (MQ=255)

      gggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCa                                                                                   >  1:4573153/1‑93 (MQ=255)

       ggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAg                                                                                  <  1:1579066/93‑1 (MQ=255)

        gTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGa                                                                                 <  1:5332455/93‑1 (MQ=255)

          cAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATc                                                                               <  1:5662910/93‑1 (MQ=255)

           aGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATcc                                                                              <  1:4127871/93‑1 (MQ=255)

           aGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATcc                                                                              >  1:2383385/1‑93 (MQ=255)

            ggCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCg                                                                             >  2:5732137/1‑93 (MQ=255)

            ggCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCg                                                                             <  2:4473440/93‑1 (MQ=255)

             gCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGt                                                                            >  1:4338844/1‑93 (MQ=255)

              cTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTc                                                                           <  2:2020495/93‑1 (MQ=255)

              cTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTc                                                                           <  2:4747808/93‑1 (MQ=255)

               tAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTca                                                                          >  1:3271583/1‑93 (MQ=255)

                aaaCGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTcag                                                                         >  1:6187259/1‑93 (MQ=255)

                aaaCGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTcag                                                                         >  1:222728/1‑93 (MQ=255)

                aaaCGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTcag                                                                         >  2:276679/1‑93 (MQ=255)

                aaaCGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTcag                                                                         <  1:6192611/93‑1 (MQ=255)

                 aaCGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTcagc                                                                        <  1:2567610/93‑1 (MQ=255)

                    gtgtTGTGATCAACAAAGGCACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGa                                                                     <  2:2172698/93‑1 (MQ=255)

                    gtgtTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGa                                                                     >  2:3034816/1‑93 (MQ=255)

                        tgtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGTTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTga                                                                 <  1:2189930/93‑1 (MQ=255)

                        tgtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTga                                                                 <  2:2377473/93‑1 (MQ=255)

                        tgtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTga                                                                 >  2:5664633/1‑93 (MQ=255)

                        tgtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTga                                                                 <  2:2232026/93‑1 (MQ=255)

                         gtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaa                                                                <  1:6048804/93‑1 (MQ=255)

                         gtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaa                                                                >  1:1257794/1‑93 (MQ=255)

                           gATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaaga                                                              <  1:5875693/93‑1 (MQ=255)

                           gATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTgaaga                                                              >  1:2536972/1‑93 (MQ=255)

                             tCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGCTTGaagaag                                                            >  2:4603897/1‑93 (MQ=255)

                             tCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGaagaag                                                            <  2:2244479/93‑1 (MQ=255)

                               aacaaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCa                                                          <  2:3593169/93‑1 (MQ=255)

                                  aaaGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACt                                                       >  2:148395/1‑93 (MQ=255)

                                   aaGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACtt                                                      >  2:3290112/1‑93 (MQ=255)

                                   aaGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACtt                                                      <  2:346758/93‑1 (MQ=255)

                                    aGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                                                     >  2:4666199/1‑93 (MQ=255)

                                    aGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                                                     <  1:3444183/93‑1 (MQ=255)

                                      accccaccacTATCATCGATGGCGTGGGTTAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTg                                                   <  1:2052492/90‑1 (MQ=255)

                                        accaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                                                 <  2:844653/93‑1 (MQ=255)

                                        accaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                                                 <  1:3940264/93‑1 (MQ=255)

                                          caccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTa                                               <  1:202464/93‑1 (MQ=255)

                                          caccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTa                                               >  1:4799706/1‑93 (MQ=255)

                                          caccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTa                                               <  1:468979/93‑1 (MQ=255)

                                           accacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTAc                                              >  2:3864930/1‑93 (MQ=255)

                                           accacTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTAc                                              >  1:5051996/1‑93 (MQ=255)

                                              acTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGAc                                           >  1:5550787/1‑93 (MQ=255)

                                              acTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGAc                                           >  2:5068533/1‑93 (MQ=255)

                                                tATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCg                                         >  1:1445938/1‑93 (MQ=255)

                                                  tcatcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGCCCGTg                                       >  1:3362322/1‑93 (MQ=255)

                                                  tcatcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTg                                       >  2:1962569/1‑93 (MQ=255)

                                                    atcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaa                                     >  1:2865251/1‑93 (MQ=255)

                                                     tcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaa                                    >  2:5274121/1‑93 (MQ=255)

                                                     tcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaa                                    <  2:6115298/93‑1 (MQ=255)

                                                     tcGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaa                                    >  2:1649549/1‑93 (MQ=255)

                                                      cGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaaa                                   <  1:1080501/93‑1 (MQ=255)

                                                      cGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGaaaa                                   <  1:3000734/93‑1 (MQ=255)

                                                        gtGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAAc                                 <  2:4251703/92‑1 (MQ=255)

                                                         tGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACt                                <  2:6067980/93‑1 (MQ=255)

                                                         tGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACt                                <  2:1734153/93‑1 (MQ=255)

                                                          ggCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTg                               <  2:6252361/93‑1 (MQ=255)

                                                           gCGTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGc                              <  1:3109897/93‑1 (MQ=255)

                                                             gTGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAg                            <  2:4842147/93‑1 (MQ=255)

                                                              tGGGTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAgg                           >  1:37885/1‑93 (MQ=255)

                                                               gggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGa                          <  2:2909792/93‑1 (MQ=255)

                                                               gggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGa                          >  1:1203241/1‑93 (MQ=255)

                                                                ggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGaa                         <  1:2821307/93‑1 (MQ=255)

                                                                ggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGaa                         <  2:4400551/93‑1 (MQ=255)

                                                                ggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGaa                         <  2:4539426/93‑1 (MQ=255)

                                                                ggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGaa                         <  1:4888956/93‑1 (MQ=255)

                                                                ggTGAAGAAGCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGaa                         >  1:578942/1‑93 (MQ=255)

                                                                    aagaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAcgcg                     <  1:964040/93‑1 (MQ=255)

                                                                     agaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGt                    >  1:741205/1‑93 (MQ=255)

                                                                     agaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGt                    <  2:2871327/93‑1 (MQ=255)

                                                                      gaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa                   <  1:1888248/93‑1 (MQ=255)

                                                                      gaagCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa                   >  1:1110037/1‑93 (MQ=255)

                                                                         gCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCg                <  2:5398945/93‑1 (MQ=255)

                                                                         gCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCg                >  1:3619485/1‑93 (MQ=255)

                                                                         gCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCg                >  1:2773587/1‑93 (MQ=255)

                                                                         gCTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCg                >  2:3472056/1‑93 (MQ=255)

                                                                          cTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa               >  1:3045752/1‑93 (MQ=255)

                                                                          cTGCAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa               <  1:374608/93‑1 (MQ=255)

                                                                             cAATCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAAc            <  1:2417833/93‑1 (MQ=255)

                                                                              aaTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACt           >  2:4190947/1‑93 (MQ=255)

                                                                              aaTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACt           >  2:4134796/1‑93 (MQ=255)

                                                                               aTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg          <  2:5898400/93‑1 (MQ=255)

                                                                               aTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg          <  2:3733958/93‑1 (MQ=255)

                                                                               aTCCAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg          >  2:1931332/1‑93 (MQ=255)

                                                                                  cAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCa       <  2:4388613/93‑1 (MQ=255)

                                                                                  cAGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCa       <  1:2269322/93‑1 (MQ=255)

                                                                                   aGGGCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAg      <  2:710347/93‑1 (MQ=255)

                                                                                     ggCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAggc    <  2:3267381/93‑1 (MQ=255)

                                                                                     ggCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAggc    <  2:1366571/93‑1 (MQ=255)

                                                                                      gCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGCCTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAggcg   <  2:4089329/93‑1 (MQ=255)

                                                                                      gCCGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAggcg   >  2:2592412/1‑93 (MQ=255)

                                                                                       ccGCGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAggcgg  <  1:3595696/93‑1 (MQ=255)

AGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGG  >  NC_000913/4369547‑4369726

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