BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,400,322	T→C	5.2%	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	5.2%	671.1 / 17.2	211	L12P (CTG→CCG)	hfq	RNA‑binding protein Hfq
Reads supporting (aligned to +/- strand): ref base T (103/97); new base C (9/2); total (112/99)
Fisher's exact test for biased strand distribution p-value = 6.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

Only 100 of 213 total aligned reads displayed.

TACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTG  >  NC_000913/4400174‑4400471

tACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTg                                                                                                                                                       <  2:2199422/150‑1 (MQ=255)

tACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTg                                                                                                                                                       >  2:2449230/1‑150 (MQ=255)

tACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTg                                                                                                                                                       <  1:2158490/150‑1 (MQ=255)

   aGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGaat                                                                                                                                                    <  1:81361/150‑2 (MQ=255)

    ggTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                                                                                                                                                   <  2:742071/150‑1 (MQ=255)

        gttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACt                                                                                                                                               <  1:832110/150‑1 (MQ=255)

        gttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACt                                                                                                                                               <  1:1805953/150‑1 (MQ=255)

             tGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                                                                                                                                                   <  2:1244338/141‑1 (MQ=255)

             tGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                                                                                                                                                   >  1:1244338/1‑141 (MQ=255)

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

              gCTATCGCGGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCg                                                                                                                                        >  2:853977/1‑150 (MQ=255)

              gCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCg                                                                                                                                        <  2:175599/150‑1 (MQ=255)

               cTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCgg                                                                                                                                       <  2:946470/150‑1 (MQ=255)

                   cGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAAc                                                                                                                                   >  1:448452/1‑150 (MQ=255)

                   cGCAGGCTGAATGTGAACAAATGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAAc                                                                                                                                   >  2:9215/1‑150 (MQ=255)

                       ggCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACgtgt                                                                                                                               <  2:44605/150‑1 (MQ=255)

                       ggCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCCG‑ACGCAGTGCGTTCAGGAACggat                                                                                                                               <  2:1426208/150‑4 (MQ=255)

                        gCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGtt                                                                                                                              <  1:416242/150‑1 (MQ=255)

                         cTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTc                                                                                                                             >  1:319108/1‑150 (MQ=255)

                             aTGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCATGAACGCAATGAG‑TCGGGAACGTGTTCCAGt                                                                                                                         >  1:2214771/1‑150 (MQ=255)

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

                                tgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTc                                                                                                                      <  1:2262306/150‑1 (MQ=255)

                                tgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTc                                                                                                                      <  1:2094438/150‑1 (MQ=255)

                                tgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTc                                                                                                                      <  2:2169925/150‑1 (MQ=255)

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

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

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

                                  tACTATTGAGACGTTTCTTGCGCAATTTTTTAAGAATGGATATGTTAAAAGTACAAATTAGAATATATGGTAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCAGTGCG‑TCGGTAACGGGTTCCAGTTTCta                                                                                                                    <  2:1120551/150‑1 (MQ=255)

                                    cAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtatt                                                                                                                  >  1:1193845/1‑150 (MQ=255)

                                    cAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTACAGTTTCtatt                                                                                                                  >  1:1287186/1‑150 (MQ=255)

                                     aaTTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtattt                                                                                                                 <  2:1268509/150‑1 (MQ=255)

                                      aTTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtattta                                                                                                                >  1:1356276/1‑150 (MQ=255)

                                       ttGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtatttat                                                                                                               >  1:1986258/1‑150 (MQ=255)

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

                                          agaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTg                                                                                                            <  2:1368239/150‑1 (MQ=255)

                                           gaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTgg                                                                                                           >  2:2022300/1‑150 (MQ=255)

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

                                                aTCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAACTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAt                                                                                                      <  1:673606/150‑1 (MQ=255)

                                                 tCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATg                                                                                                     <  2:463904/150‑1 (MQ=255)

                                                   gTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATATTTACAAGATCCGTTCCTGAAAGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGt                                                                                                   >  1:345299/1‑150 (MQ=255)

                                                       gcatattttaCAGAAACAAAATTTTAAAATTAAAAAAAAGCATATAAGCAAAAGAGAAAATAACTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTa                                                                                               >  2:1335245/6‑150 (MQ=255)

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

                                                        cAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTaa                                                                                              >  2:364883/1‑150 (MQ=255)

                                                        cAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTaa                                                                                              >  2:1212635/1‑150 (MQ=255)

                                                        cAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTaa                                                                                              >  2:2070902/1‑150 (MQ=255)

                                                          aTTTTTTCAGAATCGATAGGTTCAAAGTACAAATAATCATATAAGGAAAATAGATAATGGCTATGGGGCTATCTTTACAAGTTGCGTTCCTGAACGCACTGCT‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGc                                                                                            <  1:1033947/150‑1 (MQ=255)

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

                                                            tttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTg                                                                                          >  1:786798/1‑150 (MQ=255)

                                                            tttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCAATGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTg                                                                                          >  1:2177725/1‑150 (MQ=255)

                                                                 cAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAgg                                                                                     <  2:1042734/150‑1 (MQ=255)

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

                                                                     aTCGAAAGGTTAAAAGTACAAATAAGTATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATATTTACAAGATCCGTTCATGAACGCACTGAG‑TCGGGAACGTTTTCAAGATTCTAATTATTTGCTGAAAGGTATTTAGGTGCAAtggaga                                                                                 >  2:1079048/1‑147 (MQ=255)

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

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

                                                                           aGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAg                                                                           <  1:2022300/150‑1 (MQ=255)

                                                                               tCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCtt                                                                       >  2:1291804/1‑150 (MQ=255)

                                                                               tCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGCAACGAGTTCCAGTTACTATTTATTTGGTCAATGGTATTAAGCTCCAAGGGCAAATCGAGTCta                                                                       >  1:941848/1‑149 (MQ=255)

                                                                                  aaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTg                                                                    <  1:260397/150‑1 (MQ=255)

                                                                                  aaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑ACGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGAAAATCGAGTAttttt                                                                    >  1:2239375/1‑149 (MQ=255)

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

                                                                                        aaaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGt                                                              >  2:3096/1‑150 (MQ=255)

                                                                                           tAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCg                                                           <  2:412820/150‑1 (MQ=255)

                                                                                           tAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCg                                                           >  1:1111508/1‑150 (MQ=255)

                                                                                               cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTGCGTGAt                                                       <  2:1942948/150‑1 (MQ=255)

                                                                                                atatAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATc                                                      <  2:788871/150‑1 (MQ=255)

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

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

                                                                                                    aaGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGt                                                  >  1:1123471/1‑150 (MQ=255)

                                                                                                     aGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                     >  1:1686495/1‑146 (MQ=255)

                                                                                                      ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTg                                                >  1:1676163/1‑150 (MQ=255)

                                                                                                       gAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGa                                               <  2:712988/150‑1 (MQ=255)

                                                                                                            gagagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAac                                          <  2:971941/150‑1 (MQ=255)

                                                                                                            gagagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAac                                          <  1:1731109/150‑1 (MQ=255)

                                                                                                            gagagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAac                                          <  2:1397331/150‑1 (MQ=255)

                                                                                                                gaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACgg                                      <  2:916458/150‑1 (MQ=255)

                                                                                                                    ggCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTTCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAg                                  <  2:1897010/150‑1 (MQ=255)

                                                                                                                    ggCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAg                                  >  2:554012/1‑150 (MQ=255)

                                                                                                                        aaGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAg                              >  2:165441/1‑150 (MQ=255)

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

                                                                                                                             gCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGt                         >  1:1748151/1‑150 (MQ=255)

                                                                                                                              cAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGtt                        >  1:1640481/1‑150 (MQ=255)

                                                                                                                                aTCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTa                      <  1:2336213/150‑1 (MQ=255)

                                                                                                                                aTCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTa                      <  1:1136176/150‑1 (MQ=255)

                                                                                                                                 tCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                     >  1:1188320/1‑150 (MQ=255)

                                                                                                                                  cTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                    <  1:9215/150‑1 (MQ=255)

                                                                                                                                  cTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                    <  2:2138982/150‑1 (MQ=255)

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

                                                                                                                                       cAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAc               >  1:1508199/1‑150 (MQ=255)

                                                                                                                                       cAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAAAGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTTATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAAAAGCaa               >  2:842182/1‑149 (MQ=255)

                                                                                                                                        ggaacacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcg              >  2:1001921/7‑150 (MQ=255)

                                                                                                                                         gaacacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATACTGTTGAAAAACACGGTCAGCCAGAGGGTTTACAAGCCcgc             >  2:2113289/6‑150 (MQ=255)

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

                                                                                                                                           aTCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa           >  2:1668330/1‑150 (MQ=255)

                                                                                                                                            tCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAt          >  2:254031/1‑150 (MQ=255)

                                                                                                                                            cacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAt          >  1:991944/3‑150 (MQ=255)

                                                                                                                                               gTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTc       <  2:1779165/150‑1 (MQ=255)

                                                                                                                                                  ccTGAACGCTCTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc    <  1:195556/150‑1 (MQ=255)

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

                                                                                                                                                   cTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt   >  1:37067/1‑150 (MQ=255)

                                                                                                                                                    tGAAGGCAGTGCG‑TCGGGAACGTGTTCCAGTTTCTTTTTATTTGGTGAATGGTATTAAGCTGGAAGGGCAAATGGAGTCTTTTGATGTGTTGGTGATCCTGTTGAAAAAGACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg  <  2:941848/150‑1 (MQ=255)

TACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTG  >  NC_000913/4400174‑4400471

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