BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,056,882	T→C	Y99C (TAC→TGC)	torR ←	TorR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,056,882	0	T	C	98.6%	455.8 / ‑3.0	146	Y99C (TAC→TGC)	torR	TorR
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (68/76); total (70/76)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 152 total aligned reads displayed.

TGCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCT  >  NC_000913/1056792‑1056971

  cGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATc                                                                                       <  1:1834582/93‑1 (MQ=255)

     ttGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTc                                                                                    <  2:4057419/93‑1 (MQ=255)

     gtGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTc                                                                                    <  2:4057418/92‑1 (MQ=255)

      tGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCt                                                                                   <  1:686623/93‑1 (MQ=255)

       gAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTg                                                                                  >  1:3624086/1‑93 (MQ=255)

        aGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTgc                                                                                 <  1:4671124/93‑1 (MQ=255)

         gCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCTTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTgcg                                                                                <  1:3613324/93‑1 (MQ=255)

         gCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTgcg                                                                                >  1:533602/1‑93 (MQ=255)

          cTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTgcgc                                                                               >  1:1441319/1‑93 (MQ=255)

          cTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTgcgc                                                                               <  1:1535418/93‑1 (MQ=255)

           ttGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGcc                                                                              <  1:1197099/93‑1 (MQ=255)

           ttGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGcc                                                                              >  2:3932361/1‑93 (MQ=255)

           ttGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGCTTGGTGACGCAATCGTCTGCGcc                                                                              >  2:280858/1‑93 (MQ=255)

             gTCGCGCGAGGTCGATTCGCCAGAGCATATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCa                                                                            <  1:942698/93‑1 (MQ=255)

             gTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCa                                                                            <  2:668165/93‑1 (MQ=255)

              tCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCAt                                                                           <  1:2772268/93‑1 (MQ=255)

                gcgcgAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCAttt                                                                         >  1:807847/1‑93 (MQ=255)

                 cgcgAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTc                                                                        >  1:3934783/1‑93 (MQ=255)

                  gcgAGTTCGATTCGCCAGAGCAGTTTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTcc                                                                       <  2:3218311/93‑1 (MQ=255)

                  gcgAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTcc                                                                       >  1:3114062/1‑93 (MQ=255)

                  gcgAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGGTTGGTGACGTAATCGTCTGCGCCCATTTca                                                                       >  1:2211624/1‑92 (MQ=255)

                     gggTCGATTCGACAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGc                                                                    <  1:2408893/92‑1 (MQ=255)

                      ggTCGATTCGACAGAGCAGATTTTTCGCCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCTTCTGCGCCCATTTCCAGcc                                                                   <  2:4190586/93‑1 (MQ=255)

                       gTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGccc                                                                  <  2:128000/93‑1 (MQ=255)

                         cGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaa                                                                <  1:3982440/93‑1 (MQ=255)

                          gATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaac                                                               <  2:1614592/93‑1 (MQ=255)

                           aTTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaaca                                                              >  1:719636/1‑93 (MQ=255)

                           aTTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaaca                                                              >  1:1111501/1‑93 (MQ=255)

                           aTTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaaca                                                              >  2:2091011/1‑93 (MQ=255)

                            ttCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaacaa                                                             >  2:3470583/1‑93 (MQ=255)

                             tCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAAt                                                            >  2:1902856/1‑93 (MQ=255)

                              cGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATa                                                           <  1:3847143/93‑1 (MQ=255)

                              cGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATa                                                           <  1:4057332/93‑1 (MQ=255)

                                 cAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACgg                                                        <  2:121812/93‑1 (MQ=255)

                                  agagCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACACTACGGt                                                       >  1:3564945/1‑93 (MQ=255)

                                  agagCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGCGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGt                                                       >  1:4203345/1‑93 (MQ=255)

                                  agagCAGATTTTTCACCCGTACTACCAGTTCGCCCAGTTCCACCGGTTTGGTGACGAAATCGTCTGCGCCCATTTCCAGCCCACAAATACggc                                                       >  2:3438293/1‑92 (MQ=255)

                                   gagCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTc                                                      >  2:2501770/1‑93 (MQ=255)

                                   gagCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTc                                                      >  2:2977940/1‑93 (MQ=255)

                                   gagCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTc                                                      >  2:3066849/1‑93 (MQ=255)

                                   gagCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTc                                                      <  1:783566/93‑1 (MQ=255)

                                   gagCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTc                                                      >  1:2205638/1‑93 (MQ=255)

                                     gCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCaa                                                    >  2:2083152/1‑93 (MQ=255)

                                     gCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCaa                                                    <  1:3854130/93‑1 (MQ=255)

                                     gCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCaa                                                    >  1:1985812/1‑93 (MQ=255)

                                       aGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATc                                                  <  1:2087802/93‑1 (MQ=255)

                                        gATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATcc                                                 <  1:3325393/93‑1 (MQ=255)

                                         aTTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCg                                                <  2:2717528/93‑1 (MQ=255)

                                          tttttCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGa                                               <  1:346077/93‑1 (MQ=255)

                                            tttCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATc                                             <  2:1116383/93‑1 (MQ=255)

                                               cACCCGCCCTCCCAGTTCGCGCAATTCCAGCGGTTTGGTTACTCAATCGTCTGCGCCCATTTCCAGCCCAACAATTCGGTCAATCCGATCGCt                                          <  1:941518/93‑1 (MQ=255)

                                                aCCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTg                                         >  2:4485519/1‑93 (MQ=255)

                                                 cccGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGc                                        <  1:2915777/93‑1 (MQ=255)

                                                  ccGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCg                                       <  2:4181406/93‑1 (MQ=255)

                                                  ccGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCg                                       >  2:760020/1‑93 (MQ=255)

                                                   cGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCgc                                      >  1:1233553/1‑92 (MQ=255)

                                                   cGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGt                                      <  2:2694728/93‑1 (MQ=255)

                                                    gTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTc                                     >  1:1575777/1‑93 (MQ=255)

                                                     tactacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTcc                                    >  2:1906215/1‑93 (MQ=255)

                                                      actacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCg                                   <  1:1795797/93‑1 (MQ=255)

                                                      actacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCg                                   <  1:679150/93‑1 (MQ=255)

                                                       ctacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACTGTCAATCCGATCGCTGCGTCCgg                                  >  1:506136/1‑93 (MQ=255)

                                                          ccAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTaa                               <  2:4575801/93‑1 (MQ=255)

                                                             gTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCa                            <  2:2540543/93‑1 (MQ=255)

                                                              ttCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAg                           <  2:870708/93‑1 (MQ=255)

                                                                cgcgCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGaa                         <  1:2749079/93‑1 (MQ=255)

                                                                cgcgCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGaa                         <  1:1873194/93‑1 (MQ=255)

                                                                cgcgCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGaa                         <  1:1012305/93‑1 (MQ=255)

                                                                cgcgCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGaa                         <  2:2154955/93‑1 (MQ=255)

                                                                cccgcAGTGGCAGCTGTTTAGCGAGGTAATCGTATGCGCCCTTTTCCAGCAGATTAATACGGTCAATCCGCTCGCTGAGCCCGGTAACCAGaa                         <  2:3919820/91‑1 (MQ=255)

                                                                 gcgcAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGaat                        >  2:158133/1‑93 (MQ=255)

                                                                  cgcAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGaata                       <  1:1821212/93‑1 (MQ=255)

                                                                  cgcAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGaata                       >  2:2714127/1‑93 (MQ=255)

                                                                   gcAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGCCAATCCGATCGCTGCGTCCGGTAACCAGaataa                      >  2:3405069/1‑93 (MQ=255)

                                                                      gTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATcc                   >  1:2699716/1‑93 (MQ=255)

                                                                      gTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATcc                   <  2:2669441/93‑1 (MQ=255)

                                                                      gTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGCCCGGGAACCAGAATAATcc                   >  1:3383255/1‑93 (MQ=255)

                                                                       ttCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATccc                  >  2:1997955/1‑93 (MQ=255)

                                                                       ttCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATccc                  >  2:3584607/1‑93 (MQ=255)

                                                                         ccAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCa                >  2:4519463/1‑93 (MQ=255)

                                                                         ccAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCa                <  1:3111730/93‑1 (MQ=255)

                                                                          cAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCAc               <  2:1574819/93‑1 (MQ=255)

                                                                           aGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCAcc              >  1:3699832/1‑93 (MQ=255)

                                                                           aGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCAcc              >  1:3699831/1‑93 (MQ=255)

                                                                           aGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCAcc              <  2:1370165/93‑1 (MQ=255)

                                                                            gCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCg             <  2:3760434/93‑1 (MQ=255)

                                                                              ggTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGtt           <  1:1892772/93‑1 (MQ=255)

                                                                              ggTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGtt           >  2:3740673/1‑93 (MQ=255)

                                                                               gTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTg          <  2:460839/93‑1 (MQ=255)

                                                                               gTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTg          >  1:2698802/1‑93 (MQ=255)

                                                                                tttGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGa         <  1:229347/93‑1 (MQ=255)

                                                                                  tGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGc       <  1:4661011/93‑1 (MQ=255)

                                                                                   ggTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCg      >  1:1325044/1‑93 (MQ=255)

                                                                                   ggTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCg      <  2:1636278/93‑1 (MQ=255)

                                                                                     tGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGtt    <  2:2622837/93‑1 (MQ=255)

                                                                                     tGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGtt    <  2:216728/93‑1 (MQ=255)

                                                                                     ggACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGtt    <  2:2041231/92‑1 (MQ=255)

                                                                                      gACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTtc   <  1:1746113/93‑1 (MQ=255)

                                                                                      gACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTtc   <  1:4030733/93‑1 (MQ=255)

                                                                                       aCGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTtct  >  1:4267065/1‑93 (MQ=255)

TGCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCT  >  NC_000913/1056792‑1056971

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