BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	161,942	T→C	5.9%	P30P (CCA→CCG)	ligT ←	2'‑5' RNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	161,942	0	T	C	5.9%	300.2 / 6.4	99	P30P (CCA→CCG)	ligT	2'‑5' RNA ligase
Reads supporting (aligned to +/- strand): ref base T (47/45); new base C (3/3); total (51/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

Only 100 of 102 total aligned reads displayed.

CGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATGTTAACC  >  NC_000913/161794‑162090

cGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGt                                                                                                                                                      >  4:707974/1‑149 (MQ=255)

cGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGt                                                                                                                                                      <  2:237450/149‑1 (MQ=255)

  gCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTgg                                                                                                                                                    >  3:382963/1‑149 (MQ=255)

   cGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTggg                                                                                                                                                   >  3:477417/1‑149 (MQ=255)

   cGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTggg                                                                                                                                                   <  1:504526/149‑1 (MQ=255)

     gcatcaAGGGTGACGGTGAATCAAGGTGGAAGAACTCGCCCGGCTAAAAGAGAAAGCCCCCTGCTTTTATTGGCGGTCGCTGAGCTTAAAAATGCCCGAGTCAGATGCACATTCTCGGCGAAGACCCGAAGTTCCGCTTCAGGTGTGaa                                                                                                                                                 <  2:354410/148‑1 (MQ=255)

      catcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGacg                                                                                                                                                <  1:318046/149‑3 (MQ=255)

       atcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgt                                                                                                                                               >  3:567591/1‑149 (MQ=255)

       atcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgt                                                                                                                                               >  3:510265/1‑149 (MQ=255)

         caAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgt                                                                                                                                             >  3:743656/1‑149 (MQ=255)

         caAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgt                                                                                                                                             <  1:433150/149‑1 (MQ=255)

           aGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgtgt                                                                                                                                           <  1:787522/149‑1 (MQ=255)

           aGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgtgt                                                                                                                                           <  1:507602/149‑1 (MQ=255)

            gCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGtgtgtg                                                                                                                                          >  1:184385/1‑149 (MQ=255)

               tGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGa                                                                                                                                                  <  4:7153/138‑1 (MQ=255)

               tGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGa                                                                                                                                                  >  3:7153/1‑138 (MQ=255)

                   tgtgAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGcc                                                                                                                                   <  4:510265/149‑1 (MQ=255)

                     tgAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCaa                                                                                                                                 >  1:59428/1‑149 (MQ=255)

                      gAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAAt                                                                                                                                >  4:333889/1‑149 (MQ=255)

                        aaCCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGacgt                                                                                                                                               >  4:247396/1‑129 (MQ=255)

                        aaCCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGACgt                                                                                                                                               <  3:247396/132‑1 (MQ=255)

                          ccAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGAt                                                                                                                            >  1:494851/1‑149 (MQ=255)

                          ccAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGAt                                                                                                                            <  3:560214/149‑1 (MQ=255)

                           cAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATa                                                                                                                           <  3:707974/149‑1 (MQ=255)

                            aGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATaa                                                                                                                          <  4:567591/149‑1 (MQ=255)

                             ggTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAAt                                                                                                                         <  1:625231/149‑1 (MQ=255)

                                 gACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCTCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGt                                                                                                                     >  1:654257/1‑149 (MQ=255)

                                   cGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTc                                                                                                                   <  3:68742/149‑1 (MQ=255)

                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAAt                                                                                                                                >  3:246425/1‑129 (MQ=255)

                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAAt                                                                                                                                <  4:246425/129‑1 (MQ=255)

                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAtt                                                                                                            >  2:249428/1‑149 (MQ=255)

                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAtt                                                                                                            >  4:420590/1‑149 (MQ=255)

                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGACGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAtt                                                                                                            >  2:249767/1‑149 (MQ=255)

                                             cgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCt                                                                                                         <  1:310504/149‑1 (MQ=255)

                                                cTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCa                                                                                                      <  2:157086/149‑1 (MQ=255)

                                                    aaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTa                                                                                                  >  2:268437/1‑149 (MQ=255)

                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                                >  3:36169/1‑149 (MQ=255)

                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                                >  4:787750/1‑149 (MQ=255)

                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                                >  2:344129/1‑149 (MQ=255)

                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGCGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                                >  4:544490/1‑149 (MQ=255)

                                                           aGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGAt                                                                                           <  1:590494/149‑1 (MQ=255)

                                                            gcgcCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATa                                                                                          >  3:527919/1‑149 (MQ=255)

                                                                   ctctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCa                                                                                                      <  2:301016/130‑1 (MQ=255)

                                                                   ctctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCa                                                                                                      >  1:301016/1‑130 (MQ=255)

                                                                    tctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATgg                                                                                                                              >  3:696000/1‑105 (MQ=255)

                                                                    tctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATgg                                                                                                                              <  4:696000/105‑1 (MQ=255)

                                                                    tctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCAGAGACCGCACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCAATTCGCGTATTTATGCATGTAAGTCATTAGCAAAGaa                                                                                  <  2:584481/149‑1 (MQ=255)

                                                                      tttCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACa                                                                                <  3:787750/149‑1 (MQ=255)

                                                                        tctcGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGa                                                                              <  2:59428/149‑1 (MQ=255)

                                                                            ggCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGtt                                                                          <  1:344129/149‑1 (MQ=255)

                                                                                cTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTc                                                                                                                   <  3:9959/104‑1 (MQ=255)

                                                                                cTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCgg                                                                      <  1:684640/149‑1 (MQ=255)

                                                                                   aCTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTc                                                                   >  1:376532/1‑149 (MQ=255)

                                                                                     ttCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAg                                                                 >  1:608972/1‑149 (MQ=255)

                                                                                      tCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGa                                                                <  3:333889/149‑1 (MQ=255)

                                                                                         ccTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAt                                                             >  4:621999/1‑149 (MQ=255)

                                                                                          cTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAtg                                                            >  1:545857/1‑149 (MQ=255)

                                                                                           tAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGAGGGAAGTGCGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAtgt                                                           >  3:201665/1‑149 (MQ=255)

                                                                                              aaaTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAg                                                        <  4:743656/149‑1 (MQ=255)

                                                                                               aaTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAgg                                                       >  2:222614/1‑149 (MQ=255)

                                                                                                aTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGc                                                      >  3:681433/1‑149 (MQ=255)

                                                                                                   ccAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGa                                                                                            <  3:50989/108‑1 (MQ=255)

                                                                                                   ccAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGa                                                                                            >  4:50989/1‑108 (MQ=255)

                                                                                                     agagTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGt                                                                                               >  3:567858/1‑103 (MQ=255)

                                                                                                     agagTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGt                                                                                               <  4:567858/103‑1 (MQ=255)

                                                                                                         tCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTa                                             <  2:641883/149‑1 (MQ=255)

                                                                                                          cAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTAt                                            <  1:268437/149‑1 (MQ=255)

                                                                                                          cAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTAt                                            <  1:25982/149‑1 (MQ=255)

                                                                                                             aTGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGaa                                         >  1:84981/1‑149 (MQ=255)

                                                                                                              tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                        >  4:427685/1‑149 (MQ=255)

                                                                                                              tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                        <  3:420590/149‑1 (MQ=255)

                                                                                                              tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                        >  4:51093/1‑149 (MQ=255)

                                                                                                              tGCAAATTATCGACGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                        <  1:249767/149‑1 (MQ=255)

                                                                                                                 aaaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAc                                     >  2:670714/1‑149 (MQ=255)

                                                                                                                  aaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAcc                                    >  3:88748/1‑149 (MQ=255)

                                                                                                                  aaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAcc                                    >  3:309850/1‑149 (MQ=255)

                                                                                                                    ttATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACcgc                                  >  3:684914/1‑149 (MQ=255)

                                                                                                                    ttATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACcgc                                  <  2:376532/149‑1 (MQ=255)

                                                                                                                       tCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGAt                               <  4:201665/149‑1 (MQ=255)

                                                                                                                                  cGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAAtgtg                    <  1:670714/149‑1 (MQ=255)

                                                                                                                                  cGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAAtgtg                    <  3:621999/149‑1 (MQ=255)

                                                                                                                                   ggACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTTGGTTATGAATTACCGCGATGCTACAATGTgg                   <  2:545857/149‑1 (MQ=255)

                                                                                                                                    gACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGgc                  >  1:453296/1‑149 (MQ=255)

                                                                                                                                       gTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCgg                                                <  2:624150/116‑1 (MQ=255)

                                                                                                                                       gTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCgg                                                >  1:624150/1‑116 (MQ=255)

                                                                                                                                       gTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGcgcg               >  2:209556/1‑149 (MQ=255)

                                                                                                                                         cccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaa             >  4:507464/1‑149 (MQ=255)

                                                                                                                                         cccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaa             >  1:498776/1‑149 (MQ=255)

                                                                                                                                          ccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAg                                                        <  2:496373/105‑1 (MQ=255)

                                                                                                                                          ccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAg                                                        >  1:496373/1‑105 (MQ=255)

                                                                                                                                          ccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaag            >  4:32663/1‑149 (MQ=255)

                                                                                                                                            gCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          <  3:544490/149‑1 (MQ=255)

                                                                                                                                             ccTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAAt         <  4:82966/149‑1 (MQ=255)

                                                                                                                                              cTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATg        <  2:216519/149‑1 (MQ=255)

                                                                                                                                              cTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATg        <  1:249428/149‑1 (MQ=255)

                                                                                                                                                   gTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          >  2:604012/1‑142 (MQ=255)

                                                                                                                                                   gTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          >  2:310483/1‑142 (MQ=255)

                                                                                                                                                   gTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          <  1:604012/142‑1 (MQ=255)

                                                                                                                                                   gTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          <  1:310483/142‑1 (MQ=255)

                                                                                                                                                    tGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATGTTAAcc  >  3:79723/1‑149 (MQ=255)

CGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATGTTAACC  >  NC_000913/161794‑162090

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`