BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,054,888	A→T	7.5%	V186E (GTA→GAA)	ygfB ←	UPF0149 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,054,888	0	A	T	7.5%	1677.7 / 15.3	532	V186E (GTA→GAA)	ygfB	UPF0149 family protein
Reads supporting (aligned to +/- strand): ref base A (157/324); new base T (46/3); total (205/327)
Fisher's exact test for biased strand distribution p-value = 1.17e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 545 total aligned reads displayed.

TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037

 aCTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTGTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAc                                                                                                                                                      >  2:4255117/1‑150 (MQ=255)

   ttCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACtt                                                                                                                                                    <  1:3046883/150‑1 (MQ=255)

      tGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTg                                                                                                                                                 <  1:5529361/150‑1 (MQ=255)

       ggTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTATgg                                                                                                                                                >  1:258803/1‑150 (MQ=255)

        gTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGgc                                                                                                                                               <  2:2910117/150‑1 (MQ=255)

          gcagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTATGgagc                                                                                                                                             <  2:7591896/150‑1 (MQ=255)

           cagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGcgcg                                                                                                                                            <  2:3190077/150‑1 (MQ=255)

             gcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCTGGCGcggt                                                                                                                                          >  2:173559/1‑150 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCATGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtc                                                                                                                                         <  1:4525641/150‑1 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtc                                                                                                                                         <  2:5157659/150‑1 (MQ=255)

               aaaaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCATTACGTTTTTTGTTTTTAATGTAGAGTCGGTTTTTTTTATTCTTGAGcggtag                                                                                                                                        >  2:9257531/1‑148 (MQ=255)

                aaaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtcgg                                                                                                                                       <  2:2755914/150‑1 (MQ=255)

                  aaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGTTATCTCACTCATAAAACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTCTGGCGCGGTCGGtt                                                                                                                                     <  2:5442279/150‑1 (MQ=255)

                   aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTCTTGTGGGGCGGTCGGGTg                                                                                                                                    >  2:5582470/1‑150 (MQ=255)

                     cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                  <  2:4498494/150‑1 (MQ=255)

                     cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                  <  2:6993415/150‑1 (MQ=255)

                     cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                  <  2:6473802/150‑1 (MQ=255)

                      aGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTCCGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTACTTGTGGCGCgggcgggggcgg                                                                                                                                 >  1:1362916/1‑150 (MQ=255)

                      aGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTATGGCGCGGTCGGTTGCgg                                                                                                                                 >  1:318309/1‑150 (MQ=255)

                         gcggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATg                                                                                                                              <  2:8844564/150‑1 (MQ=255)

                          cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  2:4753514/150‑1 (MQ=255)

                          cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  1:1916693/150‑1 (MQ=255)

                          cggcgCTGCCGGGGTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTGGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  2:9046289/150‑1 (MQ=255)

                             cgcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTATGGCGCGGTCGGTTGCGGGTGAGTa                                                                                                                          >  1:3417525/1‑150 (MQ=255)

                              gcTGCCGGGTTGCCTGTGCTCCACCAGGGCCTGCCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATCACCCTCTCCTCACGTTTTTTGATTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTaa                                                                                                                         <  1:8487790/150‑1 (MQ=255)

                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  2:2338719/150‑1 (MQ=255)

                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  2:3658469/150‑1 (MQ=255)

                                  ccGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGgt                                                                                                                     <  1:6100848/150‑1 (MQ=255)

                                  ccGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGgt                                                                                                                     <  1:2380314/150‑1 (MQ=255)

                                   cGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGgtg                                                                                                                    <  1:8263245/150‑1 (MQ=255)

                                      gTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCg                                                                                                                 <  1:9249381/150‑1 (MQ=255)

                                       ttGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCAGGCGCGGTCGGTTGCGGATGAGTAAACGTGTCGt                                                                                                                >  2:6825246/1‑150 (MQ=255)

                                       ttGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGt                                                                                                                <  1:5180582/150‑1 (MQ=255)

                                         gCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTgg                                                                                                              <  1:3644071/150‑1 (MQ=255)

                                          cATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGc                                                                                                             <  2:1080431/150‑1 (MQ=255)

                                            tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCAt                                                                                                           <  2:4637146/150‑1 (MQ=255)

                                            tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCAt                                                                                                           >  2:8434329/1‑150 (MQ=255)

                                              tGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATaa                                                                                                         <  2:2220173/150‑1 (MQ=255)

                                                cTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACa                                                                                                       <  2:4017111/150‑1 (MQ=255)

                                                   caccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAAGACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAgcg                                                                                                    <  2:82385/150‑1 (MQ=255)

                                                    accaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAgcgc                                                                                                   <  1:7241256/150‑1 (MQ=255)

                                                     ccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGcgcg                                                                                                  <  1:5012043/150‑1 (MQ=255)

                                                      caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCgg                                                                                                 <  1:4307851/150‑1 (MQ=255)

                                                      caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGGTGCCGATGTGTAAAGGTGTCGTGGCATAACAGCGCgg                                                                                                 >  2:6264220/1‑150 (MQ=255)

                                                         ggCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATAATTGCGCGGTCGGTTGGGGCTGATTAAAGGTGTCGTGGCAAAAAACCGCGGCaa                                                                                              >  2:7053010/1‑150 (MQ=255)

                                                             tGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTCTTGTTGCGGGGTAGGGTGCGGATGTGTAAAGGTGTAGTGTCATAAAAGCGCGGCaaaaac                                                                                          >  1:5645899/1‑146 (MQ=255)

                                                             tGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTGATTAAGGGGCGGTGGGTTTCGGTTGAGGAAAGGTGTGGTGGGATAAAAGCCCGGCAacacc                                                                                          >  1:4101412/1‑149 (MQ=255)

                                                               aCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACATCTGGCGCGGTCGGTTGCGGATGATTAAAAGTGTCGTGGCCTAACAGCGCGGCAACACGaa                                                                                        >  2:8377207/1‑150 (MQ=255)

                                                                cGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAAc                                                                                       <  2:9334203/150‑1 (MQ=255)

                                                                    gACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTACATATTGCGCGGTCGGGTGCGGATGTGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTAt                                                                                   >  1:6178113/1‑150 (MQ=255)

                                                                    gACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTG‑TTTTAGTGTAGAGTCGGGTTTTTTTGTTGGTGGGCGGGCGGGTGCGGGTGGGTTAAGGGGTCGTGGCATAAAAGCGCGGCCACACCAAcctgtt                                                                                  >  1:4754104/1‑145 (MQ=255)

                                                                      cGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATTACAGCGCGGCAACACGAACGTATTc                                                                                 <  1:2911863/150‑1 (MQ=255)

                                                                       gCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAGTTATTGCGGGGGCGGGTGGGGGGGGGTAAATGTGTGGTGGGATAACAGAGCGGCACCACCAAcagattag                                                                                >  1:3288554/1‑143 (MQ=255)

                                                                        cTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTCTGGGGCGGTCGGGTGCGGGTGAGTAAAGGTGTGGTGGCCTAAAAGCGCGGCAACACCAAAGGATTCga                                                                               >  1:9140933/1‑150 (MQ=255)

                                                                         tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgat                                                                              <  1:6498340/150‑1 (MQ=255)

                                                                              ctctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGAtct                                                                         <  1:6319079/150‑1 (MQ=255)

                                                                               tctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGAtctc                                                                        <  2:3319948/150‑1 (MQ=255)

                                                                                ctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATATCATCGCGGATACACGAACGTATTCGTTGATctct                                                                       <  2:8146002/150‑1 (MQ=255)

                                                                                ctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATctct                                                                       <  2:7932701/150‑1 (MQ=255)

                                                                                ctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATctct                                                                       <  1:8703782/150‑1 (MQ=255)

                                                                                ctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGGAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATctct                                                                       <  1:6928103/150‑1 (MQ=255)

                                                                                    ccGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTTTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCaa                                                                   <  1:7425498/150‑1 (MQ=255)

                                                                                      gggATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGc                                                                 <  2:6704837/150‑1 (MQ=255)

                                                                                         atatCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGAc                                                              <  2:6980141/150‑1 (MQ=255)

                                                                                            tctcACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTCTTCTGGCGCGGTCGGGTGCGGGTGAGTAAAGGTGTCGTGGCATAAAACCGCGGCAACACGAACGTATTCGATGATATCTTGAAGCGACagt                                                           >  2:2039509/1‑148 (MQ=255)

                                                                                            tctcACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACAtt                                                           <  1:1009692/150‑1 (MQ=255)

                                                                                            tctcACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTCTTATGGGGCGGGGGGGTGCGGGGGTGTTAAGGTGTCGTGGCCTAAAAGCGCGGCCACACCAACGGATTCGATAAACACATAAAGCGACaat                                                           >  2:4214082/1‑148 (MQ=255)

                                                                                             ctcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACAttt                                                          <  2:5916873/150‑1 (MQ=255)

                                                                                             ctcactcaTAACACTCTCCGTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACAttt                                                          <  2:5242307/150‑1 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTc                                                         <  2:7658003/150‑1 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTc                                                         <  1:5078580/150‑1 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGGTGCGGCTGTGTAAAGGGGTCGTGGCCTAACAGCGCGGCAACACGAACGTATTCGATGATCGCTTCCAGCGCCATTTc                                                         >  2:4970042/1‑150 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGCGGTCGGTTGCGGGTGTGGAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACCTATTCGATGATCTCTTCAAGCGACATTTc                                                         >  2:538303/1‑150 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACATATGGCGCGGTCGGTTGCGGCTGAGTAAAGGTGTCGTGGCATAAAAGCGCGGCAAAACGAACGTATTAGATGATCTATTAAAGCGACAatta                                                         >  2:3087747/1‑149 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTTCGGATTAGTAAATGTGTAGTGGCATAAAAGCGCGGCAACACCAACGTATTAGATTATATCTTAAAACGACATTTCaa                                                       >  2:7131144/1‑150 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGCGGTCGGGTGCGGCTGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       >  2:5757804/1‑150 (MQ=255)

                                                                                                 ctcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAg                                                      <  1:7618159/150‑1 (MQ=255)

                                                                                                          ctctCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTg                                             <  2:192926/150‑1 (MQ=255)

                                                                                                             tcCTTACGTTTTTTGTTTTTAGTTTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTAAAGCGACATTTCAAGCTCTTCCTGAtc                                          <  1:7121378/150‑1 (MQ=255)

                                                                                                                ttACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttc                                       <  1:2137501/150‑1 (MQ=255)

                                                                                                                             ttttAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTaa                          <  2:3515795/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  2:1258483/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  2:3784661/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  2:734379/150‑1 (MQ=255)

                                                                                                                                        agtcttattttGTACTTCTGGCTCGGTCTGTTTCGGATTAGTAAAGGTGTCGTGGCATAACAGCTCGGCAACACGAACGTATTAGATGATCACTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTgcgc               <  1:7000933/143‑1 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATTAGTAAAGGTGTAGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAATCTCTTACAGATCTTCGTCTTCGTCGTAACCCAGTTTCGCaa             >  1:149218/1‑150 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatg                                                                             <  1:8894454/86‑1 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatg                                                                             >  2:8894454/1‑86 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAAAAGCGCGGCAACACGAACGTATTCGATGATCTCATAAAGCGACAATTAAAGCTCTTACATATCTTAGTCTTCGTCGTAACCCAGTTTCGCaa             >  1:3827725/1‑150 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCAAAACAGCGCCGCAACACGAACCTATTCGATGATCACTTCAAGCGACAATTCAAGCTCTTCCAGATATTCGTCTTCGTCGTAACCCACTTTCGCaa             >  1:7720750/1‑150 (MQ=255)

                                                                                                                                           cGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAAt            <  2:157391/150‑1 (MQ=255)

                                                                                                                                           cGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAAt            <  1:1716938/150‑1 (MQ=255)

                                                                                                                                             gTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGt          >  2:5841827/1‑150 (MQ=255)

                                                                                                                                             gTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGt          >  1:8100086/1‑150 (MQ=255)

                                                                                                                                              tttttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGtt         <  1:2241256/150‑1 (MQ=255)

                                                                                                                                                atagTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGAGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTACAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTAc       <  2:6119363/147‑1 (MQ=255)

                                                                                                                                                 ttGTACTTCTGGCGCGGTCGGTTGCGGATGAGGAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACg      <  2:2577857/150‑1 (MQ=255)

                                                                                                                                                    tACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAg   >  2:9116538/1‑150 (MQ=255)

                                                                                                                                                     aCTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGa  <  1:3957947/150‑1 (MQ=255)

                                                                                                                                                     aCTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGa  <  1:539114/150‑1 (MQ=255)

TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037

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