BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,423,621	C→G	6.9%	noncoding (35/76 nt)	thrV ←	tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,621	0	C	G	6.9%	481.4 / 19.8	144	noncoding (35/76 nt)	thrV	tRNA‑Thr
Reads supporting (aligned to +/- strand): ref base C (70/63); new base G (3/7); total (74/70)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 152 total aligned reads displayed.

GGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGC  >  NC_000913/3423472‑3423770

ggCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTAc                                                                                                                                                       <  2:1642168/150‑1 (MQ=255)

 gcgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCTCCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTAcc                                                                                                                                                      <  2:722815/150‑1 (MQ=255)

 gcgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTAcc                                                                                                                                                      <  2:1557094/150‑1 (MQ=255)

  cgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCa                                                                                                                                                     >  1:2431059/1‑150 (MQ=255)

    cTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAg                                                                                                                                                   <  2:862943/150‑1 (MQ=255)

      aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCTAACTGGGGACCTCACCCTTACCAAggg                                                                                                                                                 <  1:644767/150‑1 (MQ=255)

       cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGt                                                                                                                                                >  1:890360/1‑150 (MQ=255)

       cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGt                                                                                                                                                <  2:585174/150‑1 (MQ=255)

       cGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGt                                                                                                                                                >  1:424733/1‑150 (MQ=255)

        ggCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGCTACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTg                                                                                                                                               <  1:653856/150‑1 (MQ=255)

        ggCGTTTCACTTCTGAGTTCGGCATGGGGGCAGGTGGCACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTg                                                                                                                                               <  2:2488195/150‑1 (MQ=255)

         gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTga                                                                                                                                              <  2:756282/150‑2 (MQ=255)

          ggTTTCACTTCTGAGTTCGGCATGGGGGCAGGTGGGACCACCGCGGTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgag                                                                                                                                             <  2:742563/149‑3 (MQ=37)

          cGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgag                                                                                                                                             <  1:224113/150‑3 (MQ=255)

          cGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgcg                                                                                                                                             <  1:1898778/150‑1 (MQ=255)

           gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgcgc                                                                                                                                            <  2:2010476/150‑1 (MQ=255)

            tttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgaggt                                                                                                                                           <  1:1288208/150‑5 (MQ=255)

              tCACTTCTGAGTTCGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGctct                                                                                                                                         <  1:1379270/150‑1 (MQ=255)

                 cTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAcc                                                                                                                                      <  2:1258499/150‑1 (MQ=255)

                  ttCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCa                                                                                                                                     <  1:863242/150‑1 (MQ=255)

                    cTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAAc                                                                                                                                   >  1:1778072/1‑150 (MQ=255)

                    cTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAAc                                                                                                                                   <  1:352162/150‑1 (MQ=255)

                     tGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACt                                                                                                                                  <  2:1363286/150‑1 (MQ=255)

                        gTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAg                                                                                                                               <  2:939541/150‑1 (MQ=255)

                        gTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCAAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAg                                                                                                                               <  2:2076639/150‑1 (MQ=255)

                                   ggTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCa                                                                                                                    >  2:1462441/1‑150 (MQ=255)

                                   ggTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCa                                                                                                                    <  2:1729077/150‑1 (MQ=255)

                                    gTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCAc                                                                                                                   <  2:1018246/150‑1 (MQ=255)

                                    gTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTAGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCAc                                                                                                                   <  2:1508438/150‑1 (MQ=255)

                                        ggTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTa                                                                                                               >  1:945689/1‑150 (MQ=255)

                                           gggACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAAt                                                                                                            >  1:2045998/1‑150 (MQ=255)

                                             gACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAAttt                                                                                                          <  2:672593/150‑1 (MQ=255)

                                              accaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTg                                                                                                         <  1:1349331/150‑1 (MQ=255)

                                              accaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTg                                                                                                         >  1:1510154/1‑150 (MQ=255)

                                               ccaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGa                                                                                                        >  1:2008312/1‑150 (MQ=255)

                                                 accGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATg                                                                                                      <  2:2228347/150‑1 (MQ=255)

                                                  ccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGc                                                                                                     >  2:251737/1‑150 (MQ=255)

                                                       cTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGc                                                                                                <  1:251737/150‑1 (MQ=255)

                                                       cTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGc                                                                                                <  2:1723961/150‑1 (MQ=255)

                                                        tACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCa                                                                                               >  2:1019675/1‑150 (MQ=255)

                                                         aCTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAg                                                                                              <  2:2379547/150‑1 (MQ=255)

                                                         aCTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAg                                                                                              <  1:1762075/150‑1 (MQ=255)

                                                         aCTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAg                                                                                              >  2:1352637/1‑150 (MQ=255)

                                                              cgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTccc                                                                                         <  2:2410764/150‑1 (MQ=255)

                                                              cgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTccc                                                                                         <  2:2045998/150‑1 (MQ=255)

                                                                 cAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTAc                                                                                      <  1:1062324/150‑1 (MQ=255)

                                                                 cAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTAc                                                                                      >  1:70946/1‑150 (MQ=255)

                                                                 cAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTAc                                                                                      <  1:2177313/150‑1 (MQ=255)

                                                                     cAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACtctc                                                                                  >  2:370308/1‑150 (MQ=255)

                                                                         ttCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCGTATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCAt                                                                              >  2:1005422/1‑150 (MQ=255)

                                                                          tCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                             >  1:857398/1‑150 (MQ=255)

                                                                             tttCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATgggg                                                                          <  2:1970908/150‑1 (MQ=255)

                                                                               tCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGgag                                                                        >  1:1626890/1‑150 (MQ=255)

                                                                                cTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGgaga                                                                       >  1:782336/1‑150 (MQ=255)

                                                                                      tGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGCCCCcac                                                                 >  2:1808542/1‑150 (MQ=255)

                                                                                       gCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                                >  1:1984371/1‑150 (MQ=255)

                                                                                        ccGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGTTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTGTCGCATTGGGAGACCCCacac                                                               <  2:680373/150‑1 (MQ=255)

                                                                                           aaCTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCGTACTCTCGCATGGGGAGACCCCACACTAc                                                            <  2:424733/150‑1 (MQ=255)

                                                                                             cTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCa                                                          >  1:467349/1‑150 (MQ=255)

                                                                                                 aaCCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCgg                                                      <  1:1461329/150‑1 (MQ=255)

                                                                                                 aaCCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCgg                                                      <  2:890360/150‑1 (MQ=255)

                                                                                                    cTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgcgc                                                   >  2:1155625/1‑150 (MQ=255)

                                                                                                    cTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgcgc                                                   >  2:623578/1‑150 (MQ=255)

                                                                                                    cTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgcgc                                                   >  2:606856/1‑150 (MQ=255)

                                                                                                         aaGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACgg                                              >  1:1555940/1‑150 (MQ=255)

                                                                                                                  tGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                     >  2:533572/1‑150 (MQ=255)

                                                                                                                    aTACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                   <  1:1462441/150‑1 (MQ=255)

                                                                                                                     tACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                  >  1:262264/1‑150 (MQ=255)

                                                                                                                      acgcAGAGTCGAACTGGGGACGTCACCCTTACCAAGGGTGCGCTCTAGCAACTGAGCCATATGAGCACGCTATATTTGATGCGTGGCAGTTCCCTACTCTCGCTTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                 <  1:1695822/147‑1 (MQ=25)

                                                                                                                       cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                >  1:625657/1‑150 (MQ=255)

                                                                                                                       cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                >  1:265272/1‑150 (MQ=255)

                                                                                                                       cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                >  1:2437310/1‑150 (MQ=255)

                                                                                                                       cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                >  1:1889786/1‑150 (MQ=255)

                                                                                                                       cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGAACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCCTTTCACTTCTGa                                >  2:2405293/1‑150 (MQ=255)

                                                                                                                        ccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                               >  2:1869984/1‑150 (MQ=255)

                                                                                                                         cAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGt                              <  2:1778072/150‑1 (MQ=255)

                                                                                                                             gTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                          <  2:782336/150‑1 (MQ=255)

                                                                                                                             gTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                          >  1:2036358/1‑150 (MQ=255)

                                                                                                                              tCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                         <  2:1510154/150‑1 (MQ=255)

                                                                                                                               cGAACTGGGGACCTCACCCTTACCAAGTGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                        >  2:2157649/1‑150 (MQ=255)

                                                                                                                               cGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                        >  2:1560571/1‑150 (MQ=255)

                                                                                                                                gAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAt                       >  1:736365/1‑150 (MQ=255)

                                                                                                                                gAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAt                       >  1:1519979/1‑150 (MQ=255)

                                                                                                                                gAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAt                       >  1:1372542/1‑150 (MQ=255)

                                                                                                                                  aCTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                     <  1:1019675/150‑1 (MQ=255)

                                                                                                                                       ggACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCa                >  1:2054428/1‑150 (MQ=255)

                                                                                                                                          gcgcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt             >  1:578782/4‑150 (MQ=255)

                                                                                                                                          gcgcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt             >  2:438494/4‑150 (MQ=255)

                                                                                                                                          ccTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt             >  1:51143/1‑150 (MQ=255)

                                                                                                                                          ccTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt             >  2:1650487/1‑150 (MQ=255)

                                                                                                                                           cTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg            >  2:384804/1‑150 (MQ=255)

                                                                                                                                           cTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg            <  1:2123956/150‑1 (MQ=255)

                                                                                                                                            tCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg           <  2:2008312/150‑1 (MQ=255)

                                                                                                                                             cACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg          >  1:1814579/1‑150 (MQ=255)

                                                                                                                                               cccTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac        >  1:715406/1‑150 (MQ=255)

                                                                                                                                               cccTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac        >  1:293679/1‑150 (MQ=255)

                                                                                                                                                 cTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGacca      >  2:216787/1‑150 (MQ=255)

                                                                                                                                                  ttACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGaccac     >  2:1471398/1‑150 (MQ=255)

                                                                                                                                                    aCCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcg   >  2:2058579/1‑150 (MQ=255)

                                                                                                                                                     ccAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcgc  >  1:2371677/1‑150 (MQ=255)

GGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGC  >  NC_000913/3423472‑3423770

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