BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,078,857	T→C	Y143C (TAC→TGC)	sbmC ←	SbmC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,078,857	0	T	C	98.5%	399.4 / ‑4.6	133	Y143C (TAC→TGC)	sbmC	SbmC
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (61/70); total (63/70)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 135 total aligned reads displayed.

AAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAG  >  NC_000913/2078769‑2078946

aaaaaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTAccc                                                                                       >  1:3856872/1‑93 (MQ=255)

aaaaaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACAGTTCGATATCCCAGCAccc                                                                                       >  1:1963511/1‑93 (MQ=255)

 aaaaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCa                                                                                      >  1:463109/1‑93 (MQ=255)

  aaaGAGCGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCAt                                                                                     >  2:1437724/1‑93 (MQ=255)

   aaGAGTGGTGATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATAGCCCAGCACCCATc                                                                                    <  2:1871405/93‑1 (MQ=255)

   aaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATc                                                                                    >  1:2471405/1‑93 (MQ=255)

    agagTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCt                                                                                   <  2:2793348/93‑1 (MQ=255)

    agagTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCt                                                                                   <  1:618669/93‑1 (MQ=255)

      agTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTc                                                                                 <  1:138392/93‑1 (MQ=255)

       gTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTcc                                                                                >  1:1560745/1‑93 (MQ=255)

        tGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCcg                                                                               >  1:2668286/1‑93 (MQ=255)

        tGGTCATCGCGTTAACACACCGCCCTGAGACGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCcg                                                                               >  1:2649117/1‑93 (MQ=255)

          gTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCcgcg                                                                             >  1:3537924/1‑93 (MQ=255)

          gTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCcgcg                                                                             >  2:1669112/1‑93 (MQ=255)

           tCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCgcgc                                                                            >  1:4664069/1‑93 (MQ=255)

            cATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGcc                                                                           >  2:291670/1‑93 (MQ=255)

             aTCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCa                                                                          >  1:811112/1‑93 (MQ=255)

               cgcgTTAACACGCCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCAtt                                                                        <  1:2824143/93‑1 (MQ=255)

               cgcgTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCAtt                                                                        <  1:897279/93‑1 (MQ=255)

                gcgTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTg                                                                       >  1:942452/1‑93 (MQ=255)

                        acacCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCGACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATa                                                               >  1:2510549/1‑93 (MQ=255)

                          acCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATaaa                                                             >  1:2391853/1‑93 (MQ=255)

                           ccGCCCTGAGATGAATTATTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAAc                                                            >  2:691146/1‑93 (MQ=255)

                           ccGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAAc                                                            >  1:689477/1‑93 (MQ=255)

                            cGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAAcc                                                           <  2:3750717/93‑1 (MQ=255)

                            cGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAAcc                                                           <  2:119367/93‑1 (MQ=255)

                            cGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAAcc                                                           >  1:103384/1‑93 (MQ=255)

                             gCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCt                                                          <  1:3208103/93‑1 (MQ=255)

                              cccTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCCCCCATCTTCCGCGCCATTGTTCAAATAAACCTc                                                         <  2:506163/93‑1 (MQ=255)

                              cccTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTc                                                         >  1:4116240/1‑93 (MQ=255)

                              cccTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTc                                                         <  2:2098481/93‑1 (MQ=255)

                               ccTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCg                                                        >  2:4048226/1‑93 (MQ=255)

                               ccTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCg                                                        >  2:267386/1‑93 (MQ=255)

                                cTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGa                                                       >  2:246890/1‑93 (MQ=255)

                                cTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGa                                                       <  1:3349952/93‑1 (MQ=255)

                                cTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGa                                                       <  1:4185472/93‑1 (MQ=255)

                                cTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGa                                                       >  1:2678873/1‑93 (MQ=255)

                                 tGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGaa                                                      >  2:718212/1‑93 (MQ=255)

                                 tGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGaa                                                      <  2:1998884/93‑1 (MQ=255)

                                   agaTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGc                                                    >  1:4406223/1‑93 (MQ=255)

                                     aTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAg                                                  <  1:913774/93‑1 (MQ=255)

                                     aTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAg                                                  <  2:912332/93‑1 (MQ=255)

                                     aTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAg                                                  <  1:2259111/93‑1 (MQ=255)

                                     aTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAg                                                  <  1:188030/93‑1 (MQ=255)

                                     aTGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAg                                                  <  1:2505544/93‑1 (MQ=255)

                                      tGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAgg                                                 >  1:1832402/1‑93 (MQ=255)

                                      tGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAgg                                                 <  2:3370016/93‑1 (MQ=255)

                                      tGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAgg                                                 >  2:3251705/1‑93 (MQ=255)

                                        aaTTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGc                                               >  2:3751104/1‑93 (MQ=255)

                                        aaTTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGc                                               >  2:1585710/1‑93 (MQ=255)

                                         aTTAGTGATGTTTTGGCTTCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCt                                              <  2:1201829/93‑1 (MQ=255)

                                           tAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCttt                                            <  1:2495331/93‑1 (MQ=255)

                                            aGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTg                                           <  1:2169510/93‑1 (MQ=255)

                                            aGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTg                                           <  2:4538569/93‑1 (MQ=255)

                                             gTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTgg                                          >  2:1526017/1‑93 (MQ=255)

                                              tGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGt                                         <  2:2223685/93‑1 (MQ=255)

                                               gATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTa                                        >  2:153816/1‑93 (MQ=255)

                                                aTGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTaa                                       <  1:4246131/93‑1 (MQ=255)

                                                aTGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTaa                                       >  1:1527184/1‑93 (MQ=255)

                                                aTGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTaa                                       <  1:1476064/93‑1 (MQ=255)

                                                 tGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAAc                                      >  1:246453/1‑93 (MQ=255)

                                                 tGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAAc                                      >  2:3108288/1‑93 (MQ=255)

                                                    tttGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACAtt                                   <  1:3884658/93‑1 (MQ=255)

                                                     ttGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACAttt                                  >  2:279188/1‑93 (MQ=255)

                                                     ttGGCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACAttt                                  <  1:3437645/93‑1 (MQ=255)

                                                      tGGCTGCACCGCAACATACATTGCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTc                                 <  2:1384395/93‑1 (MQ=255)

                                                       ggCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAACCCTCGAAGCAGGGCTTTGGTAACATTTCa                                <  2:967906/93‑1 (MQ=255)

                                                       ggCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCa                                >  1:4517505/1‑93 (MQ=255)

                                                       ggCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCa                                >  2:34538/1‑93 (MQ=255)

                                                        gCTGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCAt                               <  2:3528674/93‑1 (MQ=255)

                                                          tGCACCGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATaa                             <  2:4673559/93‑1 (MQ=255)

                                                             cccGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCa                          <  2:898078/92‑1 (MQ=255)

                                                              ccGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCAc                         <  2:1728299/93‑1 (MQ=255)

                                                               cGCAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACt                        <  2:2575660/93‑1 (MQ=255)

                                                                 cAACATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGt                      <  2:613449/93‑1 (MQ=255)

                                                                  aaCATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTc                     >  2:824939/1‑93 (MQ=255)

                                                                  aaCATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTc                     <  2:3143750/93‑1 (MQ=255)

                                                                  aaCATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTc                     <  2:989227/93‑1 (MQ=255)

                                                                  aaCATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTc                     >  1:1034912/1‑93 (MQ=255)

                                                                  aaCATACATTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTc                     >  1:371705/1‑93 (MQ=255)

                                                                   acatacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTcc                    <  1:115820/93‑1 (MQ=255)

                                                                   acatacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTcc                    <  1:2724045/93‑1 (MQ=255)

                                                                   acatacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTcc                    >  1:3735132/1‑93 (MQ=255)

                                                                     atacatTTCGATATCCCAGCACCCATCTTCAGCGCCATTGATCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTg                  >  2:1525463/1‑93 (MQ=255)

                                                                      tacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGc                 >  1:3748772/1‑93 (MQ=255)

                                                                      tacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGc                 >  1:2822256/1‑93 (MQ=255)

                                                                      tacatTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGc                 >  2:2293936/1‑93 (MQ=255)

                                                                        catTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCaa               <  1:4622211/93‑1 (MQ=255)

                                                                        catTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCaa               <  1:3863354/93‑1 (MQ=255)

                                                                        catTTCGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCaa               <  1:4650299/93‑1 (MQ=255)

                                                                             cGATATCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGAc          <  1:2022025/93‑1 (MQ=255)

                                                                              gATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACt         >  2:193066/1‑93 (MQ=255)

                                                                               atatCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTa        >  1:3124038/1‑93 (MQ=255)

                                                                                tatCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTAt       <  1:4580812/93‑1 (MQ=255)

                                                                                tatCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTAt       <  2:2509155/93‑1 (MQ=255)

                                                                                 atCCCAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTAtt      <  2:1188806/93‑1 (MQ=255)

                                                                                   cccAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTaa    <  2:2893184/93‑1 (MQ=255)

                                                                                    ccAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTaaa   <  1:3080168/93‑1 (MQ=255)

                                                                                    ccAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTaaa   <  1:2229750/93‑1 (MQ=255)

                                                                                     cAGCACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAg  <  1:477929/93‑1 (MQ=255)

AAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAG  >  NC_000913/2078769‑2078946

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