BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,573,189	T→C	6.3%	pseudogene (776/1503 nt)	yjiT →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,189	0	T	C	6.3%	419.6 / 14.1	142	pseudogene (776/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand): ref base T (70/63); new base C (5/4); total (75/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-02

Only 100 of 145 total aligned reads displayed.

TAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCTTTG  >  NC_000913/4573042‑4573338

tAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTag                                                                                                                                                     >  2:1578176/1‑150 (MQ=255)

  gTACTTTCTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTagag                                                                                                                                                   <  1:1017807/150‑1 (MQ=255)

  gTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTagag                                                                                                                                                   <  1:794388/150‑1 (MQ=255)

   tacttactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAgg                                                                                                                                                  <  2:790909/150‑1 (MQ=255)

     cttactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGAGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCt                                                                                                                                                >  2:1356417/1‑150 (MQ=255)

      ttactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTa                                                                                                                                               >  1:957147/1‑150 (MQ=255)

        actCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAg                                                                                                                                             <  1:17286/150‑1 (MQ=255)

        actCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAg                                                                                                                                             >  2:1004103/1‑150 (MQ=255)

        actCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAg                                                                                                                                             <  2:1485258/150‑1 (MQ=255)

         ctCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGc                                                                                                                                            <  1:1612693/150‑1 (MQ=255)

         ctCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgc                                                                                                                                            <  2:671048/150‑5 (MQ=255)

          tCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGcc                                                                                                                                           >  2:329094/1‑150 (MQ=255)

          tCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGcc                                                                                                                                           >  2:35680/1‑150 (MQ=255)

           cacaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCt                                                                                                                                          >  2:1510177/1‑150 (MQ=255)

           cacaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgcat                                                                                                                                          <  1:1571648/150‑7 (MQ=255)

           cacaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgcat                                                                                                                                          <  1:187448/150‑7 (MQ=255)

             caACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                                                                        <  2:357320/150‑1 (MQ=255)

             caACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                                                                        <  1:629813/150‑1 (MQ=255)

             caACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTAACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTTGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                                                                        >  1:479621/1‑150 (MQ=255)

                ccTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACt                                                                                                                                     <  2:1229890/150‑1 (MQ=255)

                 cTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTa                                                                                                                                    <  1:816306/150‑1 (MQ=255)

                 cTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTa                                                                                                                                    >  1:57640/1‑150 (MQ=255)

                           acaAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATa                                                                                                                          <  2:443439/150‑1 (MQ=255)

                            caAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAg                                                                                                                                             >  2:973407/1‑130 (MQ=255)

                            caAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                         >  1:549425/1‑150 (MQ=255)

                             aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTTGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAag                                                                                                                        <  1:1407066/150‑1 (MQ=255)

                               aGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAagag                                                                                                                      <  1:82459/150‑1 (MQ=255)

                                gAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                                                                       <  2:1018494/116‑1 (MQ=255)

                                 aaCCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTa                                                                                                                    <  1:1528750/150‑1 (MQ=255)

                                    cAGTACAACAGCTTGAGTAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcgc                                                                                                                 <  2:1580778/150‑1 (MQ=255)

                                     aGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcgcg                                                                                                                >  1:607294/1‑150 (MQ=255)

                                         caacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                                                                            <  1:1585689/150‑1 (MQ=255)

                                            caGCTTTAGAAGTTTCACCCTACCTGGCGAAGAGAGTTTCCAATACAGTTAGATGATGAAACAGTAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                                                                         <  2:340423/150‑1 (MQ=255)

                                            caGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                                                                         >  2:15218/1‑150 (MQ=255)

                                            caGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                                                                         >  2:1138830/1‑150 (MQ=255)

                                              gCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt                                                                                                       >  1:1409812/1‑150 (MQ=255)

                                              gCTTGAGAAGGTT‑‑CCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCt                                                                                                     >  2:1378005/1‑150 (MQ=255)

                                               cTTTAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAttg                                                                                                      <  1:1047195/150‑2 (MQ=255)

                                                  gagaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAt                                                                                                   >  1:1122207/1‑150 (MQ=255)

                                                  gagaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAt                                                                                                   >  1:1064512/1‑150 (MQ=255)

                                                        gTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCa                                                                                             <  1:1619595/150‑1 (MQ=255)

                                                         ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAg                                                                                            >  1:407212/1‑150 (MQ=255)

                                                         ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAg                                                                                            <  2:1383246/150‑1 (MQ=255)

                                                          tCACCCTACCTGGCGAAGCGTGTTTCCATTACCGTTAGATGATGAATCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGTGTACGGGTTGCCAATTCTATTGGTCAGa                                                                                           <  2:1264053/150‑1 (MQ=255)

                                                             cccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                                                                        >  1:1671035/1‑150 (MQ=255)

                                                              ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCTTTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAg                                                                                       <  2:957147/150‑1 (MQ=255)

                                                                tACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCa                                                                                     <  2:1391703/150‑1 (MQ=255)

                                                                   cTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATcc                                                                                  <  2:613448/150‑1 (MQ=255)

                                                                      gCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGa                                                                               <  1:1004494/150‑1 (MQ=255)

                                                                       cGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGAt                                                                              >  2:899421/1‑150 (MQ=255)

                                                                        gAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATg                                                                             >  1:557177/1‑150 (MQ=255)

                                                                         aaGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGa                                                                            >  2:142847/1‑150 (MQ=255)

                                                                               gTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTTAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAg                                                                      <  1:609157/150‑1 (MQ=255)

                                                                               gTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAg                                                                      <  1:709200/150‑1 (MQ=255)

                                                                                tttCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAgg                                                                     <  1:975630/150‑1 (MQ=255)

                                                                                tatCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAgg                                                                     <  1:1082549/148‑1 (MQ=255)

                                                                                     aaTACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGAt                                                                >  1:43588/1‑150 (MQ=255)

                                                                                      attccGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa                                                               <  2:278600/147‑1 (MQ=255)

                                                                                       tACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGGGATAg                                                              <  1:15218/150‑1 (MQ=255)

                                                                                       tACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTCAGGGTGATAg                                                              <  1:1138830/150‑1 (MQ=255)

                                                                                          cGTTAGATGTTGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTAt                                                           >  2:726038/1‑150 (MQ=255)

                                                                                           gTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTAtc                                                          >  1:1242290/1‑150 (MQ=255)

                                                                                           gTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTAtc                                                          <  1:1578176/150‑1 (MQ=255)

                                                                                               gatgatGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCtt                                                      >  1:346680/1‑150 (MQ=255)

                                                                                                 tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTcc                                                    <  2:1671035/150‑1 (MQ=255)

                                                                                                 tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTcc                                                    <  2:1548137/150‑1 (MQ=255)

                                                                                                 tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTcc                                                    <  2:886207/150‑1 (MQ=255)

                                                                                                 tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTcc                                                    >  1:297358/1‑150 (MQ=255)

                                                                                                   atgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCgg                                                  >  2:1092480/1‑150 (MQ=255)

                                                                                                   atgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCgg                                                  >  2:1719658/1‑150 (MQ=255)

                                                                                                     gaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGAc                                                <  1:933270/150‑1 (MQ=255)

                                                                                                            gAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGttt                                         >  1:971202/1‑150 (MQ=255)

                                                                                                              aCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcc                                       >  2:614687/1‑150 (MQ=255)

                                                                                                               cGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcct                                      >  1:1474936/1‑150 (MQ=255)

                                                                                                                 cATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCACTTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGGATCTCTTCCGGACGAAGtttccgtt                                    >  2:1101250/1‑150 (MQ=255)

                                                                                                                 cATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttccttt                                    >  2:1523411/1‑150 (MQ=255)

                                                                                                                 cATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttccttt                                    >  2:1552289/1‑150 (MQ=255)

                                                                                                                 cATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttccttt                                    >  2:1255298/1‑150 (MQ=255)

                                                                                                                       ttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTa                              <  1:1181684/150‑1 (MQ=255)

                                                                                                                       ttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTa                              <  2:57640/150‑1 (MQ=255)

                                                                                                                          aaTGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACa                           <  1:35680/150‑1 (MQ=255)

                                                                                                                          aaTGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACa                           <  1:329094/150‑1 (MQ=255)

                                                                                                                            tGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAg                         <  2:376803/150‑1 (MQ=255)

                                                                                                                             ggACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGc                        >  1:1019661/1‑150 (MQ=255)

                                                                                                                             ggACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCAT‑CAGATGAGCTTAGGGTGATAGTATCTCTAACGGACTAAGTTTCCTTTCCTGTAACAAGcc                       >  2:743357/1‑149 (MQ=255)

                                                                                                                              gACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCg                       <  2:1233105/150‑1 (MQ=255)

                                                                                                                                cTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAg                     <  2:1474936/150‑1 (MQ=255)

                                                                                                                                cTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAg                     >  2:279386/1‑150 (MQ=255)

                                                                                                                                  gCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGcc                   >  1:1303332/1‑150 (MQ=255)

                                                                                                                                  gCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGcc                   >  2:1420280/1‑150 (MQ=255)

                                                                                                                                    gaggcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGt                 >  1:783775/8‑150 (MQ=255)

                                                                                                                                       gcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCaa              >  1:681700/5‑150 (MQ=255)

                                                                                                                                       gcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCaa              >  1:5678/5‑150 (MQ=255)

                                                                                                                                       gcgcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCaa              >  1:1193656/1‑150 (MQ=255)

                                                                                                                                        cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAAc             >  1:1437071/4‑150 (MQ=255)

                                                                                                                                        cgcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAAc             >  2:1327635/1‑150 (MQ=255)

                                                                                                                                        cgcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAAc             >  1:1195547/1‑150 (MQ=255)

                                                                                                                                                  gTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCttt   >  1:398016/1‑150 (MQ=255)

                                                                                                                                                  gTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCttt   >  1:1005441/1‑150 (MQ=255)

                                                                                                                                                   tAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCTTTg  <  2:785520/150‑1 (MQ=255)

TAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCTTTG  >  NC_000913/4573042‑4573338

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