BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,175,789	T→C	5.8%	noncoding (36/76 nt)	thrT →	tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,175,789	0	T	C	5.8%	503.8 / 15.9	173	noncoding (36/76 nt)	thrT	tRNA‑Thr
Reads supporting (aligned to +/- strand): ref base T (76/85); new base C (6/4); total (84/89)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 179 total aligned reads displayed.

CAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGAC  >  NC_000913/4175640‑4175938

cAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTTGAGCGCACCCTTGGt                                                                                                                                                        >  1:1102438/1‑150 (MQ=255)

    aTGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAgg                                                                                                                                                    <  2:409320/150‑1 (MQ=255)

     tGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAggg                                                                                                                                                   >  1:1514583/1‑150 (MQ=255)

         aaTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAg                                                                                                                                               <  2:1018146/150‑1 (MQ=255)

          aTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTACCAAGGGTgcgc                                                                                                                                              <  2:226397/150‑4 (MQ=255)

           tCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGt                                                                                                                                             <  2:331737/150‑1 (MQ=255)

           tCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGt                                                                                                                                             <  2:228052/150‑1 (MQ=255)

             cggTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGGTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTACCAAGGGTgcgctct                                                                                                                                           <  1:1680596/149‑7 (MQ=255)

               gTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGc                                                                                                                                         <  1:281326/150‑1 (MQ=255)

               gTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGc                                                                                                                                         >  1:549542/1‑150 (MQ=255)

                tAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCa                                                                                                                                        <  1:1468814/150‑1 (MQ=255)

                  gCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAgg                                                                                                                                              >  2:938266/1‑142 (MQ=255)

                  gCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAgg                                                                                                                                              <  1:938266/142‑1 (MQ=255)

                  gCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑t                                                                                                                                     >  1:606710/1‑150 (MQ=255)

                  gCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑t                                                                                                                                     >  1:425048/1‑150 (MQ=255)

                   ccGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑tt                                                                                                                                    >  2:1335168/1‑150 (MQ=255)

                    cGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTc                                                                                                                                   <  2:1423588/150‑1 (MQ=255)

                    cGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTc                                                                                                                                   >  1:520465/1‑150 (MQ=255)

                     gAGTTCCAGGATGCGGGCTTCGTATAATGTCTATTTCCTTAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCGATGATGTGCGGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCg                                                                                                                                  <  2:1572585/150‑1 (MQ=38)

                     gAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCg                                                                                                                                  >  2:1483604/1‑150 (MQ=255)

                       gTTCCGGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGaa                                                                                                                                >  1:231478/1‑150 (MQ=255)

                        ttCCAGGATGCGGGCATCGTATAATGGCTATTACATCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAAt                                                                                                                               <  1:322224/150‑1 (MQ=255)

                         tCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATc                                                                                                                              <  2:1398671/150‑1 (MQ=255)

                           cAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTg                                                                                                                            >  2:478953/1‑150 (MQ=255)

                             ggATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGcc                                                                                                                          <  1:1687027/150‑1 (MQ=255)

                             ggATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCACGCTGCCCGCTCCAAGAGGTGCTGATATAGCTCAGTGGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGcc                                                                                                                          <  1:805383/150‑1 (MQ=255)

                               aTGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTa                                                                                                                        >  1:1223162/1‑150 (MQ=255)

                                 gCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGATTCGAA‑CTGccgacc                                                                                                                      <  1:1390074/150‑5 (MQ=38)

                                  cGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCa                                                                                                                     >  2:946326/1‑150 (MQ=255)

                                   gggCATCGTATAATGGCTATTACCTCAGCCTTCGAAGCGGATGATGCGGGCTCGATTCCCGCTGCCCGCTCCAAGTTGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAg                                                                                                                    <  2:1321765/150‑1 (MQ=255)

                                      cATCGTATAATGGCTATTACCTCCTCCTTCCAAGTGGATGAAGCGGGTTCGATTCCCTCTGACCGCTCCAAGATCTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGcac                                                                                                                 <  2:133761/150‑1 (MQ=37)

                                       aTCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGcacc                                                                                                                <  2:1715499/150‑1 (MQ=255)

                                          gTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGcc                                                                                                                          >  2:1608319/1‑137 (MQ=255)

                                          gTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGcc                                                                                                                          <  1:1608319/137‑1 (MQ=255)

                                           tataATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCActt                                                                                                            >  1:906745/1‑150 (MQ=255)

                                             taATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGAGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCActtct                                                                                                          <  1:393734/150‑1 (MQ=255)

                                               aTGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTGCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCttt                                                                                                        <  2:183136/150‑1 (MQ=255)

                                                    tATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTctcc                                                                                                   <  2:757128/150‑1 (MQ=255)

                                                     aTTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTctcct                                                                                                  >  1:494006/1‑150 (MQ=255)

                                                      ttACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTctcctc                                                                                                 <  1:1329763/150‑1 (MQ=255)

                                                       tACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCtcctcc                                                                                                >  1:678326/1‑150 (MQ=255)

                                                           tCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGt                                                                                            >  2:1255150/1‑150 (MQ=255)

                                                           tCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGt                                                                                            <  1:946326/150‑1 (MQ=255)

                                                             aGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGttt                                                                                          >  2:1157005/1‑150 (MQ=255)

                                                             aGCCTTCCAAGCTGATGAAGCGGGGTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGCGCGCTCCCTTGGTAAGGGTGATGTCGTCAT‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGttt                                                                                          <  1:1050105/150‑1 (MQ=255)

                                                               ccTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGttttt                                                                                        >  1:723573/1‑150 (MQ=255)

                                                               ccTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGttttt                                                                                        <  2:366351/150‑1 (MQ=255)

                                                                cTTCCACGCTGATGATGCGGGATCGATTCCCGCTTCCCGCTCCAAGATCTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttttt                                                                                       <  2:1003337/150‑1 (MQ=255)

                                                                          gatgatGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGttt                                                                             >  1:1160449/1‑150 (MQ=255)

                                                                           atgatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTa                                                                            >  1:1038709/1‑150 (MQ=255)

                                                                              atgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTg                                                                         <  1:1012003/150‑1 (MQ=255)

                                                                               tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGc                                                                        >  1:1656957/1‑150 (MQ=255)

                                                                               tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGc                                                                        <  1:30355/150‑1 (MQ=255)

                                                                                  gggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCAtt                                                                     >  1:226219/1‑150 (MQ=255)

                                                                                    gtgcgcttCCCGCTGCCCGCTCCAAGCTGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTca                                                                   <  1:1203367/144‑1 (MQ=255)

                                                                                       cGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaaca                                                                <  2:706631/150‑1 (MQ=255)

                                                                                        gATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaacaa                                                               <  2:231478/150‑1 (MQ=255)

                                                                                                cTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCa                                                       <  2:451531/150‑1 (MQ=255)

                                                                                                  gCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATg                                                     >  2:640727/1‑150 (MQ=255)

                                                                                                  gCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATg                                                     <  2:499207/150‑1 (MQ=255)

                                                                                                    ccGCTCCAAGATGTGCTGATATAGCGCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGtt                                                   <  2:720074/150‑1 (MQ=255)

                                                                                                      gCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGc                                                 <  2:906745/150‑1 (MQ=255)

                                                                                                       cTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTc                                                                    <  1:569807/130‑1 (MQ=255)

                                                                                                       cTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTc                                                                    >  2:569807/1‑130 (MQ=255)

                                                                                                         ccAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTg                                              >  2:1475557/1‑150 (MQ=255)

                                                                                                           aaGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGt                                            <  1:1280258/150‑1 (MQ=255)

                                                                                                            aGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGtt                                           <  1:353414/150‑1 (MQ=255)

                                                                                                            aGATGTGCTGATATAGCTCAGTAGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTGTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGtt                                           >  2:71194/1‑150 (MQ=255)

                                                                                                               tgtgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTGCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAt                                        >  1:245718/1‑150 (MQ=255)

                                                                                                               tgtgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCGTGGTTGAt                                        >  1:685234/1‑150 (MQ=255)

                                                                                                                 tgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAtgt                                      >  2:1518654/1‑150 (MQ=255)

                                                                                                                 tgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAtgt                                      <  2:857875/150‑1 (MQ=255)

                                                                                                                 tgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAtgt                                      >  1:191451/1‑150 (MQ=255)

                                                                                                                   cTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtgg                                    >  2:1586891/1‑150 (MQ=255)

                                                                                                                   cTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtgg                                    >  2:43173/1‑150 (MQ=255)

                                                                                                                          aGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATc                             <  1:685615/150‑1 (MQ=255)

                                                                                                                           gCTCAGTTGGTAGAGCGCACGCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAACTCGGGCATGTTGCCTGGTTGATGTGGTGATATcc                            >  1:1514835/1‑149 (MQ=255)

                                                                                                                           gCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCa                            >  1:1563250/1‑150 (MQ=255)

                                                                                                                           gCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TCCGAATCTGCCTATCAGCACCACTTCTTTTCTCCCCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCa                            >  2:522968/1‑150 (MQ=255)

                                                                                                                            cTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCAc                           >  1:65474/1‑150 (MQ=255)

                                                                                                                            cTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCAc                           <  2:593152/150‑1 (MQ=255)

                                                                                                                            cTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCAc                           >  2:685436/1‑150 (MQ=255)

                                                                                                                              cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCg                         >  1:1425765/1‑150 (MQ=255)

                                                                                                                               aGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGa                        <  2:1102438/150‑1 (MQ=255)

                                                                                                                               aGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGGCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGa                        <  2:1214717/150‑1 (MQ=255)

                                                                                                                                    gTAGAGCGCACCCTTTGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTATTCCTTCTGTGTATTGCCTTCAACAAGTCGGTCAGTTTGCCAGGTTGATGTGGTGATATAACAGATTTAt                   >  2:315531/1‑150 (MQ=255)

                                                                                                                                    gTAGAGCGCACCCTTTGTAAGGGTGAGGTCGGCAG‑TTCGAAACTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTAt                   >  2:853719/1‑150 (MQ=255)

                                                                                                                                       gacctcACCCTTACCAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCg                >  2:1621808/6‑150 (MQ=255)

                                                                                                                                       gacctcACCCTTACAAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCg                >  2:530808/6‑150 (MQ=255)

                                                                                                                                        acctcACCCTTACCAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCgt               >  1:954838/5‑150 (MQ=255)

                                                                                                                                         gcgcACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCgtg              <  2:606710/150‑1 (MQ=255)

                                                                                                                                         gcgcACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCgtg              <  1:704419/150‑1 (MQ=255)

                                                                                                                                         cctcACCCTTACCAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATACgtg              >  1:1340906/4‑150 (MQ=255)

                                                                                                                                            cacaCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGAGTCt           >  2:747931/1‑150 (MQ=255)

                                                                                                                                            cACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCt           >  2:395046/1‑150 (MQ=255)

                                                                                                                                             aCCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCtt          >  1:787025/1‑150 (MQ=255)

                                                                                                                                                cTTGGTAAGGGTGAGGTCGGCAG‑TTCGCATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCGGGTTGATGTGGGGATATCACCGATTTATCCGTGTCTTaga       >  2:1065824/1‑150 (MQ=255)

                                                                                                                                                   ggTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCTGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAggg    <  2:281346/150‑1 (MQ=255)

                                                                                                                                                     tAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGAc  <  1:1475557/150‑1 (MQ=255)

                                                                                                                                                     tAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGAc  <  1:877980/150‑1 (MQ=255)

CAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGAC  >  NC_000913/4175640‑4175938

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