BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,054,888	A→T	7.9%	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.9%	1459.7 / 27.1	464	V186E (GTA→GAA)	ygfB	UPF0149 family protein
Reads supporting (aligned to +/- strand): ref base A (146/275); new base T (38/2); total (187/277)
Fisher's exact test for biased strand distribution p-value = 1.88e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 468 total aligned reads displayed.

TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037

tACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCTGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTa                                                                                                                                                       <  2:9074428/150‑1 (MQ=255)

     cTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCt                                                                                                                                                  >  2:1517758/1‑150 (MQ=255)

     cTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCt                                                                                                                                                  <  1:9025874/150‑1 (MQ=255)

     cTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCt                                                                                                                                                  <  2:7831117/150‑1 (MQ=255)

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

         tGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTATTgcg                                                                                                                                              >  1:8251749/1‑150 (MQ=255)

            agcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTCTTGCGcgg                                                                                                                                           >  1:5983493/1‑150 (MQ=255)

             gcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggt                                                                                                                                          <  2:7094941/150‑1 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTCCGTTTTTTGTCTTTAGTGTAGAGTCGGGGTTTGTTGTTCTGgcggcggg                                                                                                                                         >  2:8492378/1‑145 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtc                                                                                                                                         >  1:1367829/1‑150 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGcggtc                                                                                                                                         >  1:5366069/1‑150 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTAGTTATGGCGcggtc                                                                                                                                         >  1:3160133/1‑150 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGGACTTCTGGCGcggtc                                                                                                                                         >  2:1023263/1‑150 (MQ=255)

               aaaaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTCTTGTGGGGcggggg                                                                                                                                        >  1:2479312/1‑147 (MQ=255)

               aaaaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTCTTCTGGCGcgggcg                                                                                                                                        >  1:4618772/1‑150 (MQ=255)

                 aaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTATTGTGGGGcgggcggt                                                                                                                                      >  2:1255933/1‑150 (MQ=255)

                   aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTCTGGCGCGGTCGGGTg                                                                                                                                    >  1:6427233/1‑150 (MQ=255)

                    tCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGGTGc                                                                                                                                   >  2:5603892/1‑150 (MQ=255)

                    tCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTGTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGGCGGTTGc                                                                                                                                   >  2:5134724/1‑150 (MQ=255)

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

                      aGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCgg                                                                                                                                 <  2:8533129/150‑1 (MQ=255)

                        cgcgGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTTTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGAt                                                                                                                               <  2:6691260/150‑1 (MQ=255)

                        cgcgGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGAt                                                                                                                               <  1:7383410/150‑1 (MQ=255)

                           ggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTGTGGCGCGGTCGGGTGCGGATGAg                                                                                                                            >  1:175022/1‑150 (MQ=255)

                            gcgcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGt                                                                                                                           <  2:6821674/150‑1 (MQ=255)

                            gcgcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGt                                                                                                                           <  1:8420204/150‑1 (MQ=255)

                              gcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTaa                                                                                                                         <  2:5758950/150‑1 (MQ=255)

                               cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTaaa                                                                                                                        <  1:3572633/150‑1 (MQ=255)

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

                                     ggTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAAATATGTCGCGGGCGGGTTAGGGTTAGTAAAAgtgtg                                                                                                                  >  1:3968201/1‑149 (MQ=255)

                                       ttGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCTGGCGCGGTCGGTTGCGGGTGAGTAAAGGTGTCGt                                                                                                                >  1:6139066/1‑150 (MQ=255)

                                         gCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGGTGAGTAAAGGTGTCGTgg                                                                                                              >  2:2061044/1‑150 (MQ=255)

                                         gCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTGCTGGCGCGGTCGGTTTCGGGTGGGTAAATGTGTCGTgg                                                                                                              >  2:8905691/1‑150 (MQ=255)

                                               gCTGCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAAc                                                                                                        <  2:5962087/150‑1 (MQ=255)

                                               gCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAAc                                                                                                        <  2:7964857/150‑1 (MQ=255)

                                                 tCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTTCGGGTGTGTAAAGGTGTCGTGGCATAACAg                                                                                                      >  1:5307117/1‑150 (MQ=255)

                                                 tCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTGTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTCTTCTGGCGCGGGCGGTTGCGGATGAGTAAAGGGGTCGTGGCATAACAg                                                                                                      >  1:6844846/1‑150 (MQ=255)

                                                   caccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTCTTGTGGCGCGGTCGGTTGCGGATGTGGAAAGGTGTCGGGGCGTAAAAgcg                                                                                                    >  1:2893144/1‑150 (MQ=255)

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

                                                    accaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGATTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAgcgc                                                                                                   <  2:3513829/150‑1 (MQ=255)

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

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

                                                      caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTATGGCGCGGTCGGTTGCGGATGTGTAAAGGTGTCGTGGCATAACAGCGCgg                                                                                                 >  1:6341252/1‑150 (MQ=255)

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

                                                       aGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGc                                                                                                <  2:7322426/150‑1 (MQ=255)

                                                         ggCCTGACGTCGACGCTGAAACTCTTGCCGGGATATCTCACTCATTACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTAGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGTGTAAAGGTGTCGTGGCATAACAGCGCGGCaa                                                                                              <  2:6344584/150‑1 (MQ=255)

                                                         ggCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTAATTCTGGCGCGGGCGGGTGCGGATGTGTAAAGGTGTGGTGGGAAAACAGCGCGGCaa                                                                                              >  1:4297397/1‑150 (MQ=255)

                                                           ccTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAaca                                                                                            <  1:1404553/150‑1 (MQ=255)

                                                           ccTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAaca                                                                                            <  1:7153844/150‑1 (MQ=255)

                                                            cTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAacac                                                                                           <  1:8276853/150‑1 (MQ=255)

                                                              gACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGa                                                                                         <  1:878581/150‑1 (MQ=255)

                                                              gACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGa                                                                                         <  2:8125044/150‑1 (MQ=255)

                                                               aCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGaa                                                                                        <  1:4097131/150‑1 (MQ=255)

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

                                                                       gCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTATGGCGGGGTCGGTTTCGGCTGAGTAAAGGTTTAGTTGAAAAACAACGCCGCCAAAACAACGTAttaa                                                                                >  1:675084/1‑148 (MQ=255)

                                                                       gCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACGTCTGGCGCGGTCGGTTGCGGTTGGGTAAAGGTGTCGTGGCCTAACAGCGCGGCAACACCAACCGATTCg                                                                                >  2:4668751/1‑150 (MQ=255)

                                                                         tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgan                                                                              <  1:4975689/150‑2 (MQ=255)

                                                                         tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTAATTATTTCGAGGTGGGTTGGGGTTGATTAAATGTGTGGTGGGATAAAAACGAGGCAACACGAAAGTATTAgat                                                                              >  2:7026570/1‑150 (MQ=255)

                                                                           aaaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTCTTCTGGCGCGGTCGGGTGGGGGGGTGTTAAGGTGTGGTGGGCTAAAACCGCGGCAACACCAAAGGGTTGgatga                                                                            >  1:197352/1‑150 (MQ=255)

                                                                             aCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGAtc                                                                          <  2:1863477/150‑1 (MQ=255)

                                                                             aCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGAtc                                                                          <  2:6509356/150‑1 (MQ=255)

                                                                               tctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTATTATTGGGCGGTGGGTTTCGGGTGAGTTAAGGTGTGGGGGCATAACAGCGCGGCAAAAACAACGGATTGgattatata                                                                        >  2:3187696/1‑147 (MQ=255)

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

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

                                                                                    ccGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCaa                                                                   <  1:1598705/150‑1 (MQ=255)

                                                                                       ggATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCg                                                                <  1:51721/150‑1 (MQ=255)

                                                                                        gATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGa                                                               <  2:606046/150‑1 (MQ=255)

                                                                                            tctcACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACGTATGGCGCGGTCGGTTGCGGATGAGTAAAGGGGTCGTGGCATAACACCGCGGCAACAACAACGTATTCGAAGATCTCTTCAAGCGACAtt                                                           >  2:8198505/1‑150 (MQ=255)

                                                                                             ctcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGCGGTCGGTTGCGGGTGGGTTAAGGTGTCGTGGCATAACAGCGCGGCCACACGCACGTATTCGATGATCTCTTCAAGCGCCAgtt                                                          >  2:3660314/1‑150 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGTGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGAAGATATATTCAAGCGACAttta                                                         >  2:3724805/1‑149 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       <  2:3927184/150‑1 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       <  1:7778575/150‑1 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATACGATGATCTCTTCAAGCGACATTTCaa                                                       <  2:2985939/150‑1 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACGTCAGGCGCGGGCGGTGGCGGCGGCGGAAAGGTGTCGGGGCCTAACAGCGCCGCAACACGAACGGAATCGATGATCACGTCAAGCCACAGTTCaa                                                       >  2:5424880/1‑150 (MQ=255)

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

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

                                                                                                  tcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGTCGCGGTCGGTTGCGGATGTGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTAATCGATGTTCTCTTCAAGCGACAGTTCAAGc                                                     >  2:6003559/1‑150 (MQ=255)

                                                                                                  tcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTACTTCTGGCGCGGGCGGTTGCGGGTGAGTAAAGGTGTCGTGGGCTAACAGCGCGGCCACAACCACCGATTAGAGGATATCTACAACCGACCAATCAAGc                                                     >  2:7180442/1‑150 (MQ=255)

                                                                                                   caTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGct                                                    <  1:7695904/150‑1 (MQ=255)

                                                                                                       acacTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGATTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTc                                                <  2:8687497/150‑1 (MQ=255)

                                                                                                           tctcCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGa                                            <  2:2762094/150‑1 (MQ=255)

                                                                                                              ccTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtct                                         <  2:3458885/150‑1 (MQ=255)

                                                                                                              ccTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtct                                         <  1:2167970/150‑1 (MQ=255)

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

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

                                                                                                                     ttttttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttcgtctt                                  <  2:222679/150‑1 (MQ=255)

                                                                                                                      tttttGTTTATAGTGTAGAGTCGGATTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttcgtcttc                                 <  2:7671183/150‑1 (MQ=255)

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

                                                                                                                             ttttAGTGTAGAGTCGGTTTTTGTAAGTATGGCGCGGTCGGTTTCGGATGATTAAATGTGTCGTGGCATAAAAGCGCGGCAAAACGAACGTATTAGATTTTATATTAAAGCGAAATTTAAAGCTCATACTGATCTTGGTCTTCGTCGTaa                          >  2:414485/1‑150 (MQ=255)

                                                                                                                                 aGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCa                      <  2:2612363/150‑1 (MQ=255)

                                                                                                                                       gagTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTgcg                <  2:3876345/150‑1 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCaa             <  1:3398802/150‑1 (MQ=255)

                                                                                                                                            ggTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATg           <  2:4425006/150‑1 (MQ=255)

                                                                                                                                            ggTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATg           >  2:2421876/1‑150 (MQ=255)

                                                                                                                                            ggTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATg           >  1:2283812/1‑150 (MQ=255)

                                                                                                                                              tttttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGtt         >  2:7442904/1‑150 (MQ=255)

                                                                                                                                                 ttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACg      >  1:189337/1‑150 (MQ=255)

                                                                                                                                                  tGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAAAAGAGCGGCACCACGAACGTATTCGATGATCTCTTCAAGAGACATTTCAAGATCTTCCAGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTaaga     >  1:7082479/1‑147 (MQ=255)

                                                                                                                                                   gTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCa                                                        >  2:1792261/1‑98 (MQ=255)

                                                                                                                                                     aCTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTAGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTAAATCTCTTCCTTATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATTTTAAGCAGa  >  1:801236/1‑150 (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: `‑`