BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,373,758	T→C	S159S (TCA→TCG)	yjeJ ←	protein YjeJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,373,758	0	T	C	99.4%	1383.5 / NA	342	S159S (TCA→TCG)	yjeJ	protein YjeJ
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (159/181); total (159/183)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

AACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGTGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTGACGACTA  >  NC_000913/4373660‑4373855

 aCAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAg                                                                                                >  2:7464922/1‑101 (MQ=255)

   aaGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAgg                                                                                               >  1:7464925/1‑100 (MQ=255)

   aaGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGt                                                                                              <  2:7464353/101‑1 (MQ=255)

   aaGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGt                                                                                              <  2:7464482/101‑1 (MQ=255)

    aGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTa                                                                                             <  1:7464420/101‑1 (MQ=255)

       tATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATa                                                                                          <  2:7464241/101‑1 (MQ=255)

        aTGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATAc                                                                                         <  1:7464230/101‑1 (MQ=255)

        aTGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATAc                                                                                         <  2:7464293/101‑1 (MQ=255)

          gTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAgg                                                                                               >  2:7464930/1‑93 (MQ=255)

           tttCCATTCCGGTTAGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTa                                                                                      >  1:7464932/1‑101 (MQ=255)

            ttCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTaa                                                                                     <  1:7464281/101‑1 (MQ=255)

             tCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAAt                                                                                    >  1:7464938/1‑101 (MQ=255)

               cATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGc                                                                                  <  2:7464317/101‑1 (MQ=255)

                aTTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGcc                                                                                 >  1:7464941/1‑101 (MQ=255)

                aTTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGcc                                                                                 >  2:7464942/1‑101 (MQ=255)

                 ttCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCa                                                                                >  1:7464946/1‑101 (MQ=255)

                   ccGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGt                                                                              <  2:7464182/101‑1 (MQ=255)

                   ccGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGt                                                                              <  2:7464358/101‑1 (MQ=255)

                   ccGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGt                                                                              >  1:7464950/1‑101 (MQ=255)

                    cGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGtt                                                                             <  2:7464041/101‑1 (MQ=255)

                    cGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGtt                                                                             <  1:7464150/101‑1 (MQ=255)

                    cGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGtt                                                                             <  2:7464277/101‑1 (MQ=255)

                     ggTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTc                                                                            >  2:7464954/1‑101 (MQ=255)

                     ggTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTc                                                                            <  2:7464410/101‑1 (MQ=255)

                     ggTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTc                                                                            <  2:7464518/101‑1 (MQ=255)

                      gTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATa                                                                                          <  2:7464390/86‑1 (MQ=255)

                      gTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCa                                                                           >  1:7464955/1‑101 (MQ=255)

                       ttGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCAc                                                                          <  1:7464258/101‑1 (MQ=255)

                          ggAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCa                                                                       >  2:7464961/1‑101 (MQ=255)

                           gAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAgg                                                                                               >  2:7464965/1‑76 (MQ=255)

                            aaTAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCAtt                                                                     >  2:7464969/1‑101 (MQ=255)

                            aaTAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCAtt                                                                     >  1:7464968/1‑101 (MQ=255)

                            aaTAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCAtt                                                                     <  1:7464531/101‑1 (MQ=255)

                              tAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTc                                                                   >  1:7464974/1‑101 (MQ=255)

                                ggTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCg                                                                 <  2:7464619/101‑1 (MQ=255)

                                ggTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCg                                                                 >  1:7464978/1‑101 (MQ=255)

                                ggTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCg                                                                 <  1:7464442/101‑1 (MQ=255)

                                   gtCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCat                                                              <  1:7463973/101‑1 (MQ=255)

                                    tCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCatt                                                             <  1:7464521/101‑1 (MQ=255)

                                    tCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCatt                                                             >  1:7464983/1‑101 (MQ=255)

                                    tCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAAAAGCGAGGTAATACGTAATGCCAGTTCACGCATTTcc                                                                  >  2:7464984/1‑96 (MQ=255)

                                      atatTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGa                                                                                                 >  2:7464989/1‑63 (MQ=255)

                                      atatTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCattat                                                           <  1:7464376/101‑1 (MQ=255)

                                       tatTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCattatt                                                          <  2:7464468/101‑1 (MQ=255)

                                       tatTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCattatt                                                          >  2:7464994/1‑101 (MQ=255)

                                        atTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCAttatta                                                         >  1:7464996/1‑101 (MQ=255)

                                        atTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCAttatta                                                         <  2:7464375/101‑1 (MQ=255)

                                        atTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCAttatta                                                         <  1:7464708/101‑1 (MQ=255)

                                         ttCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTaa                                                        >  2:7465001/1‑101 (MQ=255)

                                          tCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAAt                                                       <  2:7464532/101‑1 (MQ=255)

                                           ccAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAAtt                                                      >  2:7465008/1‑101 (MQ=255)

                                           ccAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAAtt                                                      <  2:7464455/101‑1 (MQ=255)

                                            cAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTg                                                     <  1:7464346/101‑1 (MQ=255)

                                            cAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTg                                                     <  1:7464430/101‑1 (MQ=255)

                                            cAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTg                                                     >  1:7465010/1‑101 (MQ=255)

                                             aGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGc                                                    >  2:7465014/1‑101 (MQ=255)

                                               aTTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCat                                                  >  1:7465016/1‑101 (MQ=255)

                                                 tGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCatga                                                <  1:7464305/101‑1 (MQ=255)

                                                 tGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCatga                                                <  2:7464438/101‑1 (MQ=255)

                                                  gCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCatgaa                                               <  1:7464502/101‑1 (MQ=255)

                                                   ccATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCatgaat                                              <  2:7464732/101‑1 (MQ=255)

                                                    cATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCatgaatg                                             <  1:7464393/101‑1 (MQ=255)

                                                    cATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCatgaatg                                             >  1:7465020/1‑101 (MQ=255)

                                                     aTTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAA‑GCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGAt                                           <  1:7464549/101‑1 (MQ=255)

                                                        gTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGc                                                    >  2:7465026/1‑90 (MQ=255)

                                                        gTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTg                                         >  1:7465025/1‑101 (MQ=255)

                                                           cTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTc                                      <  1:7464318/101‑1 (MQ=255)

                                                            tGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTcg                                     >  1:7465030/1‑101 (MQ=255)

                                                             ggCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTcgc                                    <  1:7464475/101‑1 (MQ=255)

                                                             ggCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTcgc                                    <  2:7464556/101‑1 (MQ=255)

                                                               cAGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCgcgc                                  >  2:7465034/1‑101 (MQ=255)

                                                                aGTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGcc                                 >  1:7465038/1‑101 (MQ=255)

                                                                 gTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTg                                         <  1:7464454/92‑1 (MQ=255)

                                                                 gTCAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCa                                >  2:7465041/1‑101 (MQ=255)

                                                                   cAACGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGc                              >  1:7465045/1‑101 (MQ=255)

                                                                    aaCGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCa                             <  2:7464489/101‑1 (MQ=255)

                                                                    aaCGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCa                             >  1:7465049/1‑101 (MQ=255)

                                                                     aCGTCGTATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCAt                            <  1:7464626/101‑1 (MQ=255)

                                                                       gtcgtATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACg                                                                         <  2:7464264/54‑1 (MQ=255)

                                                                       gtcgtATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCAttt                          >  2:7465050/1‑101 (MQ=255)

                                                                       gtcgtATAACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCAttt                          <  1:7464203/101‑1 (MQ=255)

                                                                            ataACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCaa                       >  2:7465059/1‑99 (MQ=255)

                                                                            ataACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATc                     >  1:7465056/1‑101 (MQ=255)

                                                                            ataACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATc                     <  1:7464067/101‑1 (MQ=255)

                                                                             taACGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCt                    <  2:7464606/101‑1 (MQ=255)

                                                                              aaCGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTg                   >  1:7465063/1‑101 (MQ=255)

                                                                                cGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCa                 <  2:7464547/101‑1 (MQ=255)

                                                                                cGGCAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCa                 >  2:7465064/1‑101 (MQ=255)

                                                                                   cAGGAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACt              >  1:7465067/1‑101 (MQ=255)

                                                                                      gAAATCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCAt           >  2:7465074/1‑101 (MQ=255)

                                                                                       aaaTCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCAtt          >  1:7465075/1‑101 (MQ=255)

                                                                                         aTCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTga        <  1:7464581/101‑1 (MQ=255)

                                                                                          tCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCAttt                                                                    <  1:7464799/40‑1 (MQ=39)

                                                                                          tCTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTgac       >  2:7465079/1‑101 (MQ=255)

                                                                                           cTAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTgacg      <  2:7464681/101‑1 (MQ=255)

                                                                                            tAACAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTgacga     >  1:7465082/1‑101 (MQ=255)

                                                                                             aaCAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTgacgac    >  1:7465086/1‑101 (MQ=255)

                                                                                             aaCAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTgacgac    <  1:7464103/101‑1 (MQ=255)

                                                                                               cAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTGACGACTa  <  2:7464428/101‑1 (MQ=255)

                                                                                               cAGCGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTGACGACTa  <  1:7464595/101‑1 (MQ=255)

AACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTATAACGGCAGGAAATCTAACAGTGAGGTAATACGTAATGCCAGTTCACGCATTTCCGCATTATTAATTGCATGAATGATTGCTCGCGCCAGCATTTCAATCTGCAACTCATTGACGACTA  >  NC_000913/4373660‑4373855

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