BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,016,796	T→C	Y119H (TAT→CAT)	ubiE →	UbiE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,016,796	0	T	C	99.6%	954.9 / NA	272	Y119H (TAT→CAT)	ubiE	UbiE
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (144/127); total (145/127)
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.99e-01

Only 100 of 272 total aligned reads displayed.

GGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCCCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGAC  >  NC_000913/4016706‑4016885

  cAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGt                                                                                       <  1:3840154/93‑1 (MQ=255)

  cAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGt                                                                                       >  2:4539247/1‑93 (MQ=255)

  cAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGt                                                                                       >  2:4364089/1‑93 (MQ=255)

     aGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCa                                                                                    >  2:703867/1‑93 (MQ=255)

        ggTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGc                                                                                 >  1:639001/1‑93 (MQ=255)

        ggTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGc                                                                                 >  1:3257097/1‑93 (MQ=255)

        ggTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGc                                                                                 >  2:5474036/1‑93 (MQ=255)

        ggTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGc                                                                                 <  2:1358689/93‑1 (MQ=255)

        ggTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGc                                                                                 <  1:4968826/93‑1 (MQ=255)

         gTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCg                                                                                <  2:971991/93‑1 (MQ=255)

         gTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCg                                                                                <  1:5464069/93‑1 (MQ=255)

          tCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGa                                                                               >  1:1508271/1‑93 (MQ=255)

          tCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGa                                                                               <  2:4643100/93‑1 (MQ=255)

           ccTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGaa                                                                              >  1:3706357/1‑93 (MQ=255)

            cTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAAc                                                                             <  1:5771030/93‑1 (MQ=255)

                cTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTg                                                                         >  2:4443694/1‑93 (MQ=255)

                  gATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCCTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAg                                                                       >  1:533045/1‑93 (MQ=255)

                   atatCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAgg                                                                      <  1:514710/93‑1 (MQ=255)

                   atatCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGGTCAGGCGAACGCTGAgg                                                                      >  1:1836405/1‑93 (MQ=255)

                    tatCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGgc                                                                     <  1:2887846/93‑1 (MQ=255)

                    tatCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGgc                                                                     >  2:5894621/1‑93 (MQ=255)

                    tatCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGgc                                                                     <  2:1146915/93‑1 (MQ=255)

                      tCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGgcgc                                                                   >  1:5082923/1‑93 (MQ=255)

                      tCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGgcgc                                                                   >  1:1000375/1‑93 (MQ=255)

                      tCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGgcgc                                                                   <  1:2038499/93‑1 (MQ=255)

                        aaTGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTg                                                                 <  2:3181453/93‑1 (MQ=255)

                        aaTGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTg                                                                 >  2:3136805/1‑93 (MQ=255)

                         aTGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCATGCGAACGCTGAGGCGCTGc                                                                <  1:6239046/93‑1 (MQ=255)

                         aTGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGc                                                                <  2:4123721/93‑1 (MQ=255)

                          tGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGcc                                                               >  2:2078959/1‑93 (MQ=255)

                          tGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGcc                                                               <  1:2974432/93‑1 (MQ=255)

                            aaTCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAAGGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGt                                                             <  1:6150960/93‑1 (MQ=255)

                            aaTCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGt                                                             >  1:5693558/1‑93 (MQ=255)

                             aTCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGtt                                                            >  2:2879064/1‑93 (MQ=255)

                              tCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTc                                                           >  1:3141720/1‑93 (MQ=255)

                                cATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTccc                                                         >  1:2510468/1‑93 (MQ=255)

                                cATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTccc                                                         >  1:1677721/1‑93 (MQ=255)

                                cATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTccc                                                         <  1:2735496/93‑1 (MQ=255)

                                cATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTccc                                                         <  2:2701396/93‑1 (MQ=255)

                                 aTGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCg                                                        >  2:1278601/1‑93 (MQ=255)

                                   gctcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGa                                                      <  1:1262251/90‑1 (MQ=255)

                                   gctcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGa                                                      <  1:81949/90‑1 (MQ=255)

                                    ctcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGAt                                                     <  2:420662/91‑1 (MQ=255)

                                     tcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATa                                                    <  1:2355189/92‑1 (MQ=255)

                                     tcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATa                                                    <  2:4163732/92‑1 (MQ=255)

                                      cAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATaa                                                   <  1:785695/93‑1 (MQ=255)

                                        aaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAaca                                                 >  2:6002940/1‑93 (MQ=255)

                                        aaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAaca                                                 >  1:1718459/1‑93 (MQ=255)

                                         aaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAacac                                                >  1:4289670/1‑93 (MQ=255)

                                          aTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACAcc                                               >  2:5830413/1‑93 (MQ=255)

                                          aTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACAcc                                               >  2:3565427/1‑93 (MQ=255)

                                           tGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCt                                              >  2:4827994/1‑93 (MQ=255)

                                            gggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCtt                                             <  2:5182152/93‑1 (MQ=255)

                                            gggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCtt                                             <  1:4799415/93‑1 (MQ=255)

                                             ggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCttt                                            <  1:1574430/93‑1 (MQ=255)

                                             ggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCttt                                            <  1:1993329/93‑1 (MQ=255)

                                               ccGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGa                                          >  1:2083191/1‑93 (MQ=255)

                                               ccGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGa                                          <  2:2172884/93‑1 (MQ=255)

                                                cgcgAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGAt                                         <  1:229004/93‑1 (MQ=255)

                                                cgcgAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGAt                                         >  2:2289071/1‑93 (MQ=255)

                                                   gagaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGc                                      >  1:4114310/1‑93 (MQ=255)

                                                    agaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCa                                     >  2:5361033/1‑93 (MQ=255)

                                                    agaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCa                                     >  1:4157363/1‑93 (MQ=255)

                                                    agaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCa                                     >  2:3232014/1‑93 (MQ=255)

                                                     gaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGGTTGCAt                                    >  1:2241351/1‑93 (MQ=255)

                                                     gaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCAt                                    <  1:3930577/93‑1 (MQ=255)

                                                      aaGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATc                                   >  2:2513272/1‑93 (MQ=255)

                                                        gCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCAc                                 >  2:4492083/1‑93 (MQ=255)

                                                         cTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCAcc                                <  2:198946/93‑1 (MQ=255)

                                                         cTGCGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCAcc                                <  2:139462/93‑1 (MQ=255)

                                                            cGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCAtt                             >  2:4703964/1‑93 (MQ=255)

                                                            cGTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCAtt                             >  1:357519/1‑93 (MQ=255)

                                                             gTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                            >  1:5947560/1‑93 (MQ=255)

                                                             gTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                            >  2:885868/1‑93 (MQ=255)

                                                             gTAATATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                            >  2:1675757/1‑93 (MQ=255)

                                                               aaTATCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCg                          >  1:4667649/1‑93 (MQ=255)

                                                                 tatCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGtt                        <  2:4853340/93‑1 (MQ=255)

                                                                  atCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                       <  2:6193602/93‑1 (MQ=255)

                                                                   tCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTg                      >  2:46069/1‑93 (MQ=255)

                                                                   tCGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTg                      >  1:329581/1‑93 (MQ=255)

                                                                    cGGTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTgg                     >  1:4398422/1‑93 (MQ=255)

                                                                     ggTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGt                    <  2:5713827/93‑1 (MQ=255)

                                                                     ggTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGt                    <  2:6403227/93‑1 (MQ=255)

                                                                     ggTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGt                    >  1:520166/1‑93 (MQ=255)

                                                                     ggTGTGATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGt                    <  2:6349021/93‑1 (MQ=255)

                                                                        gtgATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTg                 >  1:6190032/1‑93 (MQ=255)

                                                                        gtgATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTg                 <  2:783294/93‑1 (MQ=255)

                                                                         tgATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                <  2:1785023/93‑1 (MQ=255)

                                                                          gATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCg               >  2:6237288/1‑93 (MQ=255)

                                                                          gATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCg               >  2:3949458/1‑93 (MQ=255)

                                                                          gATTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGATGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCg               >  1:749404/1‑93 (MQ=255)

                                                                           aTTGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGt              >  1:1012441/1‑93 (MQ=255)

                                                                             tGGCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTaa            <  2:5550766/93‑1 (MQ=255)

                                                                              ggCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAAc           <  2:2015616/93‑1 (MQ=255)

                                                                              ggCAACGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAAc           >  2:1526938/1‑93 (MQ=255)

                                                                                  aCGTTGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTTGTCTGCGTAACGTCa       >  1:4941219/1‑93 (MQ=255)

                                                                                     ttGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCg    <  2:4766493/93‑1 (MQ=255)

                                                                                     ttGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCg    >  2:1222251/1‑93 (MQ=255)

                                                                                     ttGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCg    >  1:1943219/1‑93 (MQ=255)

                                                                                      tGAGCATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGa   <  2:3405629/93‑1 (MQ=255)

GGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCCCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGAC  >  NC_000913/4016706‑4016885

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