BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,001,717	+C	coding (86/903 nt)	rarD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,001,717	1	.	C	99.7%	1399.6 / ‑4.0	376	coding (86/903 nt)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (139/236); total (140/236)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

Only 100 of 465 total aligned reads displayed.

AAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑G‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTT  >  NC_000913/4001622‑4001806

aaaTGCTCATCAGGACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑g                                                                                            >  1:10972613/1‑96 (MQ=255)

aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑g                                                                                            >  2:6218484/1‑96 (MQ=255)

aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑g                                                                                            <  2:10257064/96‑1 (MQ=255)

aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAAATTGAAGTA‑C‑g                                                                                            >  1:947000/1‑96 (MQ=255)

 aaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑gc                                                                                           >  1:281605/1‑95 (MQ=255)

  aTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑gct                                                                                          <  1:10665452/96‑3 (MQ=255)

  aTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑gct                                                                                          <  1:3159311/96‑3 (MQ=255)

    gCTCATCAGCACCACCATAAAGAAAAACGCCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑gctgg                                                                                        <  1:4994056/96‑5 (MQ=255)

    gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑gctgg                                                                                        >  1:7412902/1‑92 (MQ=255)

      tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGgcg                                                                                      <  2:10876417/96‑1 (MQ=255)

        atcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTATTAAATCAACTTGAAGTA‑T‑GCTGGCGCt                                                                                    >  1:7894482/1‑96 (MQ=255)

        atcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGAA‑C‑GCTGGCGCt                                                                                    >  2:6407150/1‑96 (MQ=255)

         tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCta                                                                                   >  1:13178723/1‑96 (MQ=255)

         tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCta                                                                                   <  2:1373469/96‑1 (MQ=255)

          caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCtat                                                                                  <  2:172388/96‑1 (MQ=255)

          caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATAGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCtat                                                                                  <  1:12418359/96‑1 (MQ=255)

              agcaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATAccc                                                                              >  1:10292495/3‑96 (MQ=255)

              accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATAccc                                                                              >  1:2276280/1‑96 (MQ=255)

                caccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTGGTAAATCAACTTGAAGTATG‑G‑TGGCGCTATACCCCa                                                                            >  1:3508621/1‑96 (MQ=255)

                 accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaa                                                                           <  2:4843197/96‑1 (MQ=255)

                 accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaa                                                                           <  2:530833/96‑1 (MQ=255)

                 accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaa                                                                           <  1:4810115/96‑1 (MQ=255)

                 accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaa                                                                           >  2:363339/1‑96 (MQ=255)

                  ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaa                                                                          <  1:3940548/96‑1 (MQ=255)

                  caaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATAGGCGGGCACGTAGTAAATCAACTTGAAGGA‑C‑GCTGGCGCTGTACCCCcaa                                                                          >  1:11820982/3‑96 (MQ=255)

                   caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaat                                                                         <  1:9657365/96‑1 (MQ=255)

                    aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAAACAACTTGAAGTA‑C‑GCTGGGGCTATACCCCaaata                                                                        >  2:8585912/1‑96 (MQ=255)

                      aaaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCGAGTAAATCAACTTGAACTA‑C‑GCTGGCGCTGTACCCCaaataaa                                                                      >  2:1243215/1‑96 (MQ=255)

                        aGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaataaaat                                                                    >  1:6344057/1‑96 (MQ=255)

                         gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaataaaata                                                                   <  1:8555560/96‑1 (MQ=255)

                         gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAATCAACTTTAAGGA‑CAC‑CTGGGCTATACCCCaaataaaata                                                                   >  2:5615595/1‑96 (MQ=255)

                         gAAAAACGACCAGATCACGCGATGCGTCACGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaataaaata                                                                   >  1:11234740/1‑96 (MQ=255)

                           aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAg                                                                 <  2:867995/96‑1 (MQ=255)

                            aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACATGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGc                                                                >  2:12240706/1‑96 (MQ=255)

                             aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCg                                                               <  2:13987050/96‑1 (MQ=255)

                             aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCg                                                               <  2:13404986/96‑1 (MQ=255)

                             aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCg                                                               <  1:5756871/96‑1 (MQ=255)

                             aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCg                                                               <  1:8229912/96‑1 (MQ=255)

                             aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCg                                                               >  1:4488506/1‑96 (MQ=255)

                              aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGTAGTAAATCAAATTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCgg                                                              >  2:13081703/1‑96 (MQ=255)

                               cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGc                                                             <  2:7412902/96‑1 (MQ=255)

                                 aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                           <  2:7744158/96‑1 (MQ=255)

                                 aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAAAAAACGGCaa                                                           >  1:2410124/1‑96 (MQ=255)

                                  ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                                          <  2:12543002/96‑1 (MQ=255)

                                   cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                         <  2:11818776/96‑1 (MQ=255)

                                   cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                         <  1:714827/96‑1 (MQ=255)

                                   cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                         <  1:1655604/96‑1 (MQ=255)

                                   cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                         <  1:6965101/96‑1 (MQ=255)

                                    aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGAGGAAAACAACTTGAAGTA‑C‑GCTTGGGCTATAACCCAAATAAAATAAGGGGCAagag                                                        >  2:7288059/1‑96 (MQ=255)

                                     gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                                       <  2:13313530/96‑1 (MQ=255)

                                     gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                                       <  2:10699696/96‑1 (MQ=255)

                                       tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                                     <  2:9158920/96‑1 (MQ=255)

                                       tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                                     >  1:10488306/1‑96 (MQ=255)

                                       tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTCGAAGTA‑C‑GCTGGCGCTCTACCCCAAATAAAATAAACGGCAATAGCGa                                                     >  2:11679925/1‑96 (MQ=255)

                                          cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGAGAGTa                                                  >  1:8954328/1‑96 (MQ=255)

                                           gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                                 <  2:8086912/96‑1 (MQ=255)

                                            cgATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAATCAACTTGAAGTA‑A‑GCTGGCGCTATACCCCAAAAAAAATAAGCGGCAAGAGCGAGTAAt                                                >  2:4336564/1‑96 (MQ=255)

                                                 cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑A‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGcc                                           >  2:10384385/1‑96 (MQ=255)

                                                 cGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGcc                                           >  1:1992736/1‑96 (MQ=255)

                                                  gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                          >  2:12892973/1‑96 (MQ=255)

                                                  gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                          >  2:4028625/1‑96 (MQ=255)

                                                  gTCAGGATTTCATCGGCGGGCACGTAGGAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                          >  2:9434911/1‑96 (MQ=255)

                                                    cAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTg                                        <  1:5133282/96‑1 (MQ=255)

                                                    cAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTg                                        <  2:948630/96‑1 (MQ=255)

                                                       gATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                     <  2:1961275/96‑1 (MQ=255)

                                                       gATTTCATCGGCGGGCACGTAGGAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                     >  2:3048934/1‑96 (MQ=255)

                                                        aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGGGTAATACGCCCTGCcgc                                    >  1:5909077/1‑96 (MQ=255)

                                                          ttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                                  <  1:14108961/96‑1 (MQ=255)

                                                          ttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                                  <  2:6921400/96‑1 (MQ=255)

                                                          ttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCGGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                                  >  2:2453333/1‑96 (MQ=255)

                                                           tCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGGAATACGCCCTGCCGCgtt                                 >  2:5051953/1‑96 (MQ=255)

                                                            cATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt                                >  1:9971189/1‑96 (MQ=255)

                                                               cggcggGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt                             >  1:11657581/1‑96 (MQ=255)

                                                                ggcggGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt                            <  2:2149943/96‑1 (MQ=255)

                                                                ggcggGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtat                            <  1:14269021/96‑3 (MQ=255)

                                                                 gcggGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt                           <  1:5337296/96‑1 (MQ=255)

                                                                 gcggGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt                           <  1:12892973/96‑1 (MQ=255)

                                                                  cggGCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTg                          >  2:9597393/1‑96 (MQ=255)

                                                                    ggCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                        >  1:9774981/1‑96 (MQ=255)

                                                                     gCACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                       >  2:14215837/1‑96 (MQ=255)

                                                                      cACGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTGGcatc                      >  1:4106128/1‑96 (MQ=255)

                                                                        cGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCGCTGCCGCGTTTGTTTTGcatcca                    >  1:10120083/1‑96 (MQ=255)

                                                                        cGTAGTAAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                    <  2:2865231/96‑1 (MQ=255)

                                                                         gtagtaAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                   <  1:1208278/96‑1 (MQ=255)

                                                                          tagtaAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATATGCCCTGCCGCGTTTGTTTTGcatccatc                  <  2:11438894/96‑1 (MQ=255)

                                                                          tagtaAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc                  >  1:11513039/1‑96 (MQ=255)

                                                                           agtaAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg                 <  2:5813399/96‑1 (MQ=255)

                                                                           agtaAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg                 >  2:1877270/1‑96 (MQ=255)

                                                                           agtaAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg                 <  2:7949950/96‑1 (MQ=255)

                                                                            gtaAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGTGGCAAGAGCGAGTAATACGCCCTGCCGCTTTTGTTTTGCATCCATCgg                <  2:12765324/96‑1 (MQ=255)

                                                                            gtaAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg                >  1:913647/1‑96 (MQ=255)

                                                                            gtaAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg                >  2:10179024/1‑96 (MQ=255)

                                                                             taAATCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCggg               <  2:758711/96‑1 (MQ=255)

                                                                              aaaTCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa              <  2:14409087/96‑1 (MQ=255)

                                                                                aTCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGaag            <  1:13761735/96‑2 (MQ=255)

                                                                                aTCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt            <  2:4583666/96‑1 (MQ=255)

                                                                                     cTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCa       >  1:5468600/1‑96 (MQ=255)

                                                                                      ttGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAt      <  2:327475/96‑1 (MQ=255)

                                                                                         aaGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtttt   >  1:173180/1‑96 (MQ=255)

                                                                                          aGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAttttt  <  1:13562366/96‑1 (MQ=255)

AAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTA‑C‑G‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTT  >  NC_000913/4001622‑4001806

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