BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,175,788	G→C	6.1%	noncoding (35/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,788	0	G	C	6.1%	586.4 / 20.4	175	noncoding (35/76 nt)	thrT	tRNA‑Thr
Reads supporting (aligned to +/- strand): ref base G (80/84); new base C (7/4); total (87/88)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 178 total aligned reads displayed.

GCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGG  >  NC_000913/4175639‑4175936

 cAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGCTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGt                                                                                                                                                      >  1:227860/1‑150 (MQ=255)

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

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

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

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

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

           aTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAgg                                                                                                                                            >  1:706631/1‑150 (MQ=255)

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

            tCAGGTAGCCGCGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTACCAAGGGTgcgct                                                                                                                                           <  2:139079/150‑5 (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:617324/150‑1 (MQ=255)

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

                gTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCATCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCTGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGc                                                                                                                                       <  1:1676261/150‑1 (MQ=255)

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

                 tAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCa                                                                                                                                      >  1:451531/1‑150 (MQ=255)

                  aGCCGAGGTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAg                                                                                                                                     <  1:1578936/150‑1 (MQ=255)

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

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

                     cGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGTTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTc                                                                                                                                 <  2:1716434/150‑1 (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)

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

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

                              ggATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGcc                                                                                                                        <  2:597492/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)

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

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

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

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

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

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

                                                      aTTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTctcct                                                                                                <  1:1650326/150‑1 (MQ=255)

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

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

                                                         aCCTCAGCCTTCCAAGCTGATGAGGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATTTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTccc                                                                                             <  1:336429/150‑1 (MQ=255)

                                                          ccTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCt                                                                                            >  1:688243/1‑150 (MQ=255)

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

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

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

                                                               gCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCATGtttt                                                                                       >  2:1316862/1‑150 (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)

                                                                    ccAAGCTGATGATGCGGGTTCGATTCCCGCCGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTGGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCt                                                                                  >  1:281346/1‑150 (MQ=255)

                                                                     cAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCAGCCTATCAGCACCACTTCTTTTCTCCTCCCTGATTTTTCCtt                                                                                 >  1:578760/1‑150 (MQ=255)

                                                                      aaGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTACAAGATGTGCTGATATAGCTCAGTTGGTAGAGCTCACCCTTGGTAAGGGTGAGGTCGGCTG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                                                                                >  1:203015/1‑150 (MQ=255)

                                                                        gCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTg                                                                              >  2:1418721/1‑150 (MQ=255)

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

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

                                                                              gatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTAtt                                                                        >  1:1702155/1‑150 (MQ=255)

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

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

                                                                                     gTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTca                                                                 <  2:1530511/150‑1 (MQ=255)

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

                                                                                            tCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTc                                                          <  1:841390/150‑1 (MQ=255)

                                                                                               cGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCggg                                                       <  1:1524803/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                       gCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTATTTTCTCCTCCCTGTTTTTTCCTTCTGTGTTTTGCATTCAACAAGTCGGGCATGTTgg                                               >  1:269500/1‑149 (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                                         <  2:483370/150‑1 (MQ=255)

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

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

                                                                                                               aTGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGa                                       <  1:1485336/150‑1 (MQ=255)

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

                                                                                                                  ttcTGATATAGTGCCGGTGGTAGAGCGTACCCTTGGTAAGGGGGAGGTCGGCAG‑GTCGAATCTGCCTATCCGCACCACTTCTTGTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGGTGAtgt                                    <  1:320552/148‑1 (MQ=255)

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

                                                                                                                   gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAtgtg                                   >  2:266662/1‑150 (MQ=255)

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

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

                                                                                                                      gATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtggtg                                >  1:519619/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                 gTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAt                     <  2:1656957/150‑1 (MQ=255)

                                                                                                                                  ttGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGATTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAtt                    <  2:1514583/150‑1 (MQ=255)

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

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

                                                                                                                                       agagCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATcc               <  2:1323908/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                                             cacaCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGAGTCt         >  2:747931/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)

GCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGG  >  NC_000913/4175639‑4175936

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: `‑`