BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,411,637	Δ1 bp	100%	coding (453/651 nt)	yfbT ←	sugar phosphatas

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,411,634	0	C	.	63.7%	147.2 / 263.7	146	coding (456/651 nt)	yfbT	sugar phosphatas
Reads supporting (aligned to +/- strand): ref base C (29/24); new base . (42/51); total (71/75)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

TGCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCC  >  NC_000913/2411535‑2411728

tGCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGc                                                                                                >  1:3150721/1‑100 (MQ=255)

tGCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGc                                                                                                >  2:732871/1‑100 (MQ=255)

tGCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGc                                                                                                >  1:1058619/1‑100 (MQ=255)

    ggAGCGTTACCCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCca                                                                                            <  1:134485/100‑1 (MQ=255)

      aGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGcccagca                                                                                          <  1:468202/100‑5 (MQ=255)

        cGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCcagcag                                                                                        >  2:556597/1‑100 (MQ=255)

         gTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCagcagc                                                                                       >  2:928523/1‑100 (MQ=255)

         gTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCt                                                                                      >  2:17468/1‑100 (MQ=255)

             aCCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTgcgc                                                                                  <  2:861247/100‑1 (MQ=255)

             aCCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTgcgc                                                                                  >  1:1282249/1‑100 (MQ=255)

              ccGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTgcgc                                                                                  >  2:3233911/1‑100 (MQ=255)

              ccGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTgcgc                                                                                  >  1:1634654/1‑100 (MQ=255)

              ccGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTgcgcg                                                                                 <  1:2981759/100‑1 (MQ=255)

              ccGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTgcgcg                                                                                 >  2:624686/1‑100 (MQ=255)

                gCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGcc                                                                               <  2:2557298/100‑1 (MQ=255)

                gCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGcc                                                                               <  2:1142292/100‑1 (MQ=255)

                 cAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGcc                                                                               >  2:772634/1‑100 (MQ=255)

                  aaTGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCt                                                                              <  2:3021907/100‑1 (MQ=255)

                  aaTGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTa                                                                             <  1:1335541/100‑1 (MQ=255)

                     gacatgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAAc                                                                           <  2:739637/100‑1 (MQ=255)

                     gacatgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACCCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACa                                                                          >  1:2696480/1‑100 (MQ=255)

                      acatgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACa                                                                          <  1:2433369/100‑1 (MQ=255)

                      acatgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACa                                                                          <  2:298110/100‑1 (MQ=255)

                       catgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGa                                                                        >  1:2253576/1‑100 (MQ=255)

                         tgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATa                                                                      <  1:140655/100‑1 (MQ=255)

                          gacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATAc                                                                     >  1:442453/1‑100 (MQ=255)

                          gacaACCCGCCACCAGGCCAGAAAGCCCGCCAGCGGGAGCATCTTCCACCACCACACACCCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATAc                                                                     >  1:2887399/1‑100 (MQ=255)

                           acaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACg                                                                    <  2:3350348/100‑1 (MQ=255)

                           acaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACg                                                                    >  2:1975183/1‑100 (MQ=255)

                             aaCCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGc                                                                   <  2:1634661/100‑1 (MQ=255)

                              aCCCGCCGCTAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCAt                                                                 >  1:2449811/1‑100 (MQ=255)

                                ccgccgCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATc                                                                >  1:2532567/1‑100 (MQ=255)

                                ccgccgCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCa                                                               <  1:81340/100‑1 (MQ=255)

                                 cgccgcCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAg                                                              <  1:388476/100‑1 (MQ=255)

                                     gccaggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTc                                                          >  2:1257743/1‑100 (MQ=255)

                                      ccaggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTc                                                          >  2:3116290/1‑100 (MQ=255)

                                      ccaggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCt                                                         >  2:2446799/1‑100 (MQ=255)

                                       caggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCt                                                         <  2:130197/100‑1 (MQ=255)

                                       caggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTg                                                        >  1:721852/1‑100 (MQ=255)

                                        aggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTgg                                                       <  2:2362864/100‑1 (MQ=255)

                                        aggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTgg                                                       <  2:2364533/100‑1 (MQ=255)

                                         ggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                      >  2:2194944/1‑100 (MQ=255)

                                         ggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                      <  2:2491742/100‑1 (MQ=255)

                                         ggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                      >  2:2194935/1‑100 (MQ=255)

                                         ggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                      >  1:85134/1‑100 (MQ=255)

                                          gccagAAAGCACGCCAGCGGGGGCATCTACCACCACCACCCACTCCTGCGGCGCCAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTc                                                          <  1:308289/95‑1 (MQ=255)

                                          gccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                      >  2:1932174/1‑100 (MQ=255)

                                          gccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                      >  1:720048/1‑100 (MQ=255)

                                           ccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGtt                                                     >  1:282487/1‑100 (MQ=255)

                                           ccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGttt                                                    >  1:3352862/1‑100 (MQ=255)

                                           ccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGttt                                                    >  1:177364/1‑100 (MQ=255)

                                            cagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGttt                                                    <  1:290596/100‑1 (MQ=255)

                                             agAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTc                                                  <  2:809774/100‑1 (MQ=255)

                                              gAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTcc                                                 >  1:1141074/1‑100 (MQ=255)

                                                aaGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCcg                                                >  2:2668663/1‑100 (MQ=255)

                                                 aGCACGCCAGCGGGAGCATCTTCCACCACCACTCACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCcgcg                                              <  2:2018281/100‑1 (MQ=255)

                                                 aGCACGCCAGCGGGAGCATCTTCCACCACCACACAGTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCGAACAGATACGCATCAGGTTCTGGTTTTCcgc                                               <  1:518330/99‑1 (MQ=255)

                                                   cACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCgcgc                                             >  1:1263312/1‑100 (MQ=255)

                                                    aCGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCt                                            <  2:2314345/100‑1 (MQ=255)

                                                     cGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCtt                                           <  2:2314372/100‑1 (MQ=255)

                                                     cGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCtt                                           >  2:2115458/1‑100 (MQ=255)

                                                     cGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTc                                          <  2:3146203/100‑1 (MQ=255)

                                                      gCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATGCGCATCAGGTTCTGGTTTTCCGCGCTTc                                          <  2:84722/100‑1 (MQ=255)

                                                      gCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCa                                         >  1:1987241/1‑100 (MQ=255)

                                                       ccAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCAc                                        >  2:1744753/1‑100 (MQ=255)

                                                           cGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCg                                     >  1:336882/1‑100 (MQ=255)

                                                           cGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGc                                    <  2:447447/100‑1 (MQ=255)

                                                            gggTGCATCGTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCt                                   <  1:1726601/100‑1 (MQ=255)

                                                            gggAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGc                                    <  2:3190963/100‑1 (MQ=255)

                                                            gggAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGc                                    <  2:1110814/100‑1 (MQ=255)

                                                             ggAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCt                                   >  2:333389/1‑100 (MQ=255)

                                                             ggAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCt                                   >  2:964769/1‑100 (MQ=255)

                                                             ggAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTc                                  <  1:668233/100‑1 (MQ=255)

                                                                gCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGc                               >  2:438220/1‑100 (MQ=255)

                                                                gCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGc                               >  2:1836183/1‑100 (MQ=255)

                                                                gCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGc                               <  1:257925/100‑1 (MQ=255)

                                                                 cATCTTCCACCACCACACACTCCTGCGGTGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGc                               <  2:693211/100‑1 (MQ=255)

                                                                 cATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGc                               <  2:3268924/99‑1 (MQ=255)

                                                                  aTCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCCTCAGGTTCTGGTTTTCCCCGCTTCACTATCTCAg                                >  1:204195/1‑98 (MQ=255)

                                                                   tCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCgg                             >  1:2210903/1‑100 (MQ=255)

                                                                      tCCACCACCACACCCTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTAc                         <  2:283573/100‑1 (MQ=255)

                                                                      tCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTa                          >  2:687371/1‑100 (MQ=255)

                                                                      tCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTAc                         <  2:1282253/100‑1 (MQ=255)

                                                                        caccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACaa                       >  2:971974/1‑100 (MQ=255)

                                                                         accaccacACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACaa                       >  1:3200524/1‑100 (MQ=255)

                                                                          ccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                     >  2:2104192/1‑100 (MQ=255)

                                                                          ccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                     <  2:1192981/100‑1 (MQ=255)

                                                                            accacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAacac                   >  1:581122/1‑100 (MQ=255)

                                                                             ccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACAcc                  <  2:1020417/100‑1 (MQ=255)

                                                                               acacacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACctc                >  1:563967/1‑100 (MQ=255)

                                                                                cacacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACctct               <  1:188746/100‑1 (MQ=255)

                                                                                 acacTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACctct               <  1:1925640/100‑1 (MQ=255)

                                                                                    cTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTg           <  2:3169248/100‑1 (MQ=255)

                                                                                     tCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTg           <  1:2339956/100‑1 (MQ=255)

                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCgg        >  1:680217/1‑100 (MQ=255)

                                                                                        tGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCggg       <  2:1623209/100‑1 (MQ=255)

                                                                                        tGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCggg       <  1:2843205/100‑1 (MQ=255)

                                                                                           ggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAg    >  1:1754150/1‑100 (MQ=255)

                                                                                            gcgcAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGc   <  1:209225/100‑1 (MQ=255)

                                                                                             cgcAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTGCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGc   >  1:1474314/1‑100 (MQ=255)

TGCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCC  >  NC_000913/2411535‑2411728

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