BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,112,040	T→C	A110A (GCA→GCG)	b2972 ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,112,040	0	T	C	99.7%	956.6 / NA	296	A110A (GCA→GCG)	b2972	putative peptidase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (142/153); minor base A (1/0); total (143/153)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

ACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAG  >  NC_000913/3111950‑3112129

 cTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCt                                                                                        <  2:2549279/93‑1 (MQ=255)

  tCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGTCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTg                                                                                       >  1:5876768/1‑93 (MQ=255)

  ccgcgGCAATCAGCCAGGCGGATAATATCATCGCCGCCAGCCCCCACCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTg                                                                                       <  2:104394/92‑1 (MQ=255)

   cgcgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGt                                                                                      <  2:3282643/93‑1 (MQ=255)

   cgcgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGt                                                                                      <  1:277003/93‑1 (MQ=255)

    gcgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGAGCAGACTCGCCAGCAAAAAAGCAAGCGCTGTc                                                                                     <  1:5341642/93‑1 (MQ=255)

    gcgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTc                                                                                     >  2:6277365/1‑93 (MQ=255)

    gcgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTc                                                                                     <  1:2340731/93‑1 (MQ=255)

    gcgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCACCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTc                                                                                     >  2:1941189/1‑93 (MQ=255)

       gCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAAt                                                                                  <  2:5243065/93‑1 (MQ=255)

       gCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAAt                                                                                  <  2:5390348/93‑1 (MQ=255)

       gCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAAt                                                                                  <  2:5137535/93‑1 (MQ=255)

          aTCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAAc                                                                               <  1:2363297/93‑1 (MQ=255)

          aTCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAAc                                                                               <  2:1731314/93‑1 (MQ=255)

              gCCAGGCGGATAATATCATCACCGCCAGCCCCCCTCCACTTTCTGGCCAGACCAGACTCGCCCGCAAAAAAGCAAGCGCTGTCAATAACTccc                                                                           >  1:4399677/1‑92 (MQ=255)

              gCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTcca                                                                           <  2:247038/93‑1 (MQ=255)

                 aGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTccacca                                                                        >  2:1601615/1‑93 (MQ=255)

                  ggCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAg                                                                       >  1:1008976/1‑93 (MQ=255)

                  ggCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAg                                                                       >  2:2403728/1‑93 (MQ=255)

                    cGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCg                                                                     <  2:3582640/93‑1 (MQ=255)

                        taataTCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATa                                                                 <  2:2755972/93‑1 (MQ=255)

                          atatCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGCCCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAAc                                                               <  2:5621229/93‑1 (MQ=255)

                           tatCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACg                                                              <  1:5941155/93‑1 (MQ=255)

                             tcatcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCt                                                            >  2:26760/1‑93 (MQ=255)

                             tcatcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCt                                                            <  2:5584568/93‑1 (MQ=255)

                             tcatcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCt                                                            <  2:2171628/93‑1 (MQ=255)

                             tcatcaCCGCCAGCCCCCACCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCt                                                            >  1:1612660/1‑93 (MQ=255)

                              catcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCtt                                                           <  2:3105331/93‑1 (MQ=255)

                               atcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTg                                                          >  2:4483886/1‑93 (MQ=255)

                                tcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGc                                                         >  2:3485820/1‑93 (MQ=255)

                                 ccccGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCt                                                        <  1:3307248/91‑1 (MQ=255)

                                 caCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCt                                                        <  2:2042154/93‑1 (MQ=255)

                                 caCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCt                                                        >  1:2483853/1‑93 (MQ=255)

                                  aCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTTCACCAGCAGCTAACGCTTGCta                                                       >  1:6252371/1‑92 (MQ=255)

                                  aCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTg                                                       >  1:4009587/1‑93 (MQ=255)

                                   ccGCCAGCCACCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGa                                                      >  1:1862809/1‑93 (MQ=255)

                                    cGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAt                                                     <  1:699114/93‑1 (MQ=255)

                                    cGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAt                                                     <  1:1137937/93‑1 (MQ=255)

                                     gCCAGCCCCCATCCACTTTCTGGCCCGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtt                                                    <  2:2965573/93‑1 (MQ=255)

                                     cccAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtt                                                    <  2:3806436/92‑1 (MQ=255)

                                      ccAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAttt                                                   <  2:6724095/93‑1 (MQ=255)

                                      ccAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAttt                                                   <  2:1360230/93‑1 (MQ=255)

                                       cAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtttt                                                  <  1:4416818/93‑1 (MQ=255)

                                         gCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCg                                                <  2:6149216/93‑1 (MQ=255)

                                           ccccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGcc                                              >  2:5213763/1‑93 (MQ=255)

                                           ccccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGcc                                              >  1:829669/1‑93 (MQ=255)

                                            cccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCt                                             <  1:1657087/93‑1 (MQ=255)

                                            cccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCt                                             >  1:2263690/1‑93 (MQ=255)

                                             ccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTg                                            <  2:1203449/93‑1 (MQ=255)

                                              cATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGa                                           >  2:4880654/1‑93 (MQ=255)

                                               aTCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGAc                                          <  1:1464658/93‑1 (MQ=255)

                                                tCCACTTTCTGGCTAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACa                                         >  1:6094688/1‑93 (MQ=255)

                                                tCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACa                                         >  1:3500048/1‑93 (MQ=255)

                                                  cACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGt                                       >  2:4604274/1‑93 (MQ=255)

                                                   ccTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTc                                      <  1:3168015/92‑1 (MQ=255)

                                                   aCTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTc                                      >  1:5706075/1‑93 (MQ=255)

                                                   aCTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTc                                      >  1:3424412/1‑93 (MQ=255)

                                                      ttCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTcgcg                                   <  2:3298290/93‑1 (MQ=255)

                                                        cTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGc                                 <  1:4259946/93‑1 (MQ=255)

                                                         tGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCa                                >  2:4832680/1‑93 (MQ=255)

                                                          ggCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAg                               >  1:3055960/1‑93 (MQ=255)

                                                              agaccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGcc                           >  1:4249566/1‑93 (MQ=255)

                                                                accagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCt                         <  2:4100538/93‑1 (MQ=255)

                                                                 ccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCtt                        <  1:1522476/93‑1 (MQ=255)

                                                                 ccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCtt                        <  2:419710/93‑1 (MQ=255)

                                                                   agacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTg                      <  2:599541/93‑1 (MQ=255)

                                                                    gacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGa                     >  1:4373547/1‑93 (MQ=255)

                                                                    gacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGa                     >  2:5711982/1‑93 (MQ=255)

                                                                    gacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGa                     >  1:947482/1‑93 (MQ=255)

                                                                    gacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGa                     <  2:3730189/93‑1 (MQ=255)

                                                                     acTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAg                    <  1:6607273/93‑1 (MQ=255)

                                                                     acTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAg                    <  1:5888267/93‑1 (MQ=255)

                                                                         gCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATc                <  2:6516461/93‑1 (MQ=255)

                                                                         gCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATc                <  2:186986/93‑1 (MQ=255)

                                                                          ccAGCAAAAAGGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCa               >  1:4071661/1‑93 (MQ=255)

                                                                          ccAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCa               >  2:6452930/1‑93 (MQ=255)

                                                                           cAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAg              >  1:4480786/1‑93 (MQ=255)

                                                                           cAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAg              <  1:2516096/93‑1 (MQ=255)

                                                                            aGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGc             >  1:3560735/1‑93 (MQ=255)

                                                                             gCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGcc            >  2:1497998/1‑93 (MQ=255)

                                                                              cAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCa           <  2:4671091/93‑1 (MQ=255)

                                                                               aaaaaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAg          >  1:2620176/1‑93 (MQ=255)

                                                                                aaaaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAgg         >  1:3780167/1‑93 (MQ=255)

                                                                                aaaaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAgg         >  1:549779/1‑93 (MQ=255)

                                                                                 aaaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGga        <  1:522113/93‑1 (MQ=255)

                                                                                  aaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGgag       >  1:5096181/1‑93 (MQ=255)

                                                                                   aagcaagcGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGgaga      >  2:2364986/1‑93 (MQ=255)

                                                                                    aGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGaa     <  1:4208935/93‑1 (MQ=255)

                                                                                    aGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGaa     >  1:3216626/1‑93 (MQ=255)

                                                                                    aGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGaa     >  1:1306638/1‑93 (MQ=255)

                                                                                    aGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGaa     <  2:5094712/93‑1 (MQ=255)

                                                                                     gCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAAc    >  1:2676911/1‑93 (MQ=255)

                                                                                     gCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAAc    >  1:5766404/1‑93 (MQ=255)

                                                                                     gCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTCCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAAc    <  1:2050383/93‑1 (MQ=255)

                                                                                      cAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACa   <  2:4335833/93‑1 (MQ=255)

                                                                                      cAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACa   >  2:1436572/1‑93 (MQ=255)

                                                                                       aaGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGCACAg  >  2:4743270/1‑93 (MQ=255)

                                                                                       aaGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg  >  1:3922390/1‑93 (MQ=255)

                                                                                       aaGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg  >  1:298088/1‑93 (MQ=255)

                                                                                       aaGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg  <  1:4895820/93‑1 (MQ=255)

ACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAG  >  NC_000913/3111950‑3112129

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