BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,016,822	G→A	L127L (CTG→CTA)	ubiE →	UbiE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,016,822	0	G	A	99.6%	996.7 / NA	281	L127L (CTG→CTA)	ubiE	UbiE
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (133/147); minor base C (1/0); total (134/147)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

TGAATCCATGCCCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATGTA  >  NC_000913/4016732‑4016911

  aaTCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGt                                                                                       <  1:1725315/93‑1 (MQ=255)

  aaTCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGt                                                                                       >  1:2498729/1‑93 (MQ=255)

    tCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTc                                                                                     <  1:1609569/93‑1 (MQ=255)

     ccATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTcc                                                                                    <  2:1903875/93‑1 (MQ=255)

      cATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTccc                                                                                   >  1:1416164/1‑93 (MQ=255)

       aTGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTGCCCg                                                                                  <  2:2541969/93‑1 (MQ=255)

        tGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCgg                                                                                 >  1:976403/1‑93 (MQ=255)

         gctcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGa                                                                                <  2:2362085/90‑1 (MQ=255)

            cAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATaa                                                                             <  1:1781314/93‑1 (MQ=255)

            cAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATaa                                                                             <  2:3869509/93‑1 (MQ=255)

             aaaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAac                                                                            <  2:3199504/93‑1 (MQ=255)

             aaaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAac                                                                            >  1:4308445/1‑93 (MQ=255)

              aaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAacc                                                                           >  1:640282/1‑92 (MQ=255)

              aaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAaca                                                                           >  1:1643402/1‑93 (MQ=255)

               aaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAacac                                                                          >  1:826768/1‑93 (MQ=255)

                 tGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCt                                                                        <  2:3304834/93‑1 (MQ=255)

                  gggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCTCGGATAACACCtt                                                                       >  1:2676748/1‑93 (MQ=255)

                  gggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCtt                                                                       <  2:4455804/93‑1 (MQ=255)

                     ccGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGa                                                                    <  2:4229211/93‑1 (MQ=255)

                     ccGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGa                                                                    <  1:2954346/93‑1 (MQ=255)

                       gcgAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGAtt                                                                  >  1:1498057/1‑93 (MQ=255)

                        cgAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTg                                                                 >  1:4273765/1‑93 (MQ=255)

                        cgAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTg                                                                 <  1:3459543/93‑1 (MQ=255)

                        cgAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGAGAACGCTGAGGCGCTACCGTTCCCGGATAACACATTTGATTg                                                                 >  2:1687331/1‑93 (MQ=255)

                         gagaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGc                                                                >  1:2811004/1‑93 (MQ=255)

                         gagaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGc                                                                <  2:2458531/93‑1 (MQ=255)

                         gagaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGc                                                                >  2:4343711/1‑93 (MQ=255)

                           gaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCAt                                                              >  2:222283/1‑93 (MQ=255)

                           gaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCAt                                                              <  2:680720/93‑1 (MQ=255)

                           gaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCAt                                                              >  2:2627729/1‑93 (MQ=255)

                            aaGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATc                                                             <  2:1154821/93‑1 (MQ=255)

                            aaGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATc                                                             >  1:4683761/1‑93 (MQ=255)

                            aaGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATc                                                             <  1:1533740/93‑1 (MQ=255)

                              gCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCAc                                                           <  2:1081547/93‑1 (MQ=255)

                               cTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCAcc                                                          <  1:52270/93‑1 (MQ=255)

                                tGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCa                                                         <  2:4048963/93‑1 (MQ=255)

                                tGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCa                                                         >  1:1877428/1‑93 (MQ=255)

                                 gCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAt                                                        <  1:3795396/93‑1 (MQ=255)

                                 gCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAt                                                        <  1:2482524/93‑1 (MQ=255)

                                   gTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                                                      >  1:4515478/1‑93 (MQ=255)

                                   gTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                                                      >  1:4126035/1‑93 (MQ=255)

                                   gTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                                                      <  1:2921670/93‑1 (MQ=255)

                                   gTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATCGCATCACCAttt                                                      >  2:2887832/1‑93 (MQ=255)

                                    tattatCGGTGTGATTGACAGCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACTATTCCCGGCTAACGCCTTTGATTGCATCACCATTTc                                                     <  2:810624/90‑1 (MQ=255)

                                    tAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTc                                                     <  2:2836229/93‑1 (MQ=255)

                                      atatCGGTGTGATTGTCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGt                                                   <  1:1971352/93‑1 (MQ=255)

                                      atatCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGt                                                   >  1:1397442/1‑93 (MQ=255)

                                        atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                                                 <  1:3631907/93‑1 (MQ=255)

                                         tCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTg                                                <  2:3731868/93‑1 (MQ=255)

                                         tCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTg                                                <  1:1296559/93‑1 (MQ=255)

                                            gtgtGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTc                                             >  2:19300/1‑93 (MQ=255)

                                               tgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                                          >  2:1418930/1‑93 (MQ=255)

                                               tgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                                          >  2:574266/1‑93 (MQ=255)

                                               tgATTGGAAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCTCTACCGTTCCCGGATAACGCCTTTGATTGCATCACCATTTCGTTTGGTCTGc                                          <  2:209896/93‑1 (MQ=255)

                                                 aTTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGt                                        <  2:890180/93‑1 (MQ=255)

                                                 aTTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGt                                        >  1:1056445/1‑93 (MQ=255)

                                                 aTTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGt                                        >  1:4623351/1‑93 (MQ=255)

                                                  ttGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTa                                       <  2:2155747/93‑1 (MQ=255)

                                                   tGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTaa                                      <  1:1332827/93‑1 (MQ=255)

                                                    ggCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAAc                                     >  1:3496487/1‑93 (MQ=255)

                                                     gCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACg                                    >  1:3956560/1‑93 (MQ=255)

                                                       aaCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTc                                  <  2:1119166/93‑1 (MQ=255)

                                                        aCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCa                                 >  1:3065826/1‑93 (MQ=255)

                                                        aCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCa                                 >  2:2286251/1‑93 (MQ=255)

                                                        aCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCa                                 >  2:252786/1‑93 (MQ=255)

                                                          gTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCAcc                               >  2:4373019/1‑93 (MQ=255)

                                                           ttGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCg                              >  1:750971/1‑93 (MQ=255)

                                                           ttGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCg                              <  2:1816038/93‑1 (MQ=255)

                                                            tGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGa                             >  2:1465614/1‑93 (MQ=255)

                                                             gAGTGTGTTCAGGCGAACGCTGAGTCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGAc                            >  2:1612530/1‑93 (MQ=255)

                                                              aGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACa                           >  2:3025694/1‑93 (MQ=255)

                                                                tgtGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACaaa                         >  1:2114109/3‑93 (MQ=255)

                                                                gttGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACaaa                         <  2:1234940/91‑1 (MQ=255)

                                                                 ttGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAg                        <  2:259903/92‑1 (MQ=255)

                                                                  tGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGa                       <  1:2343797/93‑1 (MQ=255)

                                                                  tGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGa                       <  2:2045280/93‑1 (MQ=255)

                                                                  tGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGa                       >  1:2205899/1‑93 (MQ=255)

                                                                   gTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGAt                      >  1:4409698/1‑93 (MQ=255)

                                                                   gTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGAt                      <  2:4414354/93‑1 (MQ=255)

                                                                    ttCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATACCACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAtata                     >  2:4663161/1‑93 (MQ=255)

                                                                     tCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATaa                    <  2:3457877/93‑1 (MQ=255)

                                                                       aGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAg                  <  2:565619/93‑1 (MQ=255)

                                                                        ggCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGc                 >  2:2508992/1‑93 (MQ=255)

                                                                           gAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACt              >  2:2504708/1‑93 (MQ=255)

                                                                           gAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACt              >  2:614629/1‑93 (MQ=255)

                                                                           gAACGCTGACGCTCTACCGCCCCCGAATAACCCCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACt              <  2:4086957/93‑1 (MQ=255)

                                                                            aaCGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTg             >  2:2726918/1‑93 (MQ=255)

                                                                             aCGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGc            >  1:1171750/1‑93 (MQ=255)

                                                                             aCGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGc            <  1:3257535/93‑1 (MQ=255)

                                                                             aCGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGc            <  2:1999615/93‑1 (MQ=255)

                                                                              cGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCg           <  1:3191659/93‑1 (MQ=255)

                                                                              cGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCg           <  2:3755962/93‑1 (MQ=255)

                                                                                  gAGGCGCTCCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCAC‑GCGCTCaa      >  1:707521/1‑93 (MQ=255)

                                                                                  gAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCa       <  1:2380859/93‑1 (MQ=255)

                                                                                    ggCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAAt     <  2:4563788/93‑1 (MQ=255)

                                                                                    ggCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAAt     <  1:4027566/93‑1 (MQ=255)

                                                                                     gcgcTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATg    <  2:3600435/93‑1 (MQ=255)

                                                                                     gcgcTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATg    >  2:62641/1‑93 (MQ=255)

                                                                                      cgcTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATGt   >  2:4643615/1‑93 (MQ=255)

                                                                                       gcTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTCGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATGTa  <  2:2857935/93‑1 (MQ=255)

TGAATCCATGCCCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATGTA  >  NC_000913/4016732‑4016911

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