BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,400,322	T→C	6.9%	L12P (CTG→CCG)	hfq →	RNA‑binding protein Hfq

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,400,322	0	T	C	6.9%	349.8 / 12.1	115	L12P (CTG→CCG)	hfq	RNA‑binding protein Hfq
Reads supporting (aligned to +/- strand): ref base T (59/48); new base C (6/2); total (65/50)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 115 total aligned reads displayed.

CAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACT  >  NC_000913/4400176‑4400470

cAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGaa                                                                                                                                                    >  2:650348/1‑150 (MQ=255)

cAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGaa                                                                                                                                                    >  2:309084/1‑150 (MQ=255)

     tgttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCAc                                                                                                                                               <  1:137829/150‑1 (MQ=255)

     tgttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCAc                                                                                                                                               <  2:350597/150‑1 (MQ=255)

       ttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAATGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTg                                                                                                                                             <  2:1440169/150‑1 (MQ=255)

       ttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTg                                                                                                                                             <  1:199238/150‑1 (MQ=255)

        tgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGc                                                                                                                                            >  1:1423049/1‑150 (MQ=255)

           tGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑Tc                                                                                                                                        >  2:92367/1‑150 (MQ=255)

            gCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCg                                                                                                                                       >  1:301122/1‑150 (MQ=255)

               aTCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGa                                                                                                                                    >  2:380443/1‑150 (MQ=255)

                  gctggCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGATAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACg                                                                                                                                 >  2:1359102/4‑150 (MQ=255)

                  gCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACg                                                                                                                                 <  2:943942/150‑1 (MQ=255)

                  gCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCCG‑ACGCAGTGCGTTCAGGAACg                                                                                                                                 <  1:620907/150‑1 (MQ=255)

                        tGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACCTGTTcc                                                                                                                           >  2:806274/1‑150 (MQ=255)

                         gAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCa                                                                                                                          >  1:104645/1‑150 (MQ=255)

                             gtgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGttt                                                                                                                      <  2:481524/150‑1 (MQ=255)

                              tgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGtttg                                                                                                                     >  2:113502/1‑149 (MQ=255)

                               gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCggg                                                                                                                                     <  1:1418365/133‑1 (MQ=255)

                               gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCggg                                                                                                                                     >  2:1418365/1‑133 (MQ=255)

                               gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCt                                                                                                                    >  2:738802/1‑150 (MQ=255)

                               gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCt                                                                                                                    <  2:902738/150‑1 (MQ=255)

                               gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGCGTTCCAGTTTCt                                                                                                                    <  2:771765/150‑1 (MQ=255)

                                tACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCta                                                                                                                   >  1:1204170/1‑150 (MQ=255)

                                  cAATTGAGACGTATCGTGCGCAATTTTTTCAGATTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtatt                                                                                                                 <  1:320298/150‑1 (MQ=255)

                                    aTTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtattta                                                                                                               <  1:48028/150‑1 (MQ=255)

                                      tGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtatttatt                                                                                                             <  1:686575/150‑1 (MQ=255)

                                       gagaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTatttattt                                                                                                            >  2:992667/1‑150 (MQ=255)

                                         gaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTgg                                                                                                          <  1:992667/150‑1 (MQ=255)

                                          aCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGt                                                                                                         <  2:1187161/150‑1 (MQ=255)

                                           cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTg                                                                                                        <  2:1512561/150‑1 (MQ=255)

                                            gTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGa                                                                                                       <  1:286302/150‑1 (MQ=255)

                                              aTCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAt                                                                                                     >  2:1288804/1‑150 (MQ=255)

                                              aTCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAt                                                                                                     >  1:675327/1‑150 (MQ=255)

                                                  tGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTa                                                                                                 <  1:315106/150‑1 (MQ=255)

                                                  tGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTa                                                                                                 <  1:380443/150‑1 (MQ=255)

                                                   gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAt                                                                                                >  1:808757/1‑150 (MQ=255)

                                                   gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAt                                                                                                <  2:359593/150‑1 (MQ=255)

                                                     gcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTa                                                                                              >  2:776469/1‑150 (MQ=255)

                                                        aTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGc                                                                                           >  2:648667/1‑150 (MQ=255)

                                                         ttttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCt                                                                                          <  2:91470/150‑1 (MQ=255)

                                                           ttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGc                                                                                        >  1:886000/1‑150 (MQ=255)

                                                            tttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCa                                                                                       >  1:729177/1‑150 (MQ=255)

                                                             ttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCaa                                                                                      >  1:1068211/1‑150 (MQ=255)

                                                              tCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAg                                                                                     <  2:289557/150‑1 (MQ=255)

                                                               cAGAATCCAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAAAGTGTTACAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAgg                                                                                    >  2:603108/1‑150 (MQ=255)

                                                                aGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAggg                                                                                   <  1:738802/150‑1 (MQ=255)

                                                                aGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAggg                                                                                   >  2:914633/1‑150 (MQ=255)

                                                                   aTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCaa                                                                                >  1:594444/1‑150 (MQ=255)

                                                                     cGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAAt                                                                              <  2:675327/150‑1 (MQ=255)

                                                                      gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATc                                                                             >  1:586570/1‑150 (MQ=255)

                                                                           gTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTc                                                                        <  2:829581/150‑1 (MQ=255)

                                                                              cAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTACTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCttt                                                                     >  1:809496/1‑150 (MQ=255)

                                                                               aaaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCtttt                                                                    >  1:316346/1‑150 (MQ=255)

                                                                                  gTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                 >  1:57319/1‑150 (MQ=255)

                                                                                   tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                >  1:1098300/1‑150 (MQ=255)

                                                                                   tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                >  1:457141/1‑150 (MQ=255)

                                                                                   tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                <  1:913416/150‑1 (MQ=255)

                                                                                   tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                >  1:622430/1‑150 (MQ=255)

                                                                                       aaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGtt                                                            >  2:1080882/1‑150 (MQ=255)

                                                                                       aaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGtt                                                            >  2:931401/1‑150 (MQ=255)

                                                                                        aTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTc                                                           >  1:1041739/1‑150 (MQ=255)

                                                                                          aaGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGt                                                         <  2:1204170/150‑1 (MQ=255)

                                                                                          aaGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGt                                                         >  2:695169/1‑150 (MQ=255)

                                                                                           aGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTg                                                        <  2:104645/150‑1 (MQ=255)

                                                                                             cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGAt                                                      <  1:776469/150‑1 (MQ=255)

                                                                                               tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                    <  2:886000/150‑1 (MQ=255)

                                                                                               tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                    >  1:250099/1‑150 (MQ=255)

                                                                                                ataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCt                                                   >  1:1121513/1‑150 (MQ=255)

                                                                                                 taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTg                                                  >  1:889960/1‑150 (MQ=255)

                                                                                                   aGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtt                                                <  1:294338/150‑1 (MQ=255)

                                                                                                      aaaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGaa                                             >  2:1461701/1‑150 (MQ=255)

                                                                                                         agagagAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTg                                               >  1:538318/1‑145 (MQ=255)

                                                                                                          gagagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAac                                         >  1:366383/1‑150 (MQ=255)

                                                                                                            gagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAacac                                       <  1:650348/150‑1 (MQ=255)

                                                                                                            gagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAacac                                       >  1:419610/1‑150 (MQ=255)

                                                                                                                aTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACAGTc                                   <  1:806274/150‑1 (MQ=255)

                                                                                                                   gCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGc                                <  2:1090320/150‑1 (MQ=255)

                                                                                                                       aGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  2:27503/150‑1 (MQ=255)

                                                                                                                        ggggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGAt                           >  1:998482/1‑150 (MQ=255)

                                                                                                                         gggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                          >  1:1312988/1‑150 (MQ=255)

                                                                                                                         aacggaTCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                          >  1:879795/6‑150 (MQ=255)

                                                                                                                          ggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATgg                         >  1:849339/1‑150 (MQ=255)

                                                                                                                          ggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATgg                         >  1:426514/1‑150 (MQ=255)

                                                                                                                             aaTCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGttt                      >  1:644735/1‑150 (MQ=255)

                                                                                                                                cTTTACAAGACCCGTTCCTGAACGCACTGCG‑TCTGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                   >  2:256172/1‑150 (MQ=255)

                                                                                                                                  ttACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAg                 <  1:1288804/150‑1 (MQ=255)

                                                                                                                                    aCAAGATCCGTTGCTGAACGCACTGGG‑TCGGGAAGGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCATGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCa               <  2:1121513/150‑1 (MQ=255)

                                                                                                                                    aCAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCa               <  1:1461701/150‑1 (MQ=255)

                                                                                                                                    aCAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCa               >  2:1518132/1‑150 (MQ=255)

                                                                                                                                      ggaacacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcg             >  1:1242568/7‑150 (MQ=255)

                                                                                                                                       aGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgc            >  2:1124517/1‑150 (MQ=255)

                                                                                                                                        aacacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgcg           >  2:392598/5‑150 (MQ=255)

                                                                                                                                         acacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa          >  1:800277/4‑150 (MQ=255)

                                                                                                                                         acacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa          >  2:1216072/4‑150 (MQ=255)

                                                                                                                                           acGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAtt        >  2:58710/2‑150 (MQ=255)

                                                                                                                                            cGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAttt       >  1:241350/1‑150 (MQ=255)

                                                                                                                                               tCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa    >  2:850014/1‑150 (MQ=255)

                                                                                                                                                ccTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc   >  2:1344937/1‑150 (MQ=255)

                                                                                                                                                ccTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc   <  1:44645/150‑1 (MQ=255)

                                                                                                                                                 cTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt  <  1:1344937/150‑1 (MQ=255)

CAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACT  >  NC_000913/4400176‑4400470

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