BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,900,410	A→T	5.4%	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	5.4%	305.7 / 6.5	112	P56P (CCT→CCA)	yieL	putative xylanase
Reads supporting (aligned to +/- strand): ref base A (54/50); new base T (4/2); total (59/53)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

Only 100 of 114 total aligned reads displayed.

TACACCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGC  >  NC_000913/3900263‑3900559

tACACCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAgg                                                                                                                                                     <  1:451680/150‑1 (MQ=255)

  cacCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGACCCCCCCCCACTACAGAGCCATTTTTTGCCCCAgggg                                                                                                                                                   >  2:1365933/1‑150 (MQ=255)

  cacCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAgggg                                                                                                                                                   <  1:1214931/150‑1 (MQ=255)

     cATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCATGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaa                                                                                                                                                <  1:1077059/150‑1 (MQ=255)

     cATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaa                                                                                                                                                >  1:401507/1‑150 (MQ=255)

      aTCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaaa                                                                                                                                               <  1:62491/150‑1 (MQ=255)

       tCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAg                                                                                                                                              >  2:462742/1‑150 (MQ=255)

        cAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGt                                                                                                                                             >  2:834905/1‑150 (MQ=255)

          aCATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGCaaaaaat                                                                                                                                           <  2:1597792/150‑5 (MQ=255)

              tGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGaa                                                                                                                                       >  2:1686094/1‑150 (MQ=255)

               gAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGCaaaaaatg                                                                                                                                          <  2:1653756/146‑6 (MQ=255)

               gAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGCaaaaaatg                                                                                                                                          >  1:1653756/1‑141 (MQ=255)

                 aacactACGCGGACAGGGTGCCTATCAGGAGAGGTGAGCGCAACGACCAGACTCCTGCTTCGTCTTTGGTCATCGTGTGAATATTGTCCGGACCCGGAACACCCACCAATACAGAGACATTTTTTGCCCCAGCTGCAAAGTAGCGAAAgg                                                                                                                                    <  2:265679/145‑1 (MQ=255)

                 aaaaaTACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAgg                                                                                                                                    <  1:478572/150‑1 (MQ=255)

                   aaaTACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTATTGCCCCAGGGGCAAAGTAGCGAAAGGTc                                                                                                                                  <  1:1561791/150‑1 (MQ=255)

                        cTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCt                                                                                                                             >  1:506627/1‑150 (MQ=255)

                              aaGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCTTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGTAACAACAACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATc                                                                                                                       <  1:647067/149‑1 (MQ=255)

                                      cTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACCCCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGa                                                                                                               >  1:1260141/1‑150 (MQ=255)

                                      cTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGGCTTTGGTCATCGGGTGAATATTGTCCGGCACCGGAACACCCACCACTACAGAGACATTTTTTGACCCCGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGa                                                                                                               >  2:758047/1‑150 (MQ=255)

                                       tttCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGAc                                                                                                              >  1:818850/1‑150 (MQ=255)

                                          cAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGATCACCCACCACTACAGAGACATTTTTTGCCCCAGGGACAAAGTAGCGAAAGGTCACGCTGTTATCGGCAAAGACCTg                                                                                                           >  2:706775/1‑150 (MQ=255)

                                          cAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTg                                                                                                           <  1:1340012/150‑1 (MQ=255)

                                           aGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATAGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGa                                                                                                          >  2:29480/1‑150 (MQ=255)

                                            ggATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAg                                                                                                         >  2:1166472/1‑150 (MQ=255)

                                            ggATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAg                                                                                                         <  2:567193/150‑1 (MQ=255)

                                              aTAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTc                                                                                                       <  1:814826/150‑1 (MQ=255)

                                               tAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAACGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTca                                                                                                      <  1:497151/150‑1 (MQ=255)

                                                aGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGGCACGCGGTTATCGGCATTGACCTGAGTcac                                                                                                     <  2:237595/150‑1 (MQ=255)

                                                 ggTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTcaca                                                                                                    >  2:1642770/1‑150 (MQ=255)

                                                 ggTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTcaca                                                                                                    >  1:1273840/1‑150 (MQ=255)

                                                 ggTGGGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGACAGGTCACGCTGTTATCGGCATTGACCTGAGTcaca                                                                                                    <  2:1161319/150‑1 (MQ=255)

                                                  gtgtGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACAt                                                                                                   >  2:1221518/1‑150 (MQ=255)

                                                   tgtgCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATa                                                                                                  <  1:1428792/150‑1 (MQ=255)

                                                        gcCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGtt                                                                                             <  1:816318/150‑1 (MQ=255)

                                                         ccACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGttt                                                                                            <  2:1378945/150‑1 (MQ=255)

                                                             gACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCAcc                                                                                        <  1:33815/150‑1 (MQ=255)

                                                               ccAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCgg                                                                                      >  2:86276/1‑150 (MQ=255)

                                                               ccAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCgg                                                                                      >  2:12137/1‑150 (MQ=255)

                                                               ccAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCgg                                                                                      >  1:1139038/1‑150 (MQ=255)

                                                                  gACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGAt                                                                                   <  1:1516171/150‑1 (MQ=255)

                                                                     tCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGt                                                                                >  1:1513739/1‑150 (MQ=255)

                                                                        tGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTggg                                                                             <  1:1166472/150‑1 (MQ=255)

                                                                             cGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCg                                                                        <  1:159009/150‑1 (MQ=255)

                                                                                cTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggc                                                                     >  1:1545867/1‑150 (MQ=255)

                                                                                cTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTATTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggc                                                                     >  1:1399455/1‑150 (MQ=255)

                                                                                     gTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCCGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggct                                                                <  2:1327990/150‑1 (MQ=255)

                                                                                        aTCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCtggctggc                                                             <  2:441025/150‑1 (MQ=255)

                                                                                         tCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCa                                                            <  1:291762/150‑1 (MQ=255)

                                                                                         tCGGGTGAATATTGGCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGGCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCa                                                            <  1:710139/150‑1 (MQ=255)

                                                                                           gggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCata                                                          <  2:811911/150‑1 (MQ=255)

                                                                                           gggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCAAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCata                                                          >  1:232748/1‑150 (MQ=255)

                                                                                            ggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCatat                                                         >  1:736817/1‑150 (MQ=255)

                                                                                              tGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATcc                                                       <  2:506627/150‑1 (MQ=255)

                                                                                                  tatTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCgctg                                                   >  1:1678503/1‑150 (MQ=255)

                                                                                                    ttGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCgctgct                                                 >  2:1658641/1‑150 (MQ=255)

                                                                                                      gTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTAt                                               >  1:406911/1‑150 (MQ=255)

                                                                                                      gTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTAt                                               >  2:1645251/1‑150 (MQ=255)

                                                                                                        ccggaaccggaaCACCCACCACTACAGCGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCGATgg                                             >  1:1068113/1‑150 (MQ=255)

                                                                                                         ctgaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCACGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGc                                            >  2:1557901/3‑150 (MQ=255)

                                                                                                         cggaaccggaacACCCACCACTACAGAGCCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGCCCTGAGTCACATACTGTGTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCTTATCCGCTGCTATGGc                                            >  1:291588/1‑150 (MQ=255)

                                                                                                         cggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGc                                            <  1:1616875/150‑1 (MQ=255)

                                                                                                         cggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGc                                            <  2:1011254/150‑1 (MQ=255)

                                                                                                          ggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGcc                                           >  1:1211452/1‑150 (MQ=255)

                                                                                                            aaccggcacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCa                                         >  1:839910/1‑150 (MQ=255)

                                                                                                            aaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGcccc                                         >  2:443787/1‑149 (MQ=255)

                                                                                                             accggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCAc                                        >  1:371127/1‑150 (MQ=255)

                                                                                                              ccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACt                                       >  2:587051/1‑150 (MQ=255)

                                                                                                               cggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTg                                      <  2:1217009/150‑1 (MQ=255)

                                                                                                                 gaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAg                                    >  2:133466/1‑150 (MQ=255)

                                                                                                                  aacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGc                                   <  1:1664144/150‑1 (MQ=255)

                                                                                                                    cacccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAg                                 <  2:1513739/150‑1 (MQ=255)

                                                                                                                    cacccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAg                                 >  1:454831/1‑150 (MQ=255)

                                                                                                                     acccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGa                                >  2:1651681/1‑150 (MQ=255)

                                                                                                                     acccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCGGAGTCACATACTGTTTCACCGGGATAGGGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGa                                <  1:1302626/150‑1 (MQ=255)

                                                                                                                          ccacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATAc                           <  1:609891/150‑1 (MQ=255)

                                                                                                                          ccacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATAc                           <  1:1291865/150‑1 (MQ=255)

                                                                                                                            acTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCg                         >  1:218719/1‑150 (MQ=255)

                                                                                                                               aCAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTg                      >  2:3303/1‑150 (MQ=255)

                                                                                                                               aCAGAGACATTTATTCCCCCAGGGGCAAAGTAGCGAAAGGCCACGCGGTTATCGGCCTTGCCCTGAGTCACATACTGTCTCACCGGGATAGTGGGGGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGcgg                      <  2:266626/150‑3 (MQ=255)

                                                                                                                                 agagACATTTTTTGCCCCAGGTGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGc                    <  1:1061305/150‑1 (MQ=255)

                                                                                                                                      cATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAg               <  1:1365933/150‑1 (MQ=255)

                                                                                                                                      cATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAATCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAg               <  1:694043/150‑1 (MQ=255)

                                                                                                                                       ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATAGc              >  1:27439/5‑150 (MQ=255)

                                                                                                                                       aTTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGc              <  1:834905/150‑1 (MQ=255)

                                                                                                                                       aTTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTATCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGc              <  2:170127/150‑1 (MQ=255)

                                                                                                                                        ttttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGcc             >  1:149176/1‑150 (MQ=255)

                                                                                                                                         tttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCa            <  2:1260141/150‑1 (MQ=255)

                                                                                                                                         actttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCa            >  1:1185180/3‑150 (MQ=255)

                                                                                                                                           tttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGc          <  1:78999/150‑1 (MQ=255)

                                                                                                                                             tGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGt        >  1:1521423/1‑150 (MQ=255)

                                                                                                                                             tGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGt        <  1:1652101/150‑1 (MQ=255)

                                                                                                                                              gCCCCAGGGGCAAAGTAGCGTCAGGTCACGCGGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCGGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGtt       <  2:114485/150‑1 (MQ=255)

                                                                                                                                              cccccAGGCGCAAAGTAGCGAAAGGGCACGCTGTTATCGGCAATGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGtt       <  2:1681019/149‑1 (MQ=255)

                                                                                                                                               ccccTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTa      >  2:398697/1‑150 (MQ=255)

                                                                                                                                                cccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTACTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaa     >  2:681672/1‑150 (MQ=255)

                                                                                                                                                 ccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaaa    >  2:1157385/1‑150 (MQ=255)

                                                                                                                                                 ccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaaa    >  2:759598/1‑150 (MQ=255)

                                                                                                                                                  cAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAag   >  1:1467718/1‑150 (MQ=255)

                                                                                                                                                   aGGTGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagc  <  2:532189/150‑1 (MQ=255)

                                                                                                                                                   aGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTTGCAATAGCCAGCGTTAAagc  <  1:758047/150‑1 (MQ=255)

TACACCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGC  >  NC_000913/3900263‑3900559

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