BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,573,188	G→A	6.3%	pseudogene (775/1503 nt)	yjiT →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,188	0	G	A	6.3%	453.7 / 11.4	142	pseudogene (775/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand): ref base G (68/64); new base A (5/4); total (74/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

Only 100 of 144 total aligned reads displayed.

TAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCTTT  >  NC_000913/4573042‑4573337

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)

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

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

        actCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAg                                                                                                                                            >  2:1004103/1‑150 (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:1092509/1‑150 (MQ=255)

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

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

            acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgcata                                                                                                                                        <  1:174451/150‑8 (MQ=255)

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

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

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

             caACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTGCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                                                                       <  2:1636033/150‑1 (MQ=255)

                ccTAATAAATCACAAAGAACCAGTACAACAGCTTGTGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCATTACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACt                                                                                                                                    <  2:1329126/150‑1 (MQ=255)

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

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

                    ataaATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaa                                                                                                                                <  2:270864/150‑1 (MQ=255)

                       aaTCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaaa                                                                                                                             >  1:1548137/1‑150 (MQ=255)

                          cacaAAGAACCAGTACAAGAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAAt                                                                                                                          <  2:406121/150‑1 (MQ=255)

                          cacaAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAAt                                                                                                                          >  1:564453/1‑150 (MQ=255)

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

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

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

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

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

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

                                     aGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcgcg                                                                                                               <  1:666706/150‑1 (MQ=255)

                                        acaacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg                                                                                                            >  2:1563588/1‑150 (MQ=255)

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

                                          aacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                                                                           <  2:1381264/149‑1 (MQ=255)

                                            caGCTTTAGAAGTTTCACCCTACCTGGCGAAGAGAGTTTCCAATACAGTTAGATGATGAAACAGTAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                                                                        <  2:340423/150‑1 (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:1064512/1‑150 (MQ=255)

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

                                                       ggTTCACCCTACCTGGCGAAGCGAGTATCCAATACAGTTAGATGATCAAACAGGAACGCATTTTTTGAATGGACTGCTATTCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTATATTGTTc                                                                                             >  2:1512795/1‑150 (MQ=255)

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

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

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

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

                                                              ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAg                                                                                      <  1:576669/150‑1 (MQ=255)

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

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

                                                                     ggCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAg                                                                               >  2:1016594/1‑150 (MQ=255)

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

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

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

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

                                                                             gAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCtt                                                                       <  1:1500362/150‑1 (MQ=255)

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

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

                                                                               gTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCATCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATCAGCTTat                                                                     >  2:338911/1‑149 (MQ=255)

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

                                                                                   ccAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTg                                                                 >  2:768053/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                  gttgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCg                                                  <  1:726038/148‑1 (MQ=255)

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

                                                                                                      aaaCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACg                                              >  2:1211369/1‑150 (MQ=255)

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

                                                                                                              aCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcc                                      <  1:811578/150‑1 (MQ=255)

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

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

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

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

                                                                                                                     ttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTg                               <  1:420516/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                  gCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTTTTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTGCCTGTCCTGTAAGAAGCGAGcc                  >  1:300770/1‑150 (MQ=255)

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

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

                                                                                                                                    tATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGt                <  2:1409812/150‑1 (MQ=255)

                                                                                                                                     aggcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTc               >  2:357088/7‑150 (MQ=255)

                                                                                                                                      tGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCa              >  2:693070/1‑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)

                                                                                                                                         gcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACt           >  2:1543545/1‑150 (MQ=255)

                                                                                                                                           aGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTAc         <  1:1175236/150‑1 (MQ=255)

                                                                                                                                            gCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTAcg        <  2:1242290/150‑1 (MQ=255)

                                                                                                                                             ccTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTAcgc       >  2:1434339/1‑150 (MQ=255)

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

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

TAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCTTT  >  NC_000913/4573042‑4573337

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