BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,919,279	A→G	6.8%	V590V (GTA→GTG)	yebT →	MCE domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,919,279	0	A	G	6.8%	2123.0 / 95.4	643	V590V (GTA→GTG)	yebT	MCE domain protein
Reads supporting (aligned to +/- strand): ref base A (305/293); new base G (21/23); total (326/317)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 648 total aligned reads displayed.

GGTGAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCA  >  NC_000913/1919130‑1919428

ggTGAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTa                                                                                                                                                       <  2:9872207/150‑1 (MQ=255)

   gAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAg                                                                                                                                                    <  2:8243818/150‑1 (MQ=255)

   gAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAg                                                                                                                                                    <  1:543707/150‑1 (MQ=255)

     aGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGc                                                                                                                                                  >  2:7530399/1‑150 (MQ=255)

            aCCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCg                                                                                                                                           <  2:11802015/150‑1 (MQ=255)

             ccGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGc                                                                                                                                          <  2:5153069/150‑1 (MQ=255)

               gTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGctc                                                                                                                                        <  2:9108616/150‑1 (MQ=255)

                     ccGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAGGCATCCCCGCTCTCCaga                                                                                                                                  <  2:1689385/150‑1 (MQ=255)

                       gcgAGCTAACGCGTTTGATANCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGc                                                                                                                                >  1:10224138/1‑150 (MQ=255)

                         gAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCAt                                                                                                                              >  2:2596169/1‑150 (MQ=255)

                             tAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTaaa                                                                                                                          <  1:10877621/150‑1 (MQ=255)

                               aCGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAgg                                                                                                                        <  2:11998966/150‑1 (MQ=255)

                               aCGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAgg                                                                                                                        <  1:5114375/150‑1 (MQ=255)

                                 gcgTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGa                                                                                                                      <  1:8935523/150‑1 (MQ=255)

                                 gcgTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGa                                                                                                                      >  2:6025748/1‑150 (MQ=255)

                                  cgTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAg                                                                                                                     <  1:2736199/150‑1 (MQ=255)

                                    tttGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGcc                                                                                                                   >  1:11735200/1‑150 (MQ=255)

                                          atCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGc                                                                                                             >  1:10607015/1‑150 (MQ=255)

                                          atCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGACGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGc                                                                                                             <  1:11868694/150‑1 (MQ=255)

                                             gATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCTTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGTGCATTAAAGGGAGCCATTAGCTTc                                                                                                          <  2:11413343/150‑1 (MQ=255)

                                             gATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTc                                                                                                          <  1:2209157/150‑1 (MQ=255)

                                             gATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTc                                                                                                          <  1:12249179/150‑1 (MQ=255)

                                             gATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTc                                                                                                          <  2:924003/150‑1 (MQ=255)

                                              aTCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCg                                                                                                         >  2:6180431/1‑150 (MQ=255)

                                               tCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGa                                                                                                        <  2:10575866/150‑1 (MQ=255)

                                               tCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGa                                                                                                        <  2:685687/150‑1 (MQ=255)

                                                 tGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACa                                                                                                      >  2:5618383/1‑150 (MQ=255)

                                                  gCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACaa                                                                                                     >  1:9841975/1‑150 (MQ=255)

                                                  gCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACaa                                                                                                     >  1:3238277/1‑150 (MQ=255)

                                                  gCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACaa                                                                                                     >  2:2813637/1‑150 (MQ=255)

                                                  gCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACaa                                                                                                     >  2:3592518/1‑150 (MQ=255)

                                                      atTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTc                                                                                                 >  1:1959105/1‑150 (MQ=255)

                                                       ttAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCa                                                                                                >  2:3159613/1‑150 (MQ=255)

                                                        tAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAg                                                                                               >  2:11552049/1‑150 (MQ=255)

                                                        tAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAg                                                                                               >  2:9680240/1‑150 (MQ=255)

                                                         aaGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCCGAGCATTAAAGGGAGCCATTAGCTGCGACAACCTCAGc                                                                                              >  2:11276886/1‑150 (MQ=255)

                                                          aGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCg                                                                                             <  1:4899032/150‑1 (MQ=255)

                                                            ccGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGt                                                                                           >  2:2157056/1‑150 (MQ=255)

                                                            ccGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGt                                                                                           >  1:2092868/1‑150 (MQ=255)

                                                            ccGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGt                                                                                           >  2:4248136/1‑150 (MQ=255)

                                                             cGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTg                                                                                          >  2:11020219/1‑150 (MQ=255)

                                                             cGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTg                                                                                          <  2:3082524/150‑1 (MQ=255)

                                                               gAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGcc                                                                                        <  2:8714904/150‑1 (MQ=255)

                                                                    tCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAgcgc                                                                                   <  1:7412677/150‑1 (MQ=255)

                                                                    tCGCAACCTTCTGACCAGCAAGAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAgcgc                                                                                   <  1:100004/150‑1 (MQ=255)

                                                                      gCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCa                                                                                 <  1:3020995/150‑1 (MQ=255)

                                                                      gCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGAATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCa                                                                                 <  1:3658083/150‑1 (MQ=255)

                                                                       cAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGACAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAg                                                                                >  2:7976389/1‑150 (MQ=255)

                                                                       cAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAg                                                                                >  2:4550087/1‑150 (MQ=255)

                                                                          ccTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCa                                                                             <  1:1124389/150‑1 (MQ=255)

                                                                            ttCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGc                                                                           >  2:8143132/1‑150 (MQ=255)

                                                                                gACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTaa                                                                       <  2:3768011/150‑1 (MQ=255)

                                                                                gACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTaa                                                                       >  1:463332/1‑150 (MQ=255)

                                                                                  ccAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAg                                                                     <  1:8370777/150‑1 (MQ=255)

                                                                                  ccAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAg                                                                     <  1:10451229/150‑1 (MQ=255)

                                                                                  ccAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTAAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAg                                                                     >  1:1507656/1‑150 (MQ=255)

                                                                                   cAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGTCAGCGCCAGTCAGCGTAAAgg                                                                    >  2:5823644/1‑150 (MQ=255)

                                                                                   cAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAgg                                                                    >  1:2504806/1‑150 (MQ=255)

                                                                                   cAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAgg                                                                    >  1:6363677/1‑150 (MQ=255)

                                                                                   cAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAgg                                                                    >  1:954621/1‑150 (MQ=255)

                                                                                   cAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAgg                                                                    <  2:5037226/150‑1 (MQ=255)

                                                                                    aGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGc                                                                   >  1:926899/1‑150 (MQ=255)

                                                                                    aGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGc                                                                   <  2:8433195/150‑1 (MQ=255)

                                                                                     gCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCg                                                                  >  1:12007942/1‑150 (MQ=255)

                                                                                     gCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCg                                                                  <  2:8875917/150‑1 (MQ=255)

                                                                                     gCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCg                                                                  <  1:5288610/150‑1 (MQ=255)

                                                                                       aaTAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGAt                                                                <  1:11542496/150‑1 (MQ=255)

                                                                                         tAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATaa                                                              <  2:4807435/150‑1 (MQ=255)

                                                                                             gtgtTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGt                                                          >  2:7909494/1‑150 (MQ=255)

                                                                                               gtTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTAt                                                        >  1:7833811/1‑150 (MQ=255)

                                                                                               gtTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGAGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTAt                                                        >  2:3910520/1‑150 (MQ=255)

                                                                                                ttCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTAtt                                                       <  2:10797690/150‑1 (MQ=255)

                                                                                                               ggggCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGaa                                        >  1:6806709/1‑150 (MQ=255)

                                                                                                                 ggCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCtctc                                                                                                                                      <  1:10501570/54‑1 (MQ=255)

                                                                                                                    gAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGc                                   >  2:10470482/1‑150 (MQ=255)

                                                                                                                     aaaGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCg                                  >  2:2345550/1‑150 (MQ=255)

                                                                                                                     aaaGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCg                                  >  2:7309181/1‑150 (MQ=255)

                                                                                                                      aaGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCgg                                 <  1:4264163/150‑1 (MQ=255)

                                                                                                                       aGTTCAGCTGATTGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGc                                <  1:9181635/150‑1 (MQ=255)

                                                                                                                        gTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGcc                               >  1:8641441/1‑150 (MQ=255)

                                                                                                                          tCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCg                             >  2:10075250/1‑150 (MQ=255)

                                                                                                                            aGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTg                           >  1:8239777/1‑150 (MQ=255)

                                                                                                                                    ggTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGc                   >  1:5806249/1‑150 (MQ=255)

                                                                                                                                     gTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCg                  <  1:6591432/150‑1 (MQ=255)

                                                                                                                                      tAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCgg                 >  1:716929/1‑150 (MQ=255)

                                                                                                                                       aGTGGTCTGACCGAGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCggg                >  1:54749/1‑150 (MQ=255)

                                                                                                                                        gTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGc               <  2:1385277/150‑1 (MQ=255)

                                                                                                                                         tGGTCTGACCGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCa              <  1:1633947/150‑1 (MQ=255)

                                                                                                                                         tGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCa              <  2:5463364/150‑1 (MQ=255)

                                                                                                                                         tGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCa              <  2:4062468/150‑1 (MQ=255)

                                                                                                                                         tGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCa              <  1:11020219/150‑1 (MQ=255)

                                                                                                                                           gTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGa            <  1:1690809/150‑1 (MQ=255)

                                                                                                                                             cTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGAtt          <  1:3901007/150‑1 (MQ=255)

                                                                                                                                              tGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTa         <  1:6153452/150‑1 (MQ=255)

                                                                                                                                               gACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTAc        >  1:805915/1‑150 (MQ=255)

                                                                                                                                                aCCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGTGGGCCGATTACg       >  2:7109410/1‑150 (MQ=255)

                                                                                                                                                aCCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACg       >  2:3311069/1‑150 (MQ=255)

                                                                                                                                                aCCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACg       <  1:5952831/150‑1 (MQ=255)

                                                                                                                                                aCCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACg       >  2:9460362/1‑150 (MQ=255)

                                                                                                                                                 ccGTGCAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGc      <  2:6402553/150‑1 (MQ=255)

GGTGAAGTGATTACCGTGCGTCCGCGAGCTAACGCGTTTGATATCGATCTGCATATTAAGCCGGAGTATCGCAACCTTCTGACCAGCAATAGCGTGTTCTGGGCAGAAGGCGGGGCGAAAGTTCAGCTGAATGGTAGTGGTCTGACCGTACAGGCATCCCCGCTCTCCAGAGCATTAAAGGGAGCCATTAGCTTCGACAACCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCA  >  NC_000913/1919130‑1919428

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