BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,231,189	A→C	S768R (AGT→CGT)	ygjK →	alpha‑glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,231,189	0	A	C	99.7%	1095.8 / ‑4.6	304	S768R (AGT→CGT)	ygjK	alpha‑glucosidase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (173/130); total (174/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 304 total aligned reads displayed.

CTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTAC  >  NC_000913/3231092‑3231285

cttcCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGt                                                                                                <  1:2397834/100‑1 (MQ=255)

  tcCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGc                                                                                              <  1:5265555/100‑1 (MQ=255)

    cGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCcg                                                                                            >  2:40342/1‑100 (MQ=255)

     ggCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCcgc                                                                                           >  1:458224/1‑100 (MQ=255)

      gCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCcgcg                                                                                          <  2:2093381/100‑1 (MQ=255)

      gCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCcgcg                                                                                          <  2:2198295/100‑1 (MQ=255)

      gCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCcgcg                                                                                          <  1:2599374/100‑1 (MQ=255)

        aCGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCa                                                                                        <  1:885722/100‑1 (MQ=255)

          gCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCACAtt                                                                                      <  1:1998224/100‑1 (MQ=255)

           ccAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCAttt                                                                                     >  1:2829208/1‑100 (MQ=255)

           ccAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCAttt                                                                                     >  1:4912345/1‑100 (MQ=255)

            cAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTg                                                                                    >  1:4788009/1‑100 (MQ=255)

             aaaGGGTTAACCGCCGATGGACCAATTGAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCAtttct                                                                                   >  2:5464700/1‑98 (MQ=255)

             aaaGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGt                                                                                   >  1:4405605/1‑100 (MQ=255)

               aGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGtat                                                                                 >  2:1323680/1‑100 (MQ=255)

                gggTTAACCGCCGATGGACCAGTTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCGCCTGGCGTGCCGCGCATTTGtata                                                                                <  1:5193303/100‑1 (MQ=255)

                gggTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGtata                                                                                <  1:2912660/100‑1 (MQ=255)

                 ggTTCACCGCCGATGGACCAATGCAGGAAAATTACAGCCCGAAGACAGACGAACAGGAAGACGCACCAAATTTATCATGGCGTGCGGAGCATTTGTatat                                                                               >  2:5433169/1‑100 (MQ=255)

                 ggTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTatat                                                                               >  2:3819686/1‑100 (MQ=255)

                 ggTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTatat                                                                               >  1:3578040/1‑100 (MQ=255)

                  gTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATAtg                                                                              >  1:5006712/1‑100 (MQ=255)

                   ttAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAATTTTCTCCTGGCGTGCCGCGCATTTGTATAtgt                                                                             <  2:4954724/100‑1 (MQ=255)

                     aaCCGCCGATGGACCAATTCAGGAAAATTACACCCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATAtgttg                                                                           <  1:2449833/100‑1 (MQ=255)

                      aCCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATAtgttgt                                                                          >  1:2162004/1‑100 (MQ=255)

                      aCCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATAtgttgt                                                                          >  1:124570/1‑100 (MQ=255)

                       ccgccgATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGta                                                                         >  2:36266/1‑100 (MQ=255)

                       ccgccgATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGta                                                                         >  1:8395/1‑100 (MQ=255)

                       ccgccgATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGta                                                                         >  2:78808/1‑100 (MQ=255)

                        cgccgATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGtat                                                                        <  1:1435163/100‑1 (MQ=255)

                        cgccgATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGtat                                                                        <  1:1602656/100‑1 (MQ=255)

                        cgccgATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGtat                                                                        <  2:4987661/100‑1 (MQ=255)

                           cgATGGACCAATTCAGGAAAATTACAACCAGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAAc                                                                     >  1:3676037/1‑100 (MQ=255)

                            gATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACg                                                                    >  1:4078831/1‑100 (MQ=255)

                            gATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACg                                                                    >  2:2557233/1‑100 (MQ=255)

                               ggACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAtt                                                                 <  2:1476305/100‑1 (MQ=255)

                                gACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAttt                                                                >  1:252926/1‑100 (MQ=255)

                                 nccAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAtttt                                                               >  2:4860890/2‑100 (MQ=255)

                                 aCCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAtttt                                                               <  1:947646/100‑1 (MQ=255)

                                 aCCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAtttt                                                               >  1:1827725/1‑100 (MQ=255)

                                 aCCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAtttt                                                               >  1:5377885/1‑100 (MQ=255)

                                  ccAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAttttt                                                              <  1:2337741/100‑1 (MQ=255)

                                  ccAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAttttt                                                              >  2:3649623/1‑100 (MQ=255)

                                   cAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAtttttt                                                             <  1:1615425/100‑1 (MQ=255)

                                   cAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGAtttttt                                                             >  2:3606764/1‑100 (MQ=255)

                                    aaTTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTc                                                            >  2:173903/1‑100 (MQ=255)

                                    aaTTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTc                                                            >  1:1019897/1‑100 (MQ=255)

                                     aTTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTcc                                                           >  1:3504732/1‑100 (MQ=255)

                                     aTTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTcc                                                           <  2:4130478/100‑1 (MQ=255)

                                      ttCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCg                                                          >  2:2578321/1‑100 (MQ=255)

                                       tctggAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGt                                                         <  2:164181/97‑1 (MQ=255)

                                        cAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTa                                                        >  2:4538436/1‑100 (MQ=255)

                                        cAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTa                                                        >  2:2642922/1‑100 (MQ=255)

                                        cAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTa                                                        >  2:1165004/1‑100 (MQ=255)

                                        cAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTa                                                        >  1:3129889/1‑100 (MQ=255)

                                         aGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGTCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTaa                                                       <  2:1905520/100‑1 (MQ=255)

                                         aGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTaa                                                       >  2:532222/1‑100 (MQ=255)

                                          gggaaaTTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAg                                                      >  1:5549251/4‑100 (MQ=255)

                                          ggAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAg                                                      >  1:3878338/1‑100 (MQ=255)

                                          ggAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAg                                                      <  1:1132052/100‑1 (MQ=255)

                                          ggAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAg                                                      <  2:3544969/100‑1 (MQ=255)

                                           gAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGc                                                     <  1:2747635/100‑1 (MQ=255)

                                                  aCAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAAtt                                              >  2:2818948/1‑100 (MQ=255)

                                                    aaCCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTg                                            >  2:4621351/1‑100 (MQ=255)

                                                     aCCCGCTGACAGGCGCACAGCAAGGCGCGCCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGc                                           <  2:5273170/100‑1 (MQ=255)

                                                     aCCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGc                                           <  1:2083754/100‑1 (MQ=255)

                                                     aCCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGc                                           >  1:3288340/1‑100 (MQ=255)

                                                     aCCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGc                                           >  2:4160726/1‑100 (MQ=255)

                                                      cccGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCa                                          <  1:5454061/100‑1 (MQ=255)

                                                       ccGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCaa                                         >  1:4758081/1‑100 (MQ=255)

                                                       ccGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCaa                                         >  2:2883940/1‑100 (MQ=255)

                                                       ccGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCaa                                         >  2:5060619/1‑100 (MQ=255)

                                                               aGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGc                                 <  2:3155142/100‑1 (MQ=255)

                                                               aGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGc                                 <  1:2579353/100‑1 (MQ=255)

                                                               aGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGc                                 <  2:3482276/100‑1 (MQ=255)

                                                                ggCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCg                                >  1:4722208/1‑100 (MQ=255)

                                                                ggCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCg                                >  2:3329697/1‑100 (MQ=255)

                                                                ggCGCACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCg                                >  2:2062920/1‑100 (MQ=255)

                                                                 gcgcACAGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCgg                               >  2:2945852/1‑100 (MQ=255)

                                                                    cacaGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATg                            >  2:3272495/1‑100 (MQ=255)

                                                                    cacaGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATg                            >  2:5137297/1‑100 (MQ=255)

                                                                     acaGCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATgc                           >  2:2445989/1‑100 (MQ=255)

                                                                        gCAAGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATgcggc                        >  1:5513113/1‑100 (MQ=255)

                                                                           aGGCGCACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTa                     <  2:2537351/100‑1 (MQ=255)

                                                                           aGGCGCACCAAATTTCTCCTGGCGTGCCGCGAATTTGTATACGTTGTATAACGATTTTTTCCGTAAGCAAGAAGTTGCAACTGCTGGCGGATGCGGCGTa                     >  2:5228912/1‑100 (MQ=255)

                                                                              cgcACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAAc                  >  1:629898/1‑100 (MQ=255)

                                                                              cgcACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAAc                  <  1:3971497/100‑1 (MQ=255)

                                                                               gcACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACg                 <  2:4693601/100‑1 (MQ=255)

                                                                                cACCAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGc                >  2:4363543/1‑100 (MQ=255)

                                                                                  ccAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCt              >  2:1595771/1‑100 (MQ=255)

                                                                                  ccAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCt              >  1:1696188/1‑100 (MQ=255)

                                                                                   cAAATTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCtt             >  2:1212888/1‑100 (MQ=255)

                                                                                      aTTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATc          <  1:26402/100‑1 (MQ=255)

                                                                                      aTTTCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATc          >  2:124298/1‑100 (MQ=255)

                                                                                       tttCTCCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATcc         >  1:1362959/1‑100 (MQ=255)

                                                                                          ctcCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCACTAATTTGCAACTGCTGGCGGATGCGGCGTACACGCCTTCTCCGcc      >  1:1867653/1‑100 (MQ=255)

                                                                                          ctcCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGcc      >  2:4774800/1‑100 (MQ=255)

                                                                                          ctcCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGcc      >  2:4835585/1‑100 (MQ=255)

                                                                                          ctcCTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGcc      >  1:3416897/1‑100 (MQ=255)

                                                                                            ccTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCt    <  2:2006344/100‑1 (MQ=255)

                                                                                             cTGGCGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTa   >  1:3472546/1‑100 (MQ=255)

CTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTAC  >  NC_000913/3231092‑3231285

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