BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,412,794	A→T	6.8%	L136F (TTA→TTT)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,794	0	A	T	6.8%	1860.0 / 26.9	573	L136F (TTA→TTT)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base A (186/339); new base T (46/2); total (232/341)
Fisher's exact test for biased strand distribution p-value = 2.23e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 581 total aligned reads displayed.

GCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGG  >  NC_000913/4412646‑4412943

  gcgcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGc                                                                                                                                                    <  1:3366239/150‑1 (MQ=255)

   cgcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGcc                                                                                                                                                   <  1:9160532/150‑1 (MQ=255)

    gcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGGTTTTTGTTTAGCCg                                                                                                                                                  >  1:11071280/1‑150 (MQ=255)

     cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCgg                                                                                                                                                 <  1:5689599/150‑1 (MQ=255)

     cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCgg                                                                                                                                                 <  1:6601584/150‑1 (MQ=255)

        ggCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGtt                                                                                                                                              <  2:1274009/150‑1 (MQ=255)

         gCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTa                                                                                                                                             <  2:8930985/150‑1 (MQ=255)

         gCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGGTTTTTGTTTAGCCCGGTTa                                                                                                                                             >  1:5920366/1‑150 (MQ=255)

          cTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTAt                                                                                                                                            <  2:7133567/150‑1 (MQ=255)

          cTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGTGGTTTTTTTTTAAGCCGGGTTAt                                                                                                                                            >  2:286678/1‑150 (MQ=255)

           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATCCGGGTTAtt                                                                                                                                           >  2:9247300/1‑150 (MQ=255)

           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTAtt                                                                                                                                           >  2:4606938/1‑150 (MQ=255)

           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCCGGTTAtt                                                                                                                                           >  1:8375660/1‑150 (MQ=255)

              aacaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTTTTAATCCGGGGTATTAGt                                                                                                                                        >  2:3031915/1‑150 (MQ=255)

               acaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTa                                                                                                                                       <  2:5426696/150‑1 (MQ=255)

                caacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTaa                                                                                                                                      <  1:6836181/150‑1 (MQ=255)

                 aacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTaaa                                                                                                                                     <  1:7478826/150‑1 (MQ=255)

                 aacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTaaa                                                                                                                                     <  1:6290057/150‑1 (MQ=255)

                   cGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGATTTTTGTTAAGCCGGGTGATTCGGAAAGa                                                                                                                                   >  2:8734394/1‑150 (MQ=255)

                       aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATcgc                                                                                                                               <  2:8221956/150‑1 (MQ=255)

                       aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTATTACGGGGTATTAGTAAATATagc                                                                                                                               >  1:8036978/1‑150 (MQ=255)

                       aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTATTGATTTTTTTTTATCCGGGGTATTCGTAAAGATagc                                                                                                                               >  1:6174715/1‑150 (MQ=255)

                        aaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTAATCCGGGGTATTAGTAAATATagcg                                                                                                                              >  1:6438923/1‑150 (MQ=255)

                          aCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGAT‑CCGGTCGTTTCTGGATTTTTGTTAAGCCGGGGTATTCGTAAAGATCGCGATg                                                                                                                           >  2:8100171/1‑150 (MQ=255)

                            ggCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGCTTTTTTTTAAGCCCGGGTATTCGTAAAGATAGCGCTGa                                                                                                                          >  2:315000/1‑150 (MQ=255)

                            ggCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTAAGCCGGGTTATTAGTAAAGATCGCGAttt                                                                                                                          >  2:8005980/1‑148 (MQ=255)

                              cATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGc                                                                                                                        <  1:3786225/150‑1 (MQ=255)

                               aTTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGcc                                                                                                                       <  2:9156177/150‑1 (MQ=255)

                               aTTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGcc                                                                                                                       <  2:4947519/150‑1 (MQ=255)

                                ttttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCg                                                                                                                      <  2:8470063/150‑1 (MQ=255)

                                 tttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGAATTTTGTTAAGCCCGGGTATTCGTCAAGATAGCGATGAGccct                                                                                                                     >  2:5201717/1‑148 (MQ=255)

                                     tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtttgt                                                                                                                 <  1:2082989/150‑1 (MQ=255)

                                       ccGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTTGATTTTTTTTTAGACGGGTTATTAGTTAATATCGCGATGAGCCgttttttt                                                                                                               >  1:4095936/1‑150 (MQ=255)

                                         gATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTtat                                                                                                             <  1:8545830/150‑1 (MQ=255)

                                         gATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTtat                                                                                                             <  2:9429481/150‑1 (MQ=255)

                                           tGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATaa                                                                                                           <  2:452140/150‑1 (MQ=255)

                                                 gaagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAgcgg                                                                                                     <  2:4679695/150‑1 (MQ=255)

                                                 gaagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGATTTTTGGTTAGCCCGGGGAGTCGGAAAGAGCGCGCGGAGCCGGTGGGTTGTAACAgcgg                                                                                                     >  1:6919128/1‑150 (MQ=255)

                                                 gaagaaTATGACTCAGCAGAACCGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATCCGGGTTATTCGTAAAGATCGCGAGGGGCCGTTTTTTTTTTAAAgcgg                                                                                                     >  1:9916125/1‑150 (MQ=255)

                                                  aagatTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTTCGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAgcggc                                                                                                    <  2:6013354/150‑1 (MQ=255)

                                                  aagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTTTTTATCCGGGGTATTGGTAAAGATAGCGATGAGCCGTTTTTTTTTTACAgcggc                                                                                                    >  1:7882437/1‑150 (MQ=255)

                                                    gaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAAAAGcggcgg                                                                                                  >  2:6257466/1‑150 (MQ=255)

                                                       tatGACTCATCAGATCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAAGTGAGTTTGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCTTTTGTTTATAACAGCGGCGGTTc                                                                                               <  2:2461882/150‑1 (MQ=255)

                                                       tatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGTATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTc                                                                                               <  1:2920465/150‑1 (MQ=255)

                                                       tatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTc                                                                                               <  2:3382030/150‑1 (MQ=255)

                                                       tatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTc                                                                                               <  2:1439303/150‑1 (MQ=255)

                                                          gACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTTTTACGGGTTTTTGGTATAGATCGCGATGTGCCGTTTGTTTTTTAAAGcggggggtgcta                                                                                            >  1:10525753/1‑144 (MQ=255)

                                                          gACTCAGCAGAACGGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGCTTTTTGTTAAGCCGGGGGGTTCGGAAAGCTCGCGCTGGGCCGGTTGGTTATAACACCGGCGGGTCctc                                                                                            >  2:5493337/1‑150 (MQ=255)

                                                            cTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGTTTTTTTTTTAGGCGGGTTTTTAGTTAAGATAGCGCTGGGCCGGTTGTTTTTAAAAGCGGCGGTTACTCtt                                                                                          >  2:121965/1‑150 (MQ=255)

                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATCCGGGTTATTCGTAAAGATAGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCttt                                                                                         >  2:6430270/1‑150 (MQ=255)

                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCttt                                                                                         <  1:10967528/150‑1 (MQ=255)

                                                                 cagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCt                                                                                     <  1:3504775/150‑1 (MQ=255)

                                                                 cagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAAAAGCGGCGGTTCCTCTTTTGCt                                                                                     >  2:835504/1‑150 (MQ=255)

                                                                 cagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATAGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCATTTGCt                                                                                     >  2:5911461/1‑150 (MQ=255)

                                                                  agAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGGTTTTTGTTTAGCCGGGTTTTTAGTAAATATCGCCATGAGGCGGTTGTTTTTAAAAGCGGCGGTTGCTCgttttatt                                                                                    >  2:2023170/1‑146 (MQ=255)

                                                                   taaCTGTCAGTATAAGTACTAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTc                                                                                   <  1:9453786/149‑1 (MQ=255)

                                                                      cTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgcg                                                                                <  2:8545648/150‑1 (MQ=255)

                                                                       tGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTTAGCCGGGTTATTAGTAAATATCGCGCTGAGCCGTTTTTTTATAACAGCGGCGGTTCCTCATTTTCGTCgcgc                                                                               >  2:5628790/1‑150 (MQ=255)

                                                                            gTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGccagac                                                                          <  2:2217564/150‑5 (MQ=255)

                                                                             tCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTTTTGCGGGGTGTTGGGTAAAAAAGCGATTAGCCCGGTGTTTTTAACAGAGGCGGGTGCACGTTTTGTTAGCGcccccgc                                                                         >  1:7796476/1‑146 (MQ=255)

                                                                                aaGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTgg                                                                      <  1:1566069/150‑1 (MQ=255)

                                                                                aaGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTgg                                                                      <  1:7696894/150‑1 (MQ=255)

                                                                                     cGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATCCGGGTTATTCGTAAAGATCCCGATGAGCCGTTTGTTTATAAAAGCGGCGGTTCCACTTTTGCTTAGCGCCCACGCTGGCgcag                                                                 >  2:3643677/1‑150 (MQ=255)

                                                                                     cGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                                 <  1:4548260/150‑1 (MQ=255)

                                                                                     cGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                                 <  1:7220794/150‑1 (MQ=255)

                                                                                       aaaaTTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcagca                                                               <  1:4823544/150‑1 (MQ=255)

                                                                                            tGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCt                                                          <  2:7280908/150‑1 (MQ=255)

                                                                                              ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg                                                        <  2:4130368/150‑1 (MQ=255)

                                                                                               tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                                       <  1:5094733/150‑1 (MQ=255)

                                                                                                 ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTTATCCGGGTTATTCGTAAAGATAGCGGTGTGCCGTTTTTTTATAAAAGCGGCGGTTCCTCTTTTGCTTAGCGCCCAGTATGGCGCAGCACGTCTGGtta                                                     >  2:1865393/1‑148 (MQ=255)

                                                                                                 ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                                     <  1:2993280/150‑1 (MQ=255)

                                                                                                 ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                                     <  1:1988764/150‑1 (MQ=255)

                                                                                                 ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                                     <  1:5455922/150‑1 (MQ=255)

                                                                                                  tacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGAt                                                    <  1:2410726/150‑1 (MQ=255)

                                                                                                     gACAATGTTGAACAGAGTTGGATCCCGGTCGTTTATGGATTTTTTTTAAGCCGGGTTATTCGTAAATATAGCGATGATCCGTTTTTTTATAAAAGCGGCGGTTACTATTTTTCTTCGCGCCCACTCATGCGCAGCACATATGGTGATTaa                                                 >  2:6932046/1‑149 (MQ=255)

                                                                                                         aTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCt                                             <  2:9810021/150‑1 (MQ=255)

                                                                                                                 cAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTcc                                     <  2:1622408/150‑1 (MQ=255)

                                                                                                                  agagTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCg                                    <  2:748282/150‑1 (MQ=255)

                                                                                                                       tGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCt                               <  2:7370391/150‑1 (MQ=255)

                                                                                                                        ggATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                              <  1:8818877/150‑1 (MQ=255)

                                                                                                                        ggATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                              <  1:1319174/150‑1 (MQ=255)

                                                                                                                               ggTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaaa                       <  1:8177192/150‑1 (MQ=255)

                                                                                                                                gTCGTTTCTTGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAAAAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaaaa                      >  2:5823661/1‑150 (MQ=255)

                                                                                                                                gTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaaaa                      <  2:4847875/150‑1 (MQ=255)

                                                                                                                                gTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaaaa                      <  1:9335606/150‑1 (MQ=255)

                                                                                                                                  cGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga                    <  1:9660361/150‑1 (MQ=255)

                                                                                                                                    tttCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGt                  >  2:11075363/1‑150 (MQ=255)

                                                                                                                                       cTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCtt               <  1:2605892/150‑1 (MQ=255)

                                                                                                                                        tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCTCGCCCAGTATGGCGCAGCACATCTGGTGATTACAACTGGCGCTCCCGATCTGGCAAAAAAAATTCTTa              >  1:10679642/1‑150 (MQ=255)

                                                                                                                                        tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa              >  2:10783/1‑150 (MQ=255)

                                                                                                                                        tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa              >  2:7848064/1‑150 (MQ=255)

                                                                                                                                          gATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTAct            >  1:10042209/1‑150 (MQ=255)

                                                                                                                                           aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc           >  2:3943595/1‑150 (MQ=255)

                                                                                                                                            tttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGTTTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActct          >  1:9981081/1‑150 (MQ=255)

                                                                                                                                            tttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActct          >  2:10589706/1‑150 (MQ=255)

                                                                                                                                             ttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCtt         >  2:10310152/1‑150 (MQ=255)

                                                                                                                                              tttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTc        <  1:3294173/150‑1 (MQ=255)

                                                                                                                                              tttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTc        <  1:10245480/150‑1 (MQ=255)

                                                                                                                                               ttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCg       >  1:1693991/1‑150 (MQ=255)

                                                                                                                                                 gTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCggg     >  2:6656398/1‑150 (MQ=255)

GCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGG  >  NC_000913/4412646‑4412943

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