BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,900,410	A→T	15.7%	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.7%	405.6 / 51.6	157	P56P (CCT→CCA)	yieL	putative xylanase
Reads supporting (aligned to +/- strand): ref base A (54/77); new base T (22/3); total (76/81)
Fisher's exact test for biased strand distribution p-value = 1.52e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

GGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCG  >  NC_000913/3900312‑3900507

ggTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAg                                                                                                  >  2:565753/1‑100 (MQ=255)

ggTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAg                                                                                                  >  1:350631/1‑100 (MQ=255)

ggTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAg                                                                                                  >  1:1399709/1‑100 (MQ=255)

ggTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAg                                                                                                  >  1:983856/1‑100 (MQ=255)

 gtgtGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAg                                                                                                  >  1:131441/1‑99 (MQ=255)

   gtgCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAgggg                                                                                               >  2:1627937/1‑100 (MQ=255)

    tgCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGc                                                                                              >  2:614042/1‑100 (MQ=255)

    tgCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGc                                                                                              >  2:3490842/1‑100 (MQ=255)

    tgCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGc                                                                                              >  1:1961420/1‑100 (MQ=255)

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

        ccACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAg                                                                                          >  1:3275234/1‑100 (MQ=255)

         cACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGTCCCAGGGGCAAAGt                                                                                         >  2:2162991/1‑100 (MQ=255)

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

            gACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGc                                                                                      <  2:2628650/100‑1 (MQ=255)

              ccAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGa                                                                                    >  2:2381274/1‑100 (MQ=255)

              ccAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGa                                                                                    <  1:2932611/100‑1 (MQ=255)

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

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

                 gACTCCAGCTTCGTCTTTGGTCATCGGGAGAATATAGTACGAAAACGGAACACCCACCACTACACAAAAATTTTCTGCCCCTGGGGaaaaaaatgtct                                                                                   <  1:3193038/98‑9 (MQ=255)

                  aCTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAgg                                                                                <  2:1211098/100‑1 (MQ=255)

                      cTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCAc                                                                            >  1:1202232/1‑100 (MQ=255)

                      cTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGACACAGGGGCAAAGTAGCGAAAGGTCAc                                                                            >  1:615060/1‑100 (MQ=255)

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

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

                             gTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGGAAAGTAGCGAAAGGTCACGCTGTTa                                                                     >  2:1075128/1‑100 (MQ=255)

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

                              tCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGAAGCGAAAGGTCACGCTGTTAt                                                                    <  1:2562209/100‑1 (MQ=255)

                               cTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATc                                                                   >  1:2376302/1‑100 (MQ=255)

                                 ttGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCgg                                                                 <  1:809912/100‑1 (MQ=255)

                                 ttGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACGACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCgg                                                                 <  2:2971583/100‑1 (MQ=255)

                                  tGGGCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTAACGCTGTTATCGGc                                                                >  1:2936103/1‑100 (MQ=255)

                                    gTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCAt                                                              <  1:3432506/100‑1 (MQ=255)

                                      cATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTg                                                            >  1:58122/1‑100 (MQ=255)

                                      cATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTg                                                            >  1:3184221/1‑100 (MQ=255)

                                      cATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTg                                                            >  1:57991/1‑100 (MQ=255)

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

                                       aTCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGa                                                           <  2:1802300/100‑1 (MQ=255)

                                        tCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTTTTATTTTTTCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGAc                                                          >  1:2281284/1‑100 (MQ=255)

                                         cGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGAcc                                                         <  1:2928674/100‑1 (MQ=255)

                                          gggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCt                                                        >  2:2674902/1‑100 (MQ=255)

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

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

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

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

                                              gAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGt                                                    >  2:78738/1‑100 (MQ=255)

                                                atttgttCCCGAACCGGAACACCCATCACTACATCTCCATTTTTTGCCCCAGGGGCAAAGTAGCTCAAGGTCACGCTGTTATCGTCATTGACCTg                                                       <  1:2703808/89‑1 (MQ=255)

                                                     gTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTTTTTTTTTTTCAAAGTAGCCAAAGGTCACGTTGTTATCGTCATTCCGCTGTGTTACATaa                                             >  1:2935901/1‑99 (MQ=255)

                                                       ccggaaccggaaCACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTg                                           >  2:3116159/1‑100 (MQ=255)

                                                        cggaaccggaacACCCACCACTACAGAGACATTTTTTTTTTTTTTTTTAAAGTAGCAAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGt                                          >  1:1290807/1‑100 (MQ=255)

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

                                                        cggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGt                                          >  2:1583290/1‑100 (MQ=255)

                                                         ggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGtt                                         >  2:2071739/1‑100 (MQ=255)

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

                                                         ggaaccggaacACCCACCACTACAGAGACATTTTGTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGtt                                         <  1:3005920/100‑1 (MQ=255)

                                                          gaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGtt                                         <  1:2506748/99‑1 (MQ=255)

                                                           aaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGtt                                         <  1:845983/98‑1 (MQ=255)

                                                           aaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTc                                       <  1:2474828/100‑1 (MQ=255)

                                                            accggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCa                                      <  2:2843932/100‑1 (MQ=255)

                                                             ccggaacACCCACCACTACAGAGACATTTTTTGCTTTTTGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGCCCTGAGTCACATACTGTTTCAc                                     >  1:2993293/1‑100 (MQ=255)

                                                             ccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCAc                                     <  2:629722/100‑1 (MQ=255)

                                                              cggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCAcc                                    >  2:331021/1‑100 (MQ=255)

                                                               ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATTGGCATTGACCTGAGTCACATACTGTTTCACCg                                   <  1:716528/100‑1 (MQ=255)

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

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

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

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

                                                                        accacTACAGAGACATTTTTTGCTCTAGTGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTgg                          >  1:71543/1‑100 (MQ=255)

                                                                        accacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTgg                          <  1:1789593/100‑1 (MQ=255)

                                                                        accacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTgg                          <  1:2702365/100‑1 (MQ=255)

                                                                         ccacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTggg                         <  2:1072724/100‑1 (MQ=255)

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

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

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

                                                                                agagACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCggg                  >  2:1608618/1‑100 (MQ=255)

                                                                                 gagaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggc                 <  2:2917206/100‑1 (MQ=255)

                                                                                   gaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTGACATACTGTTTCACCGGGATAGTGGGTGCCGggctg               <  2:2508893/100‑1 (MQ=255)

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

                                                                                   gaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctg               <  2:536603/100‑1 (MQ=255)

                                                                                   gaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctg               <  2:1202239/100‑1 (MQ=255)

                                                                                    aCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctgg              <  2:3275257/100‑1 (MQ=255)

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

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

                                                                                      ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCCCATACTGTTTCACCGGGATAGTGGGTGCCGggctggct            >  1:1541866/5‑100 (MQ=255)

                                                                                      ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggct            >  1:2635265/5‑100 (MQ=255)

                                                                                      ctacttccccccTGGAGCAAAGTAGCGAAAAAACACACAGTTACCCGCGTAGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggct            >  1:434112/9‑100 (MQ=255)

                                                                                      aTTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggct            >  2:3055508/1‑100 (MQ=255)

                                                                                       ttttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGgctggctg           >  2:1280284/1‑100 (MQ=255)

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

                                                                                        actttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACGTGAGTCACAGACTGTATCACCGGGATAGTGGGTGCCGGGctggctgg          >  1:2403250/3‑100 (MQ=255)

                                                                                        actttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg          >  2:334533/3‑100 (MQ=255)

                                                                                        actttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg          >  1:1304270/3‑100 (MQ=255)

                                                                                        actttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg          >  2:3549771/3‑100 (MQ=255)

                                                                                        actttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg          >  2:1355208/3‑100 (MQ=255)

                                                                                        actttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg          >  2:943505/3‑100 (MQ=255)

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

                                                                                          tttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCa        >  2:2690396/1‑100 (MQ=255)

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

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

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

                                                                                                ctaGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCg  <  1:2148005/98‑1 (MQ=255)

GGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCG  >  NC_000913/3900312‑3900507

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