BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,433,176	T→C	R111R (CGT→CGC)	sun →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,433,176	0	T	C	95.7%	619.3 / 7.8	209	R111R (CGT→CGC)	sun	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (5/4); new base C (72/128); total (77/132)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

TATGGTTGGTTTGTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGTCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGATGC  >  NC_000913/3433086‑3433265

tATGGTTGGTTTGTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGccc                                                                                         <  1:1621477/93‑1 (MQ=255)

 aTGGTTGGTTTGTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCg                                                                                        <  1:2902053/93‑1 (MQ=255)

   ggttggttTGTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCa                                                                                      >  2:4355244/1‑93 (MQ=255)

       ggttTGTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACtt                                                                                  <  1:3879361/93‑1 (MQ=255)

          ttGTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTaaa                                                                               <  2:4293406/93‑1 (MQ=255)

           tGTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGTCCGCAACTTAAAg                                                                              >  2:1972059/1‑93 (MQ=255)

            gTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAgg                                                                             <  2:1924308/93‑1 (MQ=255)

            gTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAgg                                                                             <  2:2809386/93‑1 (MQ=255)

            gTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAgg                                                                             <  2:3265303/93‑1 (MQ=255)

            gTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGGTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAgg                                                                             <  1:744295/93‑1 (MQ=255)

               tCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGtt                                                                          <  2:4640602/93‑1 (MQ=255)

               tCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGtt                                                                          <  2:1341420/93‑1 (MQ=255)

                cAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTg                                                                         <  1:2320696/93‑1 (MQ=255)

                 accTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGa                                                                        <  2:1665380/91‑1 (MQ=255)

                 aaCTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGa                                                                        >  2:3931119/1‑93 (MQ=255)

                  aCTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGAt                                                                       <  2:2618682/93‑1 (MQ=255)

                  aCTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGAt                                                                       >  1:1383387/1‑93 (MQ=255)

                  aCTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGAt                                                                       >  2:2132775/1‑93 (MQ=255)

                    tGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATaa                                                                     <  1:4310009/93‑1 (MQ=255)

                     gCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATaaa                                                                    <  2:330295/93‑1 (MQ=255)

                      cTTTATACCCGCGTTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAAc                                                                   <  1:1091712/93‑1 (MQ=255)

                      cTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAAc                                                                   >  1:4446557/1‑93 (MQ=255)

                        ttATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGTCCGCAACTTAAAGGGTTGATAAACgg                                                                 <  1:4518589/93‑1 (MQ=255)

                        ttATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACgg                                                                 <  1:3528785/93‑1 (MQ=255)

                         tataCCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGgt                                                                <  1:4487566/93‑1 (MQ=255)

                          ataCCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGgtg                                                               >  1:3665939/1‑93 (MQ=255)

                          ataCCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGgtg                                                               >  2:2369386/1‑93 (MQ=255)

                           taCCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGgtgt                                                              <  1:2131198/93‑1 (MQ=255)

                             cccGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTAt                                                            <  1:1944529/93‑1 (MQ=255)

                             cccGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTAt                                                            >  1:89036/1‑93 (MQ=255)

                               cGCATTCCGCCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTa                                                          <  2:659728/93‑1 (MQ=255)

                                gCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTAc                                                         <  2:3849141/93‑1 (MQ=255)

                                gCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTAc                                                         <  2:4141353/93‑1 (MQ=255)

                                gCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTAc                                                         >  1:4055327/1‑93 (MQ=255)

                                   ttCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGcc                                                      <  1:2279903/93‑1 (MQ=255)

                                   ttCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGcc                                                      <  2:4249531/93‑1 (MQ=255)

                                      cACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGt                                                   >  1:4149928/1‑93 (MQ=255)

                                      cACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGt                                                   <  2:629842/93‑1 (MQ=255)

                                      cACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGt                                                   <  2:3013853/93‑1 (MQ=255)

                                       aCCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGtt                                                  <  2:3577194/93‑1 (MQ=255)

                                       aCCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGtt                                                  >  2:1656250/1‑93 (MQ=255)

                                         cTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGTCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTcc                                                <  1:3570527/93‑1 (MQ=255)

                                         cTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTcc                                                >  2:3100391/1‑93 (MQ=255)

                                          tCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACTGTGTATTACGCCAGTTCCa                                               <  2:3547957/93‑1 (MQ=255)

                                          tCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGCTCCa                                               >  2:4242734/1‑93 (MQ=255)

                                           cATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAg                                              <  1:269574/93‑1 (MQ=255)

                                           cATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAg                                              <  1:2369386/93‑1 (MQ=255)

                                            aTGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGc                                             >  1:3015136/1‑93 (MQ=255)

                                               ctgcGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTc                                          <  1:1332211/93‑1 (MQ=255)

                                               ctgcGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTc                                          >  2:2287638/1‑93 (MQ=255)

                                                tgcGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCa                                         <  2:4310112/93‑1 (MQ=255)

                                                tgcGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCa                                         <  1:1230972/93‑1 (MQ=255)

                                                tgcGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCa                                         >  1:3555521/1‑93 (MQ=255)

                                                    ctggctgAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCaa                                     >  1:3540831/1‑93 (MQ=255)

                                                      ggctgAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCaaga                                   <  1:1773235/93‑1 (MQ=255)

                                                      ggctgAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCaaga                                   >  2:3349014/1‑93 (MQ=255)

                                                       gctgAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCaagaa                                  <  1:2727104/93‑1 (MQ=255)

                                                        ctaaaaCGGGTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCaagaag                                 <  1:1884535/90‑1 (MQ=255)

                                                         tgAAACGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAagaaga                                <  2:3137827/93‑1 (MQ=255)

                                                            aaCGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGtt                             <  2:3284049/93‑1 (MQ=255)

                                                            aaCGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAATTCCAGCGTCAGCAAGAAGAGta                             >  2:2591645/1‑92 (MQ=255)

                                                             aCGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGtta                            >  2:2522705/1‑93 (MQ=255)

                                                             aCGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGtta                            <  1:4191517/93‑1 (MQ=255)

                                                             aCGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGtta                            >  2:1844441/1‑93 (MQ=255)

                                                              cGGTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGttat                           <  2:3649805/93‑1 (MQ=255)

                                                               ggTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGttatt                          <  2:2027331/93‑1 (MQ=255)

                                                               ggTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGttatt                          <  2:419031/93‑1 (MQ=255)

                                                                gTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGttatta                         <  1:3733989/93‑1 (MQ=255)

                                                                gTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGttatta                         <  1:983782/93‑1 (MQ=255)

                                                                gTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGttatta                         >  2:171917/1‑93 (MQ=255)

                                                                gTTGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGttatta                         <  1:36357/93‑1 (MQ=255)

                                                                 ttGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAg                        >  1:1670847/1‑93 (MQ=255)

                                                                  tGAAGGCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGc                       <  1:4557713/93‑1 (MQ=255)

                                                                      ggCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAg                   <  2:1524770/93‑1 (MQ=255)

                                                                      ggCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAg                   <  2:2909332/93‑1 (MQ=255)

                                                                      ggCGCTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAg                   >  1:867640/1‑93 (MQ=255)

                                                                       gcgcTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGt                  >  1:2756294/1‑93 (MQ=255)

                                                                       gcgcTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGt                  <  1:3055422/93‑1 (MQ=255)

                                                                        cgcTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGtt                 <  1:4098175/93‑1 (MQ=255)

                                                                         gcTATCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGttt                <  2:4202327/93‑1 (MQ=255)

                                                                            aTCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAAt             <  1:2553787/93‑1 (MQ=255)

                                                                            aTCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAAt             >  2:4568177/1‑93 (MQ=255)

                                                                             tCGCAATTAAGCGTCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAAtc            >  2:2759898/1‑92 (MQ=255)

                                                                             tCGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATg            <  1:1129959/93‑1 (MQ=255)

                                                                              cGCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGc           >  1:2596316/1‑93 (MQ=255)

                                                                               gCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGcc          <  2:3690150/93‑1 (MQ=255)

                                                                               gCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGcc          >  2:34730/1‑93 (MQ=255)

                                                                               gCAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGcc          <  1:4671060/93‑1 (MQ=255)

                                                                                cAATTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCa         <  1:3363447/93‑1 (MQ=255)

                                                                                 aaTTAAGCGCCCGTAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAg        >  1:3137106/1‑93 (MQ=255)

                                                                                 aaTTAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAg        <  2:3224817/93‑1 (MQ=255)

                                                                                   ttAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTg      >  1:1253771/1‑93 (MQ=255)

                                                                                    tAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGa     <  2:3792829/93‑1 (MQ=255)

                                                                                    tAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGa     <  2:3429817/93‑1 (MQ=255)

                                                                                    tAAGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGa     <  1:149915/93‑1 (MQ=255)

                                                                                     aaGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGAt    <  2:4014814/93‑1 (MQ=255)

                                                                                     aaGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGAt    <  1:2214041/93‑1 (MQ=255)

                                                                                     aaGCGCCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGAt    <  2:677366/93‑1 (MQ=255)

                                                                                      aGCGCCCGCAACTTAAAGGGTTGATAAGCGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGATg   <  1:1790114/93‑1 (MQ=255)

                                                                                       gcgcCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGATGc  >  2:2636493/1‑93 (MQ=255)

TATGGTTGGTTTGTATCAACTGCTTTATACCCGCATTCCACCTCATGCTGCGCTGGCTGAAACGGTTGAAGGCGCTATCGCAATTAAGCGTCCGCAACTTAAAGGGTTGATAAACGGTGTATTACGCCAGTTCCAGCGTCAGCAAGAAGAGTTATTAGCCGAGTTTAATGCCAGTGATGC  >  NC_000913/3433086‑3433265

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