BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,648,130	A→T	5.5%	noncoding (68/84 nt)	agrA →	inactive antisense sRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,648,130	0	A	T	5.5%	377.0 / 9.1	126	noncoding (68/84 nt)	agrA	inactive antisense sRNA
Reads supporting (aligned to +/- strand): ref base A (40/79); new base T (4/3); total (44/82)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

Only 100 of 127 total aligned reads displayed.

CTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCT  >  NC_000913/3647982‑3648278

  ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAg                                                                                                                                                   >  2:1092837/1‑150 (MQ=255)

  ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTAAGGGCTGGAAAGTAAAAAACCGACGCAAAg                                                                                                                                                   >  2:1583663/1‑150 (MQ=255)

     tAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGATGCAAAGTCg                                                                                                                                                >  2:942795/1‑150 (MQ=255)

       aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCATTAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                                                              >  1:1809769/1‑150 (MQ=255)

       aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCTGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACGAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                                                              >  1:106977/1‑150 (MQ=255)

       aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                                                              <  2:1224316/150‑1 (MQ=255)

       aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                                                              >  1:1562730/1‑150 (MQ=255)

           tACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                                                                          <  2:1305616/150‑1 (MQ=255)

           tACAGAACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                                                                          >  1:1974792/1‑150 (MQ=255)

            aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCATCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                         <  1:1998673/150‑1 (MQ=255)

            aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGtttttt                                                                                                                                         <  2:848855/150‑1 (MQ=38)

            aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                         <  2:220883/150‑1 (MQ=255)

             cAGTACATTCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                        <  2:1298543/150‑1 (MQ=255)

             cAGTACATGCAGGTTATAAAACCCTCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                        <  2:1247460/150‑1 (MQ=255)

             cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                        <  2:1851871/150‑1 (MQ=255)

             cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                        <  1:337272/150‑1 (MQ=255)

              aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                       <  2:1381902/150‑3 (MQ=38)

              aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                                                                       <  1:472188/150‑1 (MQ=255)

              aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                                                                       >  1:2000296/1‑150 (MQ=255)

              aGTACATGCAGGTTATAAAACCAGCACGGCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                                                                       <  1:694181/150‑1 (MQ=255)

              aGTACATGCAGGGTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                       <  1:1088202/150‑3 (MQ=38)

               gTACATGCAGGTTAAAACACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTATATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTAc                                                                                                                                      <  2:953303/150‑1 (MQ=255)

                    tGCAGGTTATAAAACCAGCACGTCCTTGCATTAGTTTCAGTATGGTTATAGCATTGATGCGTTTGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCg                                                                                                                                 <  2:1530269/150‑1 (MQ=255)

                    tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCg                                                                                                                                 >  1:149460/1‑150 (MQ=255)

                          ttATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCg                                                                                                                           <  2:1219073/150‑1 (MQ=255)

                           tataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCACCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGGTTTTTTACATCAGGATTCgg                                                                                                                          >  1:1110210/1‑150 (MQ=255)

                           tataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCgg                                                                                                                          >  2:574784/1‑150 (MQ=255)

                           tataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCgg                                                                                                                          <  2:404136/150‑1 (MQ=255)

                            ataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGa                                                                                                                         <  2:1889127/150‑1 (MQ=255)

                            ataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCAACCAACTGCGGGCTGGAAAGTAAAAAACAGCCGCAAAGTCGGTTTGTTTACAGACGGATACggc                                                                                                                         >  1:1388031/1‑149 (MQ=255)

                             tataaacAGCAATTCCTTGCAATTTTTTCATTTTGGTATTAGCATTGATGTGTTAGATTATTTCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAc                                                                                                                        <  2:1080420/144‑1 (MQ=14)

                              aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGACGGTTTTTTTACATCCGGATTCGGACa                                                                                                                       >  2:323676/1‑150 (MQ=255)

                              aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGCATAGTAAAAAACCGAGGAAAAGACGGTATATTTACATCCGGAATCgcaaa                                                                                                                       >  1:683802/1‑146 (MQ=37)

                                aaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAg                                                                                                                     <  2:1815606/150‑1 (MQ=255)

                                 aCCAGCACGTCCTTGAAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCATGGATGGAAAGTAAAAAACAGACGAAAAGTCGGTTTTTTTACATCAGGATTCTGACATgg                                                                                                                    >  1:1894796/1‑150 (MQ=255)

                                    aGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCtt                                                                                                                 >  2:544960/1‑150 (MQ=255)

                                     gCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGTGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGTCTTa                                                                                                                >  1:1210004/1‑150 (MQ=255)

                                     gCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTa                                                                                                                >  2:951980/1‑150 (MQ=255)

                                       aCGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGTGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAat                                                                                                              <  1:627567/150‑1 (MQ=255)

                                          tCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCTACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATg                                                                                                           >  1:686235/1‑150 (MQ=255)

                                           ccTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGTAAATTAAAAAACCGACTGAAAGACAGTATAATTACAAACGCATACGCACAAGGCTTAatataa                                                                                                          >  1:1759215/1‑148 (MQ=25)

                                             ttGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                                                                        <  1:264646/150‑1 (MQ=255)

                                              tGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                                                                       <  2:1485604/150‑1 (MQ=255)

                                               gCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGTAAAAGGCTTAATATGACGAt                                                                                                      >  1:362838/1‑150 (MQ=255)

                                                       ttCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGt                                                                                              <  1:1340743/150‑1 (MQ=255)

                                                           ttATGGTATTAGCATTGAGGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaa                                                                                          <  2:1660851/149‑1 (MQ=255)

                                                           gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCTGATTCTGACATGGATTAATATGACGATGAACCTGTGaaa                                                                                          >  1:1946977/1‑150 (MQ=25)

                                                            tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                                                                         <  2:1931491/150‑1 (MQ=255)

                                                             aTGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGt                                                                                        <  1:355393/150‑1 (MQ=255)

                                                              tGGTATTAGGATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCATAGCCACCTACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATGCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                       <  2:1346731/150‑1 (MQ=255)

                                                              tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                       <  1:1281911/150‑1 (MQ=255)

                                                              tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                       <  1:1885495/150‑1 (MQ=255)

                                                               ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                                                                      >  1:660585/1‑150 (MQ=255)

                                                                 tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                                                                    <  2:1399563/150‑1 (MQ=255)

                                                                      gCATTGATGCGGTAGATGATGTCTATCTCACTCCAGTAAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGAGTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                                                               <  1:1363687/150‑1 (MQ=255)

                                                                         ttGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGACACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTc                                                                            <  1:1476336/150‑1 (MQ=255)

                                                                           gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCAc                                                                          >  2:1100525/1‑150 (MQ=255)

                                                                               cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                                                                      <  2:262722/150‑1 (MQ=255)

                                                                                gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGa                                                                     <  1:1490960/150‑1 (MQ=255)

                                                                                 ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGATGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                    <  1:942795/150‑1 (MQ=255)

                                                                                 ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                    <  2:1961286/150‑1 (MQ=255)

                                                                                  tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                                                                   <  1:1526169/150‑1 (MQ=255)

                                                                                   aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatat                                                                  >  2:1377443/1‑150 (MQ=255)

                                                                                      tgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                                                               >  1:861415/1‑150 (MQ=255)

                                                                                       gatgGCTAGCTCACTCCATGCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGTATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc                                                              <  1:1143238/150‑1 (MQ=255)

                                                                                        ttgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                                                             <  2:765845/149‑1 (MQ=255)

                                                                                        atgGCTATCTCTCTCCAGTCAGATCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                                                             <  1:439829/150‑1 (MQ=255)

                                                                                        atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                                                             <  2:501795/150‑1 (MQ=255)

                                                                                         tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCg                                                            <  2:564620/150‑1 (MQ=255)

                                                                                          ggCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGa                                                           <  1:1086944/150‑1 (MQ=255)

                                                                                            cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg                                                         <  1:408832/150‑1 (MQ=255)

                                                                                              aTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAATTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTg                                                       <  1:1919155/150‑1 (MQ=255)

                                                                                                 tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc                                                    <  2:445984/150‑1 (MQ=255)

                                                                                                    cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                                                                           >  1:128815/1‑124 (MQ=255)

                                                                                                    cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                                                                           <  2:128815/124‑1 (MQ=255)

                                                                                                      acAGTCAGATCCACCAACTCAGGGCTGGAAAGTAAAATACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc                                               <  1:168361/149‑1 (MQ=255)

                                                                                                       cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCg                                              <  2:544866/150‑1 (MQ=255)

                                                                                                        aGTCAGAGCCACCAACTAATGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATTCCGAAGTGCTCCCTCCGc                                             <  2:1500323/150‑1 (MQ=255)

                                                                                                          taagagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCa                                           <  1:251028/148‑1 (MQ=255)

                                                                                                              agCCACCAACTCAGTGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTg                                       <  2:1413689/150‑1 (MQ=255)

                                                                                                                 caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                    <  1:1880038/150‑1 (MQ=255)

                                                                                                                 caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                    <  2:1809769/150‑1 (MQ=255)

                                                                                                                 caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                    <  2:1142784/150‑1 (MQ=255)

                                                                                                                      aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATc                               <  1:289993/150‑1 (MQ=255)

                                                                                                                       cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCg                              <  2:350800/150‑1 (MQ=255)

                                                                                                                       atCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCg                              <  1:1643942/149‑1 (MQ=255)

                                                                                                                        tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCTGATTCTGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGc                             <  1:323676/150‑1 (MQ=255)

                                                                                                                         cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCt                            <  1:1747850/150‑1 (MQ=255)

                                                                                                                                  aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                   <  1:1385157/150‑1 (MQ=255)

                                                                                                                                  aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                   <  1:574784/150‑1 (MQ=255)

                                                                                                                                     taaaaaaaCCGACTTTGCGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                >  2:1507835/3‑150 (MQ=255)

                                                                                                                                     gTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCTATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                <  1:1477650/150‑1 (MQ=255)

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

                                                                                                                                       aaaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAATGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGt              >  2:1612441/1‑150 (MQ=255)

                                                                                                                                        aaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGtt             >  1:902503/1‑150 (MQ=255)

                                                                                                                                         aaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTa            >  1:452054/1‑150 (MQ=255)

                                                                                                                                                 cGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTgg    <  1:1152507/150‑1 (MQ=255)

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

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

                                                                                                                                                   cAAAGTCGGTTTTTTTACATCCGGATTCGGTCAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCTTCACTGCGATATATACCTAAGTGCTCCCTCCGCCAGCTGTAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCt  <  1:1025695/150‑1 (MQ=255)

CTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCT  >  NC_000913/3647982‑3648278

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