BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,648,132	A→C	8.4%	noncoding (70/84 nt)	agrA →	small RNA AgrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,648,132	0	A	C	8.4%	496.7 / 31.1	176	noncoding (70/84 nt)	agrA	small RNA AgrA
Reads supporting (aligned to +/- strand): ref base A (67/92); new base C (4/11); total (72/104)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

TGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAAT  >  NC_000913/3647983‑3648281

tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCaaa                                                                                                                                                       >  2:202807/1‑150 (MQ=255)

 tgATAATCATACAGTACATGCATGCTATATAACCAGCAAGTCCTTGCAATATTTTCAGTATGGTATTAGCATTGATGCGTTAGATTATTTCTATCTCACTCCAGGCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAg                                                                                                                                                      <  2:1241834/149‑1 (MQ=25)

 ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAg                                                                                                                                                      >  1:493129/1‑150 (MQ=255)

     aaTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCgg                                                                                                                                                  <  2:1758132/150‑1 (MQ=255)

        cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttt                                                                                                                                               <  2:2021173/150‑1 (MQ=255)

          tACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGGttttt                                                                                                                                             <  1:2071500/150‑1 (MQ=38)

          tACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttttt                                                                                                                                             >  2:2111923/1‑150 (MQ=255)

           aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGGtttttt                                                                                                                                            <  2:1966728/150‑1 (MQ=38)

           aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGGtttttt                                                                                                                                            <  2:2393302/150‑1 (MQ=38)

           aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                            <  1:2126315/150‑1 (MQ=255)

           aCAGCAGATGCACGCTATAAAAACAGCACCTCCTTGCAATCTTTTCCCCTTGGTAATAGCATTGATGCGTTAGATCATCCCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGGtttttt                                                                                                                                            <  1:2359052/150‑1 (MQ=12)

            cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                           >  2:673162/1‑150 (MQ=255)

            cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                           <  1:316105/150‑1 (MQ=255)

            cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                           <  2:855826/150‑1 (MQ=255)

            cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                            >  1:773938/1‑149 (MQ=255)

             aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                          <  1:248589/150‑3 (MQ=38)

             aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                          <  2:1428655/150‑3 (MQ=38)

             aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                          <  2:1420139/150‑3 (MQ=38)

              gatcttGCATGTTATAAAACCAGGACGTCCTTTCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTAc                                                                                                                                         <  2:791760/145‑1 (MQ=255)

              gTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                          <  2:277775/150‑3 (MQ=38)

              gTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTAc                                                                                                                                         <  2:1332479/150‑1 (MQ=255)

               tACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACa                                                                                                                                        <  2:1582609/150‑1 (MQ=255)

                 cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGGttttttactttc                                                                                                                                      <  1:1075458/150‑7 (MQ=38)

                 cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATc                                                                                                                                      >  2:421847/1‑150 (MQ=255)

                  gtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATcc                                                                                                                                     >  1:333951/2‑150 (MQ=255)

                  gtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATcc                                                                                                                                     >  2:1291244/2‑150 (MQ=255)

                   tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACagatc                                                                                                                                    >  1:791760/1‑146 (MQ=255)

                        gTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTc                                                                                                                               <  1:669139/150‑1 (MQ=255)

                         ttATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCg                                                                                                                              <  1:2472981/150‑1 (MQ=255)

                          tataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCTTAGCCAACAACTCAGGGCTGGAAAGT‑TAAATACCGACGGAAAGTCGGTTTTTTTACATGCGGAtttgt                                                                                                                             <  2:191842/150‑4 (MQ=255)

                          tataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCgg                                                                                                                             <  1:159095/150‑1 (MQ=255)

                          tataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCgg                                                                                                                             >  1:2219359/1‑150 (MQ=255)

                            taAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGAc                                                                                                                           <  1:1752739/150‑1 (MQ=255)

                            aaaaaCCAGCAATTTCTTGCAATTTTTTCCCTTTGGTATTAGTATTGATGTGTTAGATTAGTCCTATCTCACTCCAGGCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAc                                                                                                                           <  2:51354/149‑1 (MQ=11)

                             aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAACT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACa                                                                                                                          <  1:1214046/150‑1 (MQ=255)

                              aaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACaa                                                                                                                         >  2:2220655/1‑150 (MQ=255)

                                     cACGTCCTTTCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGTCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGTAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTaa                                                                                                                  <  2:748242/150‑1 (MQ=255)

                                     cACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTaa                                                                                                                  >  1:692975/1‑150 (MQ=255)

                                     cACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGATTaa                                                                                                                  >  1:582169/1‑150 (MQ=255)

                                         tCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATg                                                                                                              >  1:841045/1‑150 (MQ=255)

                                          ccTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                                                                             >  1:1959673/1‑150 (MQ=255)

                                           cTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                                                                            >  2:1797153/1‑150 (MQ=255)

                                                 aTAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                                                                      <  2:2429618/150‑1 (MQ=255)

                                                 aTAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                                                                      >  2:2474497/1‑150 (MQ=255)

                                                     tttCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                                                                                  <  1:1291244/150‑1 (MQ=255)

                                                      ttCAGTATGGTATAAGCATTGATGCGTTTGATGAGGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGt                                                                                                 <  2:1440388/150‑1 (MQ=255)

                                                        cAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGa                                                                                               >  2:2309317/1‑150 (MQ=255)

                                                        cAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGa                                                                                               >  1:1166229/1‑150 (MQ=255)

                                                          gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaa                                                                                             >  1:1743140/1‑150 (MQ=255)

                                                           tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCgg                                                                                                                                                  <  1:943598/96‑1 (MQ=255)

                                                           tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCgg                                                                                                                                                  >  2:943598/1‑96 (MQ=255)

                                                             tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                          <  2:1829042/150‑1 (MQ=255)

                                                             tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                          <  1:1571858/150‑1 (MQ=255)

                                                             tGGTATTAGCATTGATGCGTTAGATGATGGATATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                          <  2:225079/150‑1 (MQ=255)

                                                               gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                                                                        <  2:2299403/150‑1 (MQ=255)

                                                                tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                                                                       >  1:172476/1‑150 (MQ=255)

                                                                 aTTAGCATTGATGCGTTAGATGATGACTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                                                                      >  2:363579/1‑150 (MQ=255)

                                                                     gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                                                                  <  1:1765490/150‑1 (MQ=255)

                                                                      cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                                                                 <  2:493129/150‑1 (MQ=255)

                                                                            tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTg                                                                           <  2:841045/150‑1 (MQ=255)

                                                                             gCGTTAGATGATGGCCATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                                                                          <  1:997979/150‑1 (MQ=255)

                                                                                ttAGATGATGGCTATCTCACTCGATTCAGAGCATCCAACTCAGTGCTGGAAAGT‑AAAAAACCGAGGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTTATATGACGATTAGCCAGTGAAAGTATATAAATCGTCACTGctgt                                                                       <  1:62539/150‑4 (MQ=255)

                                                                                ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                       <  2:172476/150‑1 (MQ=255)

                                                                                ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                       <  1:481056/150‑1 (MQ=255)

                                                                                 tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAGCTCAGAACAGCAAAGT‑AAAAAAAAAAAGAAAAGTCGGGTTTATAACAGCCGGATTCGCAGAACGATTAATATGAAGATGACCCAGTGAAAGTATATAAATCGTCGCTGCggtt                                                                      >  2:1461595/1‑147 (MQ=255)

                                                                                   gatgatTGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                                                                    <  2:1261849/150‑1 (MQ=255)

                                                                                   gatgatGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                                                                    >  1:2446898/1‑150 (MQ=255)

                                                                                              tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGc                                                         >  1:2000676/1‑150 (MQ=255)

                                                                                               ctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTTCGGTATATATCGCAGTgac                                                        <  1:2378929/150‑3 (MQ=255)

                                                                                                 cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc                                                      <  2:24410/150‑1 (MQ=255)

                                                                                                 cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGAGATATATACCGAAGTGctcc                                                      >  2:785056/1‑150 (MQ=255)

                                                                                                     ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATACACCGAAGTGctccctcc                                                  >  1:168565/1‑150 (MQ=255)

                                                                                                      cAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCTAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg                                                            <  1:1412270/139‑1 (MQ=255)

                                                                                                        gTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGcc                                               >  2:1042670/1‑150 (MQ=255)

                                                                                                           agagCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGc                                            <  1:421847/150‑1 (MQ=255)

                                                                                                            gagCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCt                                           <  1:1012626/150‑1 (MQ=255)

                                                                                                             agCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTg                                          <  2:16390/150‑1 (MQ=255)

                                                                                                                caccaACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATCAAGATGACCCAGTCAATGAATATAAATCGTCACTGAGATATATACAGAAGTGCTCCCTCCCCAAGCTgaat                                       >  1:871690/1‑149 (MQ=25)

                                                                                                                 accaACTATTGGCTGGAAATT‑AAAAAAGCGACGCAATGTCGGTTTTTTTACATCCTGATTAGGATAAGGCTTAATATGACGATGACGAAGTTAAAGTATATAAATCGTCACTTCGATATATACGGAATTGCTCCCTCAGCCAGCTgaagt                                      <  2:2040899/150‑2 (MQ=255)

                                                                                                                       tCAGGGCTGGAAAGT‑AAAAAACCGACTCAAAGTCTGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGc                                <  1:1270463/150‑1 (MQ=255)

                                                                                                                         aGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTa                              <  1:673162/150‑1 (MQ=255)

                                                                                                                          gggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCTGATTCGGACAAGGATTAATATGAAGATGACCCAGTGAAAGTATATAAATATTCACTGAGATATATACAGAAGTGCTAAATAATCAAGATGAAGAAATCGCTaa                             >  2:1863956/1‑150 (MQ=21)

                                                                                                                          gggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTaa                             >  1:1223771/1‑150 (MQ=255)

                                                                                                                           ggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAAt                            >  2:1289239/1‑150 (MQ=255)

                                                                                                                               ggAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCtt                        >  2:1369443/1‑150 (MQ=255)

                                                                                                                               ggAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCtt                        >  1:208398/1‑150 (MQ=255)

                                                                                                                               ggAAAGT‑AAAAAACCCACGCAAAGTCGGTTTTTTTACATCCGGATTAGGACAACGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAAAGAATTGCTCACTCAGCCAGCTGAAGAAATCGCTAATTCtt                        >  2:2041709/1‑150 (MQ=38)

                                                                                                                                 gatgtaaaaaaaCCGACTTTGGGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTg                       >  1:2415806/7‑150 (MQ=255)

                                                                                                                                   aGT‑AAAAAACCGACGCAAAGTCGTTTTTTTTACATCCTGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCaa                    <  1:457097/150‑1 (MQ=255)

                                                                                                                                    gT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                   >  1:193249/1‑150 (MQ=255)

                                                                                                                                      aaaaaaaCCGACTTTGCGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATg                  >  1:112624/2‑150 (MQ=255)

                                                                                                                                        aaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGtt                <  1:2309317/150‑1 (MQ=255)

                                                                                                                                        aaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGtt                <  2:692975/150‑1 (MQ=255)

                                                                                                                                           aaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGGAATCGCTAATTCTTGCAATGTTAGc             <  2:1205472/150‑1 (MQ=255)

                                                                                                                                           aaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGc             >  1:298158/1‑150 (MQ=255)

                                                                                                                                             ccGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCa           >  1:209954/1‑150 (MQ=255)

                                                                                                                                                  gCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGc      <  2:324161/150‑1 (MQ=255)

                                                                                                                                                   cAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCt     <  1:1042670/150‑1 (MQ=255)

                                                                                                                                                      aGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTACCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAAt  >  2:1614946/1‑150 (MQ=255)

                                                                                                                                                      aGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACAGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAAt  <  1:224232/150‑1 (MQ=255)

TGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAAT  >  NC_000913/3647983‑3648281

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