BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,054,888	A→T	6.2%	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	6.2%	788.6 / 2.5	247	V186E (GTA→GAA)	ygfB	UPF0149 family protein
Reads supporting (aligned to +/- strand): ref base A (89/137); new base T (16/3); total (107/140)
Fisher's exact test for biased strand distribution p-value = 1.72e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037

tACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTa                                                                                                                                                       >  2:3678468/1‑150 (MQ=255)

tACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTATAGAGTCGGTTTTTGTa                                                                                                                                                       <  1:3647050/150‑1 (MQ=255)

tACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCAGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGATTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTa                                                                                                                                                       <  2:2598056/150‑1 (MQ=255)

 aCTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTGTTGTAc                                                                                                                                                      >  1:2396798/1‑150 (MQ=255)

   ttCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACtt                                                                                                                                                    >  2:1497590/1‑150 (MQ=255)

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

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

         tGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACGTCTGgcg                                                                                                                                              >  2:2533607/1‑150 (MQ=255)

          gcagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGgcgc                                                                                                                                             <  1:2711251/150‑1 (MQ=255)

           cagcaAAAATCAGCGCGGCGCTGCCGGGGTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGcgcg                                                                                                                                            <  2:1469763/150‑1 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTAAATGTGGCGcggtg                                                                                                                                         >  1:1656671/1‑149 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTCTTATGGCGcgggg                                                                                                                                         >  1:1771207/1‑148 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATGTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCTGGCGcggtc                                                                                                                                         >  2:2349228/1‑150 (MQ=255)

                 aaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTTCGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtcggt                                                                                                                                      <  2:3001443/150‑1 (MQ=255)

                   aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTg                                                                                                                                    <  2:606555/150‑1 (MQ=255)

                    tCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGc                                                                                                                                   <  2:2668044/150‑1 (MQ=255)

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

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

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

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

                       gcgcGGCGCTGCCGGGTTGTATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTACTTCGGGCGCGGGCGGGTGCggc                                                                                                                                >  2:2896448/1‑149 (MQ=255)

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

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

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

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

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

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

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

                           ggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTATGGCGCGGTCGGGTGCGGGTGAg                                                                                                                            >  2:3963410/1‑150 (MQ=255)

                              gcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAGTTCTGGCGCGGGCGGTTGCGGGTGGGTaa                                                                                                                         >  2:1234688/1‑150 (MQ=255)

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

                                  ccGGGGTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGgt                                                                                                                     <  2:360551/150‑1 (MQ=255)

                                   cGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTAATTATGGCGCGGGCGGGTGGGGATGTGTTAAGgtg                                                                                                                    >  1:1287015/1‑150 (MQ=255)

                                    gggTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACATCTGGCGCGGTCGGGTTCGGTTGTGTAAAGgtgt                                                                                                                   >  1:2973135/1‑150 (MQ=255)

                                    gggTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTGTTGTTTTTAGTGTAGAGTCGGTGTTTGTACTTCTGGCGCGGGCGGGTGCGGATGCGGAAAGgtgt                                                                                                                   >  2:1256083/1‑150 (MQ=255)

                                     ggTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTc                                                                                                                  <  2:3316039/150‑1 (MQ=255)

                                      gTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGGTGTGTAAAGGTGTCg                                                                                                                 >  2:1292847/1‑150 (MQ=255)

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

                                        tGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTg                                                                                                               <  2:2757029/150‑1 (MQ=255)

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

                                         gCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGATTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTgg                                                                                                              <  2:323243/150‑1 (MQ=255)

                                           aTTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTATTAGTGCGCGGTGGGTTGCGGGTGGGTGAAAGTGTGGTGGCa                                                                                                            >  2:1404415/1‑150 (MQ=255)

                                            tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTTCTTCTGGCGCGGGGGGTTTCGGGTGGGTAAAGGTGTCGTGGGAt                                                                                                           >  1:3506059/1‑150 (MQ=255)

                                                cTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTACTTATGGCGCGGTCGGTTGCGGGTGTGTAAAGGTGTCGTTGCATaaaa                                                                                                       >  1:1702436/1‑148 (MQ=255)

                                                 tCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGGCGGTTGCGGCTGTGTAAAGGTGTCGTGGCATAAAAg                                                                                                      >  2:3168377/1‑150 (MQ=255)

                                                 tCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTATTCTTGCGCGGTCGGTTGCGGATTAGGAAAGGTGGCGTGGCATAAAAg                                                                                                      >  1:2486006/1‑150 (MQ=255)

                                                  ccaccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTAATATGGCGCGGGCGGGTGGGGGTGTGTAAAGGTGTCGGGGGAAAACAgc                                                                                                     >  1:1461588/1‑150 (MQ=255)

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

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

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

                                                      caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTGTTTTTACTTATGGCGCGGTCGGGTGCGGATGTGTAAAGGTGTCGTGGCATAAAAGCGCgg                                                                                                 >  2:534919/1‑150 (MQ=255)

                                                             tGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACg                                                                                          >  2:2380417/1‑150 (MQ=255)

                                                               aCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGGTGGGGTTGAGTAAAGGTGTGGTGGCAAAACAGCGCGGCAACACGaa                                                                                        >  1:4051633/1‑150 (MQ=255)

                                                                     aCGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTAtt                                                                                  <  2:2091676/150‑1 (MQ=255)

                                                                       gCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTAATTATTGGGCGGTGGGGTGCGGGTGAGTGAAGGTGTCGTGGCATAAAAGCGCGGCAACACGAACggaatag                                                                                >  1:1389670/1‑144 (MQ=255)

                                                                        cTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTACGTATGTCGCGGTCGGGTGCGGGTGGGGAGAGGGGTCGGGGCATAACGGCGCGGCAACACGGACGTGTTCgc                                                                               >  2:1695453/1‑149 (MQ=255)

                                                                              ctctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTATGGCGCGGTCGGTTTCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGAtct                                                                         >  2:3098044/1‑150 (MQ=255)

                                                                               tctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCTGGCGCGGTCGGTTGCGGGTGAGTAAAGGTGTGGTGGCCTAACAGCGCGGCAAAACGAACGTATTAGATGAtatc                                                                        >  2:1923262/1‑150 (MQ=255)

                                                                                ctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTGAAGGTGTCGTGGCCTAAAAGCGCGGCAACACGGACGTATTCGATGATctct                                                                       >  2:2428970/1‑150 (MQ=255)

                                                                                 ttGCCGGGATATCTCACTCATAACACTCTCCTTCCGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCtt                                                                      <  2:2815119/150‑1 (MQ=255)

                                                                                  tGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTc                                                                     <  2:2878426/150‑1 (MQ=255)

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

                                                                                       ggATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGCGGTCGGTTGCGGGTGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTATTCAAGCg                                                                >  2:3785868/1‑150 (MQ=255)

                                                                                        gATATCTCACTCATCACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTGTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGa                                                               <  2:777022/150‑1 (MQ=255)

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

                                                                                            tctcACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAAAAGCGCGGCCACACGAAAGTATTAGATGTTCTCTTCAAGCGACaat                                                           >  2:2738249/1‑148 (MQ=255)

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

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTGCATCATGCGCGGTCGGTTGCGGATGGGTAAAGGTGTCGTGGCAAAAAACCGCGGCAACACGAACGGATTAGAAGATATCTTAAAGCGAccatta                                                         >  2:926701/1‑145 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAAATATGGCGCGGTGGGTTTCGGATGAGTAAAGGTGTCGTGGCAAAAAAGAGCGGCAACACGAAAGTAATCGATGATATATTCAAGCGACAttta                                                         >  2:847588/1‑149 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAAATAATTCGCCGTGGGTTGCGGATGAGTAAAAGTGTCGTGGCAAAAAAAACCGGCAACACCAAACTATTAGAATATCTAATAAAACGACAatta                                                         >  2:122205/1‑149 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACGTCTGGGGCGGTCGGGTGCGGCGGTGTCCAGGTGTCGTGGCATAACAGCGCGGCAACACGCACGGATTCGATGATCACTTCAAGCGACATTTc                                                         >  2:2255957/1‑150 (MQ=255)

                                                                                               cactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTCTGGGGCGGGCGGTTTCGGATGTGTAAAGGTTTCGTGGCATAAAAGCGCCGCAAAACCCACCGATTCGATTATCTCTTAAAGCGAAAtttaa                                                        >  2:3927033/1‑148 (MQ=255)

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

                                                                                                       acacTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTc                                                <  1:1198678/150‑1 (MQ=255)

                                                                                                       acacTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTc                                                <  1:2733020/150‑1 (MQ=255)

                                                                                                                   cGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttcgtc                                    <  1:1073657/150‑1 (MQ=255)

                                                                                                                     ttttttGTTTTTAGTGTAGAGTCGGTTTTTGTACATAAGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCAAAACAGCGCCGCAACAACAACGTATTACAAGATCACATAAAGCGACAATTCAAGCACATACAGAtattagtctt                                  >  2:2732361/1‑150 (MQ=255)

                                                                                                                      tttttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttcgtcttc                                 <  1:3486844/150‑1 (MQ=255)

                                                                                                                       ttttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttcgtcttcg                                <  1:1676101/150‑1 (MQ=255)

                                                                                                                        ttttttttTAGTGTAGAGTCGGTTTTTGTTCTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGTCATAACAGCGCGGCAACACGAACGTATTCGATGATCTGTTCAAGCGACATTTCAAGCTCTTCCTGATcttcgtcttcgt                               <  2:2430049/150‑1 (MQ=255)

                                                                                                                        tttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATcttcgtcttcgt                               <  2:3258060/150‑1 (MQ=255)

                                                                                                                          tGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTtcgtcg                             <  1:3871212/150‑1 (MQ=255)

                                                                                                                           gTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTcgtcgt                            <  1:1416342/150‑1 (MQ=255)

                                                                                                                            tttttAGTGTAGATACGGGTTTTGTAAGTATGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAAAAGCGCGGCCACACGAACGTATTCGGTGTTGGCTTCAAGCGACATTTCAAGCGCTTGCTGATCTTAGtcgtcgtcgga                           >  2:4094512/1‑148 (MQ=255)

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

                                                                                                                                tAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAAccc                       <  2:2135767/150‑1 (MQ=255)

                                                                                                                                  gtgtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAg                     <  2:1767475/150‑1 (MQ=255)

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

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

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCaa             >  2:1135345/1‑150 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTAAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCaa             >  1:2201829/1‑150 (MQ=255)

                                                                                                                                           cGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGGTCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAAt            <  2:3503278/150‑1 (MQ=255)

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

                                                                                                                                             gTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACAATTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGt          >  2:3117108/1‑150 (MQ=255)

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

                                                                                                                                               ttttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTa        <  1:1753100/150‑1 (MQ=255)

                                                                                                                                                tttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTAc       >  2:677743/1‑150 (MQ=255)

                                                                                                                                                tttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTAc       >  2:2682737/1‑150 (MQ=255)

                                                                                                                                                 ttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACg      <  1:3685272/150‑1 (MQ=255)

                                                                                                                                                   gTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGTTCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCCGTTGCGCAATGTTACGCa    >  1:1504649/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: `‑`