BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,900,410	A→T	15.1%	P56P (CCT→CCA)	yieL ←	putative xylanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,900,410	0	A	T	15.1%	536.7 / 75.3	211	P56P (CCT→CCA)	yieL	putative xylanase
Reads supporting (aligned to +/- strand): ref base A (81/98); new base T (25/7); total (106/105)
Fisher's exact test for biased strand distribution p-value = 8.76e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

Only 100 of 213 total aligned reads displayed.

GGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGC  >  NC_000913/3900312‑3900508

 gtgtGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAgg                                                                                                  >  2:57253/1‑100 (MQ=255)

      cgcCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaa                                                                                             >  2:101737/1‑100 (MQ=255)

      cgcCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaa                                                                                             <  1:3511957/100‑1 (MQ=255)

           cGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGCaaaaaat                                                                                        <  2:3894666/100‑5 (MQ=255)

            gACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATAGTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGCaaaaaatg                                                                                       <  2:3284417/100‑6 (MQ=255)

              ccAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGCaaaaaatgtc                                                                                     <  2:2705041/100‑8 (MQ=255)

               cataCTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGGCCTGAACCGGACCACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAACTAGCGaa                                                                                    <  1:218647/97‑1 (MQ=255)

               cAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGaa                                                                                    <  2:138108/100‑1 (MQ=255)

                aGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGaaa                                                                                   <  2:2655335/100‑1 (MQ=255)

                 gACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAg                                                                                  <  1:2997606/100‑1 (MQ=255)

                 gACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACGACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAg                                                                                  >  1:981738/1‑100 (MQ=255)

                   cTCCTGTTTCGTCGTTGGTCATCGGGTTTATATTGTCCGGAACCTGGACAACCACCACTCCATAGACATGTTTTGCCTCTGGTGCaaaa                                                                                           <  1:240770/89‑2 (MQ=255)

                   cTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGt                                                                                <  2:3938587/100‑1 (MQ=255)

                   cTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAATAGCGAAAGGt                                                                                <  1:4320904/100‑1 (MQ=255)

                    tgcTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCAGTACAGAGACATTTTTTGCCTCAGGGGCAAAGTAGCGAAAgg                                                                                 <  1:3813993/96‑1 (MQ=255)

                     ccTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCa                                                                              >  1:3652491/1‑100 (MQ=255)

                     ccTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCa                                                                              >  2:1991619/1‑100 (MQ=255)

                        gCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGc                                                                           <  2:2145203/100‑1 (MQ=255)

                         cTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCt                                                                          <  2:1754323/100‑1 (MQ=255)

                         cTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCt                                                                          >  1:1901287/1‑100 (MQ=255)

                          ttCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTg                                                                         <  2:3506153/100‑1 (MQ=255)

                             gTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTa                                                                      >  2:2208225/1‑100 (MQ=255)

                                    gTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCAt                                                               >  2:267912/1‑100 (MQ=255)

                                    gTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCAt                                                               >  1:1774063/1‑100 (MQ=255)

                                    gTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCAt                                                               >  1:535869/1‑100 (MQ=255)

                                    gTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCAt                                                               >  2:48893/1‑100 (MQ=255)

                                      cATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCTTCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTg                                                             >  1:834578/1‑100 (MQ=255)

                                      cATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTg                                                             <  2:2583273/100‑1 (MQ=255)

                                       aTCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGa                                                            >  1:3447301/1‑100 (MQ=255)

                                        tCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGAc                                                           <  2:3122793/100‑1 (MQ=255)

                                         cGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTTTTTCAGTGTCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGAcc                                                          >  1:3026711/1‑100 (MQ=255)

                                           ggTGAATATTGTCCGGAACCGGAACACCCACCAGTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTg                                                        >  1:2642441/1‑100 (MQ=255)

                                           ggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTg                                                        <  2:3556846/100‑1 (MQ=255)

                                           ggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTg                                                        <  2:3642095/100‑1 (MQ=255)

                                             tGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTTCTTTTTTTTCAAAGTAACGAAAGGTCACGCTGTTATCGGCATTGACCTGAg                                                      >  1:1602855/1‑100 (MQ=255)

                                             tGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAg                                                      <  2:2192448/100‑1 (MQ=255)

                                             tGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGCAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAg                                                      <  2:3098103/100‑1 (MQ=255)

                                              gAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGt                                                     <  2:3835442/100‑1 (MQ=255)

                                              gAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGt                                                     <  2:2540734/100‑1 (MQ=255)

                                               aaTATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTc                                                    <  1:2040196/100‑1 (MQ=255)

                                                atatTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTca                                                   <  2:2266218/100‑1 (MQ=255)

                                                 tatTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTcac                                                  <  2:29004/100‑1 (MQ=255)

                                                  atTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGTCATTGACCTGAGTcaca                                                 >  2:3952584/1‑100 (MQ=255)

                                                   ttGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTTCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTCTCGGCATTGACCTGAGTCACAt                                                >  1:3545314/1‑100 (MQ=255)

                                                    tGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATa                                               >  1:3121088/1‑100 (MQ=255)

                                                     gTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATAc                                              <  2:444715/100‑1 (MQ=255)

                                                        cggaaccggaacACCCACCACTACAGAGACATTTTTTTTTTTTTTTTCAAAGTAGCGAAAGGTCACGCTATTATGGGCATTGACCTGAGTCACATACTGt                                           >  1:400512/1‑100 (MQ=255)

                                                        cggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGt                                           <  2:870507/100‑1 (MQ=255)

                                                         ggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGtt                                          <  2:2205853/100‑1 (MQ=255)

                                                         ggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGtt                                          <  1:1235662/100‑1 (MQ=255)

                                                         ggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGtt                                          <  2:1458143/100‑1 (MQ=255)

                                                               ggaacACCCACCACTACAGAGACATTTTTTTTTTTTTTTTCAAAGTAGCGAAAGGTCACGCTGTTATCGGGATTGCCCTGAGTCACATACTGTTTTACCg                                    >  1:1035185/1‑100 (MQ=255)

                                                               ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTCGCGAAAGCTCCCGCTGTTATCGCCATTGACCTGAGTCACATACTGTTTCACCg                                    >  1:428933/1‑100 (MQ=255)

                                                               ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCg                                    <  1:2760412/100‑1 (MQ=255)

                                                               ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCg                                    >  1:1504079/1‑100 (MQ=255)

                                                               ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCg                                    >  2:1158451/1‑100 (MQ=255)

                                                               ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCg                                    <  1:4071287/100‑1 (MQ=255)

                                                                 aacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCggg                                  <  2:4199353/100‑1 (MQ=255)

                                                                  acacCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGa                                 <  1:2289933/100‑1 (MQ=255)

                                                                  acacCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGa                                 >  1:2249697/1‑100 (MQ=255)

                                                                  acacCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGa                                 <  2:2818814/100‑1 (MQ=255)

                                                                     cccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAg                              >  2:1488845/1‑100 (MQ=255)

                                                                     cccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAg                              >  2:1412237/1‑100 (MQ=255)

                                                                      ccaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGt                             <  2:3768442/100‑1 (MQ=255)

                                                                       caccacTACAGAGACATTTTTTTTTTTTTTTTCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGTTAGt                             >  1:1352033/1‑99 (MQ=255)

                                                                       caccacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTg                            <  2:1931910/100‑1 (MQ=255)

                                                                             tACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGTGATAGTGGGTGcc                      <  2:2961695/100‑1 (MQ=255)

                                                                              aCAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCg                     >  2:1840377/1‑100 (MQ=255)

                                                                               cAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCgg                    <  2:941532/100‑1 (MQ=255)

                                                                               cAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCgg                    >  1:1532318/1‑100 (MQ=255)

                                                                                 gagaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTACCGggc                  >  2:34014/1‑100 (MQ=255)

                                                                                  agaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGATCTTAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggct                 >  1:3086255/1‑100 (MQ=255)

                                                                                   gaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctg                >  2:2756625/1‑100 (MQ=255)

                                                                                     gctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggc              >  2:3916231/6‑100 (MQ=255)

                                                                                     gctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggc              >  2:2942062/6‑100 (MQ=255)

                                                                                     gctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggc              >  1:3412839/6‑100 (MQ=255)

                                                                                     cATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggc              <  2:2417852/100‑1 (MQ=255)

                                                                                     cATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggc              <  2:383052/100‑1 (MQ=255)

                                                                                      ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCATGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggct             >  1:1377146/5‑100 (MQ=255)

                                                                                      ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggct             >  2:20467/5‑100 (MQ=255)

                                                                                      ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggct             >  2:1098724/5‑100 (MQ=255)

                                                                                      ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggc              >  1:2795088/5‑99 (MQ=255)

                                                                                      aTTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggct             <  2:3176237/100‑1 (MQ=255)

                                                                                       ttttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGgctggctg            <  2:3025780/100‑1 (MQ=255)

                                                                                       tactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGgctggctg            >  1:2917191/4‑100 (MQ=255)

                                                                                        tttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg           <  2:791945/100‑1 (MQ=255)

                                                                                        tttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg           <  2:4225709/100‑1 (MQ=255)

                                                                                         ttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCtggctggc          <  2:1251456/100‑1 (MQ=255)

                                                                                            tGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTTTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCata       >  1:708145/1‑100 (MQ=255)

                                                                                             gCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCatat      >  1:2968183/1‑100 (MQ=255)

                                                                                             gCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCatat      >  1:4239889/1‑100 (MQ=255)

                                                                                              ccccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATc     >  2:148031/1‑100 (MQ=255)

                                                                                              ccccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATc     >  2:1288569/1‑100 (MQ=255)

                                                                                               cccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATcc    >  1:3931378/1‑100 (MQ=255)

                                                                                               cccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATcc    <  1:2472646/100‑1 (MQ=255)

                                                                                               cccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATcc    <  2:1556058/100‑1 (MQ=255)

                                                                                               cccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATcc    >  2:735708/1‑100 (MQ=255)

                                                                                               cccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATcc    <  2:1590063/100‑1 (MQ=255)

                                                                                               cccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATc     <  2:4133431/99‑1 (MQ=255)

                                                                                                ccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCg   <  1:3761452/100‑1 (MQ=255)

GGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGC  >  NC_000913/3900312‑3900508

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