BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,511,632	T→C	Q84R (CAG→CGG)	fecB ←	FecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,511,632	0	T	C	97.1%	515.7 / 0.4	173	Q84R (CAG→CGG)	fecB	FecB
Reads supporting (aligned to +/- strand): ref base T (2/3); new base C (75/93); total (77/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

Only 100 of 198 total aligned reads displayed.

CCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCC  >  NC_000913/4511542‑4511724

 ccGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGcc                                                                                           >  2:2006435/1‑93 (MQ=255)

   gCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCAc                                                                                         >  2:1010250/1‑93 (MQ=255)

    cATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCCCg                                                                                        >  2:2349742/1‑93 (MQ=255)

     aTGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGCCCGCCACgg                                                                                       >  1:3077987/1‑93 (MQ=255)

      tGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGt                                                                                      <  1:1373195/93‑1 (MQ=255)

         cGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGATTCCAGGCTCGGCTGCGCGCGCGGTTCGGCGGACTGTCACGGgttg                                                                                   >  2:3828608/1‑92 (MQ=255)

           aCTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAg                                                                                 <  1:588952/93‑1 (MQ=255)

           aCTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAg                                                                                 >  1:926469/1‑93 (MQ=255)

           aCTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGGTTCAg                                                                                 >  2:1183366/1‑93 (MQ=255)

              gctgTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTg                                                                              >  1:2522012/1‑93 (MQ=255)

              gctgTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTg                                                                              >  2:2247335/1‑93 (MQ=255)

                tgTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTgcg                                                                            >  1:4125003/1‑93 (MQ=255)

                  tCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGcgca                                                                          >  1:2668537/1‑93 (MQ=255)

                     gCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGcgcacgc                                                                       >  2:470989/1‑93 (MQ=255)

                       aaTGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCgcacgcac                                                                     >  2:1118113/1‑93 (MQ=255)

                       aaTGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGAGTTACGACGGACCGCCACGGTTTCAGGTGCgcacgcac                                                                     >  2:851227/1‑93 (MQ=255)

                        aTGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACt                                                                    >  2:3413374/1‑93 (MQ=255)

                        aTGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACt                                                                    <  2:2408575/93‑1 (MQ=255)

                        aTGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACt                                                                    <  1:3017576/93‑1 (MQ=255)

                         tGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACtt                                                                   <  1:422461/93‑1 (MQ=255)

                          gATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTc                                                                  <  1:856363/93‑1 (MQ=255)

                          gATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTc                                                                  <  1:2236883/93‑1 (MQ=255)

                           aTCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCg                                                                 >  2:1649706/1‑93 (MQ=255)

                           aTCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCg                                                                 >  1:4008736/1‑93 (MQ=255)

                           aTCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCg                                                                 >  1:3011263/1‑93 (MQ=255)

                           aTCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCg                                                                 <  2:209712/93‑1 (MQ=255)

                           aTCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGCTGCAGGTGCGCTAGCACTTCg                                                                 >  2:4086420/1‑93 (MQ=255)

                           aTCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCCTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCg                                                                 <  2:3029354/93‑1 (MQ=255)

                            tCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCgg                                                                >  2:205689/1‑93 (MQ=255)

                             cAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCggg                                                               <  1:1605899/93‑1 (MQ=255)

                             cAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCggg                                                               <  1:1925944/93‑1 (MQ=255)

                              aGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGc                                                              <  1:1398151/93‑1 (MQ=255)

                                 tCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAgg                                                           >  1:3656955/1‑93 (MQ=255)

                                 tCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCCGCTGCGAGCGCGTTACGACGGACCGACACGGTTTCAGGTGCGCACGCACTTCGGGCAgg                                                           >  2:3296348/1‑93 (MQ=255)

                                  cTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGa                                                          <  1:3356237/93‑1 (MQ=255)

                                  cTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGa                                                          <  1:1446680/93‑1 (MQ=255)

                                  cTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGa                                                          <  2:2708244/93‑1 (MQ=255)

                                  cTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGa                                                          <  2:881804/93‑1 (MQ=255)

                                   tGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGAt                                                         >  2:3635276/1‑93 (MQ=255)

                                    ggTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATg                                                        <  2:1259990/93‑1 (MQ=255)

                                      tttCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCg                                                      <  1:1185830/93‑1 (MQ=255)

                                      tttCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACCGTTTCAGGTGCGACCGCCCTTACGGCCGGATGCg                                                      >  1:3223531/1‑93 (MQ=255)

                                        tCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGtt                                                    <  2:677288/93‑1 (MQ=255)

                                        tCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGtt                                                    <  2:1331713/93‑1 (MQ=255)

                                         cAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGttt                                                   >  1:2728949/1‑93 (MQ=255)

                                          agagCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGtttt                                                  <  2:30904/93‑1 (MQ=255)

                                           gagCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTg                                                 >  2:3785260/1‑93 (MQ=255)

                                             gCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCa                                               <  2:4524455/93‑1 (MQ=255)

                                             gCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGCGCGCACGCACTTCGGGCAGGATGCGTTTTGCa                                               >  1:2047645/1‑93 (MQ=255)

                                              cGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCAt                                              <  2:1632532/93‑1 (MQ=255)

                                               ggCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATc                                             <  1:1829253/93‑1 (MQ=255)

                                                gCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCg                                            <  1:2323852/93‑1 (MQ=255)

                                                gCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCg                                            >  2:1711612/1‑93 (MQ=255)

                                                 cAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGt                                           >  1:2117295/1‑93 (MQ=255)

                                                 cAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGt                                           >  2:1056143/1‑93 (MQ=255)

                                                 cAATGGCTGCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACCTCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGt                                           <  1:1865380/93‑1 (MQ=255)

                                                  aaTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGtt                                          <  1:1594524/93‑1 (MQ=255)

                                                   aTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGTAGGATGCGTTTTGCATCGttt                                         >  2:4186874/1‑92 (MQ=255)

                                                   aTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTa                                         >  2:2416204/1‑93 (MQ=255)

                                                   aTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTa                                         <  2:175139/93‑1 (MQ=255)

                                                     ggCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAgg                                                                                >  1:1151322/1‑52 (MQ=255)

                                                     ggCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtc                                       <  1:447649/93‑1 (MQ=255)

                                                        ttCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgtc                                    >  1:2800494/1‑93 (MQ=255)

                                                         tCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgtcg                                   <  1:1495089/93‑1 (MQ=255)

                                                         tCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgtcg                                   <  2:2785945/93‑1 (MQ=255)

                                                           cAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGc                                 <  2:2808791/93‑1 (MQ=255)

                                                            aGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCa                                >  1:4313763/1‑93 (MQ=255)

                                                            aGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCa                                >  2:3181862/1‑93 (MQ=255)

                                                            aGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCa                                <  2:2504701/93‑1 (MQ=255)

                                                             ggCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCaa                               <  2:282807/93‑1 (MQ=255)

                                                             ggCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCaa                               >  1:4025833/1‑93 (MQ=255)

                                                             ggCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCaa                               >  1:2972476/1‑93 (MQ=255)

                                                               cTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATa                             <  1:1263508/93‑1 (MQ=255)

                                                                     tGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATc                       >  2:767529/1‑93 (MQ=255)

                                                                     tGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATc                       <  2:339346/93‑1 (MQ=255)

                                                                       cgcgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCgg                     <  2:4426766/93‑1 (MQ=255)

                                                                       cgcgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCgg                     >  1:4630197/1‑93 (MQ=255)

                                                                        gcgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCggg                    >  2:537583/1‑93 (MQ=255)

                                                                        gcgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCggg                    <  2:2838772/93‑1 (MQ=255)

                                                                         cgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGa                   <  1:1416215/93‑1 (MQ=255)

                                                                         cgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGa                   <  1:4413648/93‑1 (MQ=255)

                                                                          gcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGAt                  <  2:2386773/93‑1 (MQ=255)

                                                                           cgcgTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATg                 >  1:4511317/1‑93 (MQ=255)

                                                                           cgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATg                 >  2:1447937/1‑93 (MQ=255)

                                                                            gcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATgc                >  2:941055/1‑92 (MQ=255)

                                                                            gcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATgc                >  2:4394234/1‑92 (MQ=255)

                                                                               ttCCGAGAGACCACCACGGTTTCAGCTGAGCACGCTTTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGt             <  2:4135112/93‑1 (MQ=255)

                                                                               ttCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACga             >  2:3162688/1‑92 (MQ=255)

                                                                               ttCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGt             <  2:3441503/93‑1 (MQ=255)

                                                                               ttCCGACGGACCGCCAAGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGt             <  1:4198079/93‑1 (MQ=255)

                                                                                tCCGACGGACCTCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTc            <  2:1958817/93‑1 (MQ=255)

                                                                                 ccGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTcc           <  2:230499/93‑1 (MQ=255)

                                                                                 ccGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTcc           <  2:154539/93‑1 (MQ=255)

                                                                                  cGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTcgc          >  1:3840972/1‑91 (MQ=255)

                                                                                  cGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCa          >  1:3562986/1‑93 (MQ=255)

                                                                                   gACGGGCCGCCACGGTTTCAGGTGCGCATGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAc         <  2:2349738/93‑1 (MQ=255)

                                                                                     cGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAcgg       <  1:72684/93‑1 (MQ=255)

                                                                                       gACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAcggcg     <  2:4613455/93‑1 (MQ=255)

                                                                                       gACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAcggcg     >  1:2483579/1‑93 (MQ=255)

                                                                                          tGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGcc  <  1:3808423/93‑1 (MQ=255)

CCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCC  >  NC_000913/4511542‑4511724

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