BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,966	T→C	S241G (AGC→GGC)	lplA ←	lipoate‑protein ligase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,966	0	T	C	82.8%	581.0 / 120.7	273	S241G (AGC→GGC)	lplA	lipoate‑protein ligase A
Reads supporting (aligned to +/- strand): ref base T (20/27); new base C (111/115); total (131/142)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

Only 100 of 275 total aligned reads displayed.

ACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGG  >  NC_000913/4620875‑4621058

aCGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTa                                                                                             >  2:257885/1‑93 (MQ=255)

 cGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCAc                                                                                            <  2:6487041/93‑1 (MQ=255)

  gTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCAct                                                                                           <  2:4856239/93‑1 (MQ=255)

  gTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCAct                                                                                           >  1:4194335/1‑93 (MQ=255)

  gTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCAct                                                                                           >  2:5139297/1‑93 (MQ=255)

   tCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCActc                                                                                          <  2:314068/93‑1 (MQ=255)

   tCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCActc                                                                                          <  2:2291049/93‑1 (MQ=255)

     gAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTg                                                                                        >  2:907434/1‑93 (MQ=255)

      aaaTGCAGTTCCACGCCGCTCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTgg                                                                                       <  2:3693441/93‑1 (MQ=255)

      aaaTGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTgg                                                                                       <  1:7520693/93‑1 (MQ=255)

       aaTGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGc                                                                                      >  2:4602557/1‑93 (MQ=255)

        aTGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCg                                                                                     >  1:237006/1‑93 (MQ=255)

         tGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCGTGACCGAAGTTCCATTCCCAGCCACTCTGGCgg                                                                                    <  2:3813224/93‑1 (MQ=255)

         tGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCgg                                                                                    >  1:6344440/1‑93 (MQ=255)

         tGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCgg                                                                                    <  2:4490032/93‑1 (MQ=255)

          gCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCggg                                                                                   >  1:1642374/1‑93 (MQ=255)

             gTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCaa                                                                                >  2:6065501/1‑93 (MQ=255)

              ttCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCaaa                                                                               <  2:4740501/93‑1 (MQ=255)

              ttCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCaaa                                                                               <  1:4909464/93‑1 (MQ=255)

               tCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAg                                                                              >  1:5513704/1‑93 (MQ=255)

                ccACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAgg                                                                             <  1:6817943/93‑1 (MQ=255)

                  aCGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGtt                                                                           >  1:5330387/1‑93 (MQ=255)

                  aCGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGtt                                                                           <  2:2463087/93‑1 (MQ=255)

                   cgccgcCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGttt                                                                          <  2:4248406/93‑1 (MQ=255)

                   cgccgcCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGttt                                                                          <  1:5567758/93‑1 (MQ=255)

                      cgccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCgg                                                                       <  1:3692792/93‑1 (MQ=255)

                      cgccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCgg                                                                       >  2:3127384/1‑93 (MQ=255)

                       gccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGc                                                                      >  1:1955323/1‑93 (MQ=255)

                          ccAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGaa                                                                   <  2:2687182/93‑1 (MQ=255)

                          ccAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGaa                                                                   <  2:3656528/93‑1 (MQ=255)

                           cAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAg                                                                  >  1:7085698/1‑93 (MQ=255)

                             ggTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGtt                                                                >  2:391197/1‑93 (MQ=255)

                              gTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGttt                                                               >  1:6470211/1‑93 (MQ=255)

                               tAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTg                                                              <  2:6909144/93‑1 (MQ=255)

                                aaaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTgg                                                             <  1:726133/93‑1 (MQ=255)

                                 aaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGc                                                            <  1:7020758/93‑1 (MQ=255)

                                 aaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGc                                                            <  1:2827978/93‑1 (MQ=255)

                                   gCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCaa                                                          <  2:4632201/93‑1 (MQ=255)

                                   gCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCaa                                                          >  2:2268863/1‑93 (MQ=255)

                                   gCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCaa                                                          >  2:1774995/1‑93 (MQ=255)

                                   gCGTTCACCCAGCAGATGCGAGAAGGCCGGAGCCTGACCGAAGTTCCATCCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCaa                                                          <  1:2693925/93‑1 (MQ=255)

                                    cGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAg                                                         <  2:4135439/93‑1 (MQ=255)

                                     gTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGt                                                        <  2:1395605/93‑1 (MQ=255)

                                      ttCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTc                                                       <  1:1753174/93‑1 (MQ=255)

                                       tCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCt                                                      <  2:792419/93‑1 (MQ=255)

                                        cATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTg                                                     <  2:1295135/93‑1 (MQ=255)

                                        cATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTg                                                     >  2:1924895/1‑93 (MQ=255)

                                        aaTCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTg                                                     <  1:5129286/92‑1 (MQ=255)

                                         aTCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTgg                                                    <  1:963276/93‑1 (MQ=255)

                                         aTCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCAGCGAAGTTTGGCAAGTCTgg                                                    >  2:438016/1‑93 (MQ=255)

                                          tCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGc                                                   <  2:1562983/93‑1 (MQ=255)

                                          tCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGc                                                   >  2:7024019/1‑93 (MQ=255)

                                            cagcagATGCGAGAATGCCGGAGCCTGACCGCAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTTGGGGTAGCTTGGTTAGTCTGGCGt                                                 >  1:5238964/1‑93 (MQ=255)

                                            cagcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGt                                                 >  1:4669134/1‑93 (MQ=255)

                                            cagcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGt                                                 >  2:3953888/1‑93 (MQ=255)

                                             agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGtt                                                >  2:5726869/1‑93 (MQ=255)

                                              gcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGttt                                               >  1:3533291/1‑93 (MQ=255)

                                               cagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGtttt                                              <  2:5573797/93‑1 (MQ=255)

                                                agATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTg                                             <  2:560696/93‑1 (MQ=255)

                                                agATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTg                                             >  2:3068697/1‑93 (MQ=255)

                                                 gATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGt                                            >  2:2797710/1‑93 (MQ=255)

                                                   tGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTc                                          >  2:6612716/1‑93 (MQ=255)

                                                   tGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTc                                          >  2:2093449/1‑93 (MQ=255)

                                                   tGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTc                                          <  2:1522779/93‑1 (MQ=255)

                                                     cGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCgg                                        <  2:6641570/93‑1 (MQ=255)

                                                     cGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCgg                                        <  2:4568109/93‑1 (MQ=255)

                                                      gagaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCggg                                       <  1:3967078/93‑1 (MQ=255)

                                                       agaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCgggg                                      <  1:3040961/93‑1 (MQ=255)

                                                       agaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCgggg                                      <  1:18202/93‑1 (MQ=255)

                                                        gaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGga                                     <  1:2262761/93‑1 (MQ=255)

                                                        gaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGga                                     >  2:6662527/1‑93 (MQ=255)

                                                          aTGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGgaga                                   >  2:6161300/1‑93 (MQ=255)

                                                          aTGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGgaga                                   <  2:7190375/93‑1 (MQ=255)

                                                          aTGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGgaga                                   <  2:4915279/93‑1 (MQ=255)

                                                           tGCCGGAGCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAgat                                  >  1:5828939/1‑93 (MQ=255)

                                                               ggAGCCTGACCGAATTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGAtt                              <  1:4537392/93‑1 (MQ=255)

                                                                 aGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTc                            <  2:2149215/93‑1 (MQ=255)

                                                                 aGCCTGACCGAAGTTCCATCCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTc                            <  1:4046597/93‑1 (MQ=255)

                                                                  gCCTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTcc                           >  1:1242568/1‑93 (MQ=255)

                                                                   ccTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCg                          <  2:2968126/93‑1 (MQ=255)

                                                                    cTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGc                         >  2:1646395/1‑93 (MQ=255)

                                                                    cTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGc                         >  1:140772/1‑93 (MQ=255)

                                                                    cTGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGc                         >  1:3299555/1‑93 (MQ=255)

                                                                     tGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCt                        >  1:4608174/1‑93 (MQ=255)

                                                                     tGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCt                        >  2:4520373/1‑93 (MQ=255)

                                                                     tGACCGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCt                        >  2:6332869/1‑93 (MQ=255)

                                                                        ccGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTcc                     >  2:5333715/1‑93 (MQ=255)

                                                                         cGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCa                    >  1:6170532/1‑93 (MQ=255)

                                                                         cGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCa                    <  2:3426978/93‑1 (MQ=255)

                                                                         cGAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCa                    <  2:3986645/93‑1 (MQ=255)

                                                                          gAAGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCAc                   <  2:5024281/93‑1 (MQ=255)

                                                                            aGTTCCATTCCCAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCAcgc                 <  1:1155444/93‑1 (MQ=255)

                                                                              ttCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCCAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc               >  1:6584116/1‑93 (MQ=255)

                                                                                      ccAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATa       >  1:892895/1‑93 (MQ=255)

                                                                                       cAGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATaa      <  1:2059532/93‑1 (MQ=255)

                                                                                        aGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAt     <  2:5068041/93‑1 (MQ=255)

                                                                                        aGCCACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAt     >  1:838940/1‑93 (MQ=255)

                                                                                          cTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATgg   >  1:3377957/1‑93 (MQ=255)

                                                                                           tACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATggg  >  1:2874829/1‑93 (MQ=255)

                                                                                           tACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATggg  >  2:3262956/1‑93 (MQ=255)

ACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGG  >  NC_000913/4620875‑4621058

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