BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,113,474	T→A	9.3%	L61* (TTG→TAG) ‡	uspD →	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,474	0	T	A	9.3%	318.3 / 15.2	117	L61* (TTG→TAG) ‡	uspD	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base T (58/48); new base A (3/8); total (61/56)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

CGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGT  >  NC_000913/4113325‑4113622

cggttaTGAAGAAGATGCCTTACTGGTGAATTAAGCCATGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAtt                                                                                                                                                      <  2:292955/145‑1 (MQ=255)

cGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAtt                                                                                                                                                      <  2:1063034/150‑1 (MQ=255)

cGGGAATGAAGAAGATGCCTTACTGGTGAACAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGGTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAtt                                                                                                                                                      <  1:665769/150‑1 (MQ=255)

 gggAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                                                                     <  1:1501871/150‑1 (MQ=255)

 gggAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                                                                     >  1:615912/1‑150 (MQ=255)

  ggAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                                                                    <  2:634144/150‑1 (MQ=255)

  ggAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                                                                    <  1:381307/150‑1 (MQ=255)

  ggAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                                                                    <  1:1110746/150‑1 (MQ=255)

  ggAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTAACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                                                                    <  1:100687/150‑1 (MQ=255)

   gAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGAGAGTTGTATCCGGGTATCTACTTCCCTTCATCAGAAGATATTCTTCAAttgtt                                                                                                                                                   >  2:1207999/1‑150 (MQ=255)

    aaTGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCTCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                                                                  <  1:382557/150‑1 (MQ=255)

       taagaagATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                                                                               <  1:591451/149‑1 (MQ=255)

       gaagaaGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                                                                               >  1:549161/1‑150 (MQ=255)

           aagaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                                                                           <  1:776846/150‑1 (MQ=255)

            agaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATaa                                                                                                                                          <  1:1657786/150‑1 (MQ=255)

            acatcccTTACTGTTGAATAAAGAAATGTAGCTCGCCAGACATAATGACGCTCACCTTAAGTTAATTCATATTGATTATTGCTTAAGGGAGTTGTATCCGGGTATCTACTTCGCTGCAACAGAAGATATTCTTCAACAATTGAAGAATat                                                                                                                                          <  2:267805/145‑2 (MQ=255)

              aTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatct                                                                                                                                        <  1:756359/150‑4 (MQ=255)

              aTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatct                                                                                                                                        <  2:524456/150‑4 (MQ=255)

                gCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg                                                                                                                                      >  2:162633/1‑150 (MQ=255)

                cccTAACTGCTGAATAAAGCCCAGGAGCTCGCCAGACATTATCACCCTCACCTGACGATATTTCATATTGATGAAGGCTTATGCGAGTTGTATCCGGCTTTCTACTTCCCTGCAACAGAAGATATTCTTCTATTGTTGAAGAATAAGTCg                                                                                                                                      <  1:837719/149‑1 (MQ=255)

                 ccTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                                                                     <  2:1295760/150‑7 (MQ=255)

                 ccTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACCCTCACCTGACGTTAATTCATTTTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGTAGAATAAGTCgg                                                                                                                                     >  2:106962/1‑150 (MQ=255)

                  cTTACTGGTTAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctg                                                                                                                                    <  2:341015/150‑8 (MQ=255)

                       tGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                                                                               >  1:1626365/1‑150 (MQ=255)

                        ggTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTACAATTGTTGAAGAATAAGTCGGATAACaa                                                                                                                              >  2:259798/1‑150 (MQ=255)

                        ggTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCGTCAATTGTTGAAGAATAAGTCGGATAACaa                                                                                                                              >  1:1497611/1‑150 (MQ=255)

                        ggTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGAATCTACTTCCATGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACaa                                                                                                                              >  1:1605952/1‑150 (MQ=255)

                           gAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCt                                                                                                                           >  2:712535/1‑150 (MQ=255)

                            aaTAAAGCCCTGGAGCTCGCCAGACATAATGACTCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                                                                          <  1:584755/150‑1 (MQ=255)

                             aTAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGt                                                                                                                         >  1:203279/1‑150 (MQ=255)

                                aaGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                                                                      <  1:62361/150‑1 (MQ=255)

                                   cccTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                                                                   >  2:1375172/1‑150 (MQ=255)

                                      tgtaGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGAAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGa                                                                                                                <  1:299716/147‑1 (MQ=255)

                                       ttagttCGCCAGAGATAATGACGCTCTCCTGACGTTAATTAATATTGATGATGGCTTAAGGGAGTTGTATCCGGGTATCTACTTCCCTGCAACATAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                                                                               <  2:1486014/145‑1 (MQ=255)

                                        gAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACg                                                                                                              >  1:1563175/1‑150 (MQ=255)

                                             cGCCAGACTTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaaa                                                                                                         >  1:1685894/1‑150 (MQ=255)

                                             cGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCTACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaaa                                                                                                         >  1:1601427/1‑150 (MQ=255)

                                              gCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAat                                                                                                        >  1:750311/1‑150 (MQ=255)

                                                  gACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                                                                    >  1:1157360/1‑150 (MQ=255)

                                                    cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                                                                                  >  1:1005475/1‑150 (MQ=255)

                                                     aTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAAt                                                                                                 <  1:879982/150‑1 (MQ=255)

                                                      tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                                                                                >  2:611464/1‑150 (MQ=255)

                                                      tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                                                                                >  1:260372/1‑150 (MQ=255)

                                                          gACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                                                                            <  2:615912/150‑1 (MQ=255)

                                                            cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGaa                                                                                          >  2:1258893/1‑150 (MQ=255)

                                                            cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGaa                                                                                          >  1:1264888/1‑150 (MQ=255)

                                                             gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                                                                         >  2:1018240/1‑150 (MQ=255)

                                                                cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATACTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACa                                                                                      <  1:277752/150‑1 (MQ=255)

                                                                cACCTGACGTTAATTCATATTGATGATGGCTTAAGAGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTATTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACa                                                                                      >  2:401556/1‑150 (MQ=255)

                                                                 aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                                                                   <  2:1682610/136‑1 (MQ=255)

                                                                 aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                                                                   >  1:1682610/1‑136 (MQ=255)

                                                                  ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaa                                                                                    <  2:756699/150‑1 (MQ=255)

                                                                    tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                                                                  >  2:1409988/1‑150 (MQ=255)

                                                                     gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACt                                                                                 >  1:569480/1‑150 (MQ=255)

                                                                       cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                                                                               >  2:279434/1‑150 (MQ=255)

                                                                         ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGc                                                                                              >  2:139626/1‑133 (MQ=255)

                                                                           aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                                                                           >  2:389532/1‑150 (MQ=255)

                                                                           aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                                                                           <  1:1018240/150‑1 (MQ=255)

                                                                                  atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                                                                    <  2:203279/150‑1 (MQ=255)

                                                                                    tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                                                                  <  2:549161/150‑1 (MQ=255)

                                                                                      atgatgGCTTAAGCGAGTTGTATCCTGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAATCTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag                                                                <  1:1052402/150‑1 (MQ=255)

                                                                                        gatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa                                                              >  2:1604575/1‑150 (MQ=255)

                                                                                             cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc                                                         <  2:260372/150‑1 (MQ=255)

                                                                                             cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc                                                         <  2:1626365/150‑1 (MQ=255)

                                                                                              ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg                                                        <  2:580680/150‑1 (MQ=255)

                                                                                              ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAAGCGGAGAAATGCCg                                                        >  2:164282/1‑150 (MQ=255)

                                                                                                aaGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa                                                      <  2:1001314/150‑1 (MQ=255)

                                                                                                  gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaaa                                                    >  1:1014266/1‑150 (MQ=255)

                                                                                                  gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCAGAAGACAAAACTGCGTATTCAACGCGGAGAAATGCCGGaaa                                                    >  2:549641/1‑150 (MQ=255)

                                                                                                   cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGACGAAGAAAAAACTGCGTATTGAACGAGTAGAAATGCCGGAAac                                                   >  2:1101662/1‑150 (MQ=255)

                                                                                                     aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAGGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAacac                                                 <  1:1174708/150‑1 (MQ=255)

                                                                                                      gTTGTATCCGGGTATCTACTTCCCTGCTACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACAct                                                <  2:1601427/150‑1 (MQ=255)

                                                                                                      gTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                                                                    <  1:660265/130‑1 (MQ=255)

                                                                                                          tATCCGGGTATCTACTTCCCGGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATACCAAGCTGTATAAACTGACGAAAAATATTCAAAGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACActgctg                                            >  1:1105422/1‑150 (MQ=255)

                                                                                                           aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGTATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTgg                                           <  1:1254543/150‑1 (MQ=255)

                                                                                                            tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGAAAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGa                                          <  2:805059/150‑1 (MQ=255)

                                                                                                              cGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGaaa                                        <  1:1375172/150‑1 (MQ=255)

                                                                                                                 gTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTa                                     <  2:1264888/150‑1 (MQ=255)

                                                                                                                  tATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCGGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTAt                                    <  1:1409988/150‑1 (MQ=255)

                                                                                                                      tACTTCACTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTAAATTGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGTCGGAAACACTGCTGGAAATTATGCaa                                >  2:1203842/1‑150 (MQ=255)

                                                                                                                          tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATTCCGTAAACACTGCTGGAAATTATGCaaaaaa                            >  2:310909/1‑149 (MQ=255)

                                                                                                                          tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAag                            >  1:296185/1‑150 (MQ=255)

                                                                                                                             cTGCTACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGc                         >  1:601021/1‑150 (MQ=255)

                                                                                                                             cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACTAAAAATATTCAATGGCCGAAGACAAAACTGCGTCTTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGc                         <  1:1008404/150‑1 (MQ=255)

                                                                                                                             cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGc                         >  1:358629/1‑150 (MQ=255)

                                                                                                                              tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCa                        >  1:1059671/1‑150 (MQ=255)

                                                                                                                                   cAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCg                   >  2:1299118/1‑150 (MQ=255)

                                                                                                                                     gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGCCGAAAAATATTCAATGGCAGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCgcc                 >  2:304059/1‑148 (MQ=255)

                                                                                                                                     gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAACCAGTATAAACTGACGAAAAATATTCAATGGACGAAGACAAAACTGCGTATTGAAAGCGGAGAAATGCCGGAAACACTGCTGGAAATTATTCAAAAAGAGCAGTGCGAc                 >  2:132702/1‑150 (MQ=255)

                                                                                                                                      gacttatTCTTCAACAATTGAAGAATAACTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGAcc                >  2:1536713/5‑150 (MQ=255)

                                                                                                                                         atatTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGAcctcc             <  2:1563175/150‑1 (MQ=255)

                                                                                                                                          tatTCTTCAACTGTTGAAGAATAAGGCGGATACCAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTCATGGAAATTATGCAAAAAGAGCCGTGCGAcctcct            >  1:41944/1‑150 (MQ=255)

                                                                                                                                           atTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCtt           <  1:963296/150‑1 (MQ=255)

                                                                                                                                             tcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGt         >  2:1659298/1‑150 (MQ=255)

                                                                                                                                             tcttcAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGt         >  1:1147211/1‑150 (MQ=255)

                                                                                                                                              cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGta        >  2:1529845/1‑149 (MQ=255)

                                                                                                                                              cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTc        <  1:1359871/150‑1 (MQ=255)

                                                                                                                                                  aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCtgtg    <  1:279434/150‑1 (MQ=255)

                                                                                                                                                  aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCtgtg    >  1:1175268/1‑150 (MQ=255)

                                                                                                                                                    ttgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGt  <  1:1604575/150‑1 (MQ=255)

CGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGT  >  NC_000913/4113325‑4113622

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