BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,125	T→G	11.0%	Y416S (TAC→TCC)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,125	0	T	G	11.0%	394.5 / 12.0	155	Y416S (TAC→TCC)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (37/98); new base G (18/2); total (55/100)
Fisher's exact test for biased strand distribution p-value = 1.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

AGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727980‑3728274

 gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGACGGtcctcg                                                                                                                                                  >  2:958807/1‑145 (MQ=255)

   ggCATCACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCACTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAg                                                                                                                                                <  2:2545818/150‑1 (MQ=255)

   ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAg                                                                                                                                                >  2:1119371/1‑150 (MQ=255)

       tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGGACTCGCgcggc                                                                                                                                            >  2:2689812/1‑150 (MQ=255)

         aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGctgctg                                                                                                                                          <  1:2499492/150‑1 (MQ=255)

         aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTCAGCGAGctgctg                                                                                                                                          >  2:914357/1‑150 (MQ=255)

              gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGc                                                                                                                                     <  2:852558/150‑1 (MQ=255)

               cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                    <  1:2662917/150‑1 (MQ=255)

               cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                    <  1:2730051/150‑1 (MQ=255)

               cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                    <  1:773232/150‑1 (MQ=255)

                   tCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGata                                                                                                                                <  1:1528754/150‑1 (MQ=255)

                     tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGTCGGTGATCGGGCTGCTGGCCGGTggtgtc                                                                                                                              >  2:266005/1‑145 (MQ=255)

                       gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGGGAGCGCGCGGCTGACCGCGGaggggcg                                                                                                                            >  2:2678279/1‑144 (MQ=255)

                       gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGTCCGGGAATCGAGCTGCTGCCCGCGGGTATCCg                                                                                                                            >  2:2678087/1‑150 (MQ=255)

                         aGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTTCGGTAATCGAGCTGCTGACCGCTGATATCCGcc                                                                                                                          <  1:232389/150‑1 (MQ=255)

                         aGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGcc                                                                                                                          <  2:1044958/150‑1 (MQ=255)

                               gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATc                                                                                                                    <  1:14309/150‑1 (MQ=255)

                                aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTGGTCGCGCTGCTGACCGCTGATAGccgcccgccgct                                                                                                                   >  2:2181166/1‑146 (MQ=255)

                                  tGTTCTATCGGTAGTTGCGGCAACAATTCAATGAGCGAGTTCTCTGGCTTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                                 <  1:473059/150‑1 (MQ=255)

                                  tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                                 <  1:1339828/150‑1 (MQ=255)

                                  tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                                 <  1:346709/150‑1 (MQ=255)

                                  tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                                 <  1:1339823/150‑1 (MQ=255)

                                       tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTACTCGGGCTGCTGACCGCTGAGATCCGCCAGCGTCTGACGccc                                                                                                            >  2:1935823/1‑149 (MQ=255)

                                         aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGCCGGGCATCGAGCTGCTGGCCGCGGATAGCCGCCCGCAGCTGACGccggg                                                                                                          >  2:2746193/1‑147 (MQ=255)

                                          cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                                                         <  2:2531793/150‑1 (MQ=255)

                                            gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACt                                                                                                       <  1:1611458/150‑1 (MQ=255)

                                            gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTGCGGGAGTCGAGCTGCTGGCCGCTGCTCTCCGCCGGCCTCTGACGCCCGGACt                                                                                                       >  2:514912/1‑150 (MQ=255)

                                             tAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTc                                                                                                      <  1:1472007/150‑1 (MQ=255)

                                                   cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcg                                                                                                <  1:1614064/150‑1 (MQ=255)

                                                   cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcg                                                                                                <  1:842397/150‑1 (MQ=255)

                                                    ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTGCGGGCAGCGAGCTGCTGACCGCGGGTATCCGCCGGCGGCTGGCGCCCGTCCTCACGAcgc                                                                                               >  2:2359664/1‑150 (MQ=255)

                                                      caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTACTCGCGCTGCTGTCCGCTGGTCTGCGCCCGCGGCTGACGGCAGTACTCGCTACgcgc                                                                                             >  2:1623518/1‑150 (MQ=255)

                                                      caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCGGATATCCGCCAGCATCTGACGCCAGTACTCGCTCCgcgc                                                                                             >  2:120936/1‑150 (MQ=255)

                                                      caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACgcgc                                                                                             <  1:1182985/150‑1 (MQ=255)

                                                           aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGGGCGGGGATCGGGCTGCTGGGCGCTGGGAGCCGCCGGCCGCTGGCGCCGGGGGTCGCTGCGcgccccgc                                                                                        >  1:2443061/1‑150 (MQ=255)

                                                           aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCGGTACTCACTACGcgccccgc                                                                                        >  2:1756184/1‑150 (MQ=255)

                                                           aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTACGGTAATCGGGCTGCTGACCGCTGATATCCGCCCGCATCTGTCGCCAGTACTCCCTACGcgccccgc                                                                                        >  2:1291045/1‑150 (MQ=255)

                                                              cAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCgccccgcccc                                                                                     <  2:2460916/150‑1 (MQ=255)

                                                                  gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCCGGGAGCGCGCTGCGGGCCGCGGAGAGCCGCCAGCAGCAGCCGCCCGGCCGCACGACGCGcccccccccccac                                                                                 >  2:1733449/1‑147 (MQ=255)

                                                                   aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGTGCGGTGCTGGGGCTGGGGACCGCGGTTGTCCGCCTGCGTCTGACGCCCGTGCGCACTCGGCGCCCCGgccccgctc                                                                                >  1:1251323/1‑146 (MQ=255)

                                                                   aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCa                                                                               >  1:2110766/1‑150 (MQ=255)

                                                                      gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                             <  2:2167058/150‑1 (MQ=255)

                                                                      gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                             <  2:2230745/150‑1 (MQ=255)

                                                                      gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                             <  2:2386945/150‑1 (MQ=255)

                                                                       aTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTGATCGAGCTGCTGGCCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTCCGCGCCCCGCCCCCCCTCAACg                                                                            >  2:2911001/1‑150 (MQ=255)

                                                                          tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTaa                                                                         <  1:582241/150‑1 (MQ=255)

                                                                             ctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCTGTACTCACTACGCGCCCCGCCCCCAATCAACGTAacac                                                                      <  2:2609723/150‑1 (MQ=255)

                                                                                  ttCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGAATCGGGCTGCTGTCCGCGGATATCCGCCAGCGTCTGACGCCAGTCCTCACTGCGCGCCCCGCCCCCAATCCACGTAGCACTGTGc                                                                 >  2:386525/1‑150 (MQ=255)

                                                                                   tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCg                                                                <  1:16071/150‑1 (MQ=255)

                                                                                    cgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCgg                                                               <  2:2724726/150‑1 (MQ=255)

                                                                                    cgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCgg                                                               <  2:1023540/150‑1 (MQ=255)

                                                                                      ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACCCTTTGCGGtt                                                             <  2:406735/150‑1 (MQ=255)

                                                                                      ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                             <  2:1505388/150‑1 (MQ=255)

                                                                                      ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                             <  2:57076/150‑1 (MQ=255)

                                                                                      ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGGAATCGAGCTGCTGACCGCTGATATCCGCCCGCAGCGGACGCCAGTACGCACTCCGCGCCCCGCCCCCCCTCCACGTAACCCGTTGCgggg                                                             >  2:756518/1‑148 (MQ=255)

                                                                                       cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                            <  1:2007577/150‑1 (MQ=255)

                                                                                       cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                            <  1:1605849/150‑1 (MQ=255)

                                                                                       cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                            <  1:476839/150‑1 (MQ=255)

                                                                                       cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                            <  2:732175/150‑1 (MQ=255)

                                                                                           aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGACTCGGGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCCGTACTCCCTACGCGCCCCGCCCCCATTCAACGTACCACGTTGCGGttttcta                                                        >  2:2343318/1‑146 (MQ=255)

                                                                                                cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCa                                                   <  1:262350/150‑1 (MQ=255)

                                                                                                cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAGTCGAGCTGCTGACCGCTGATATCCGCCAGCGTCTGACGCCAGTACGCGCTTCGCGCCCCGCCCCCACTCACCGGGCCCCGTTGCGGTTTAGTACCGct                                                   >  2:2501637/1‑149 (MQ=255)

                                                                                                cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                   <  1:2778820/150‑1 (MQ=255)

                                                                                                cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                   <  1:1751121/150‑1 (MQ=255)

                                                                                                cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                   <  2:578786/150‑1 (MQ=255)

                                                                                                 gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaa                                                  <  1:1731960/150‑2 (MQ=255)

                                                                                                  ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAgg                                                 <  1:2689812/150‑1 (MQ=255)

                                                                                                         tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc                                          <  2:1629965/150‑1 (MQ=255)

                                                                                                         tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc                                          <  1:2851885/150‑1 (MQ=255)

                                                                                                         tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc                                          <  2:2389889/150‑1 (MQ=255)

                                                                                                         tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc                                          <  2:210304/150‑1 (MQ=255)

                                                                                                           ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGa                                        <  1:1504883/150‑1 (MQ=255)

                                                                                                            tgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGAc                                       <  1:1201234/150‑1 (MQ=255)

                                                                                                               gcCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGAGCTGCTGGCCGCTGATGTCCGCCCGCAGCTGACGCCCGTACTCACTACGCGCCCCGCCCCCCATCAACGTGACACGTTGCGGTTTAATACCGCAGGCATGCACGACATc                                    >  1:1714302/1‑150 (MQ=255)

                                                                                                                 ccAGTGCTGGCCCCACATCCCCCCCCGGACGGTAATCGAGCTGCTGACCGCTGGTATCCGCCAGCATCTGTCGCCAGTACTCACTACGCGCCCCGCCCCCCCTCAACGTCACCCTTTGCGGGTTAAGACCGCCGGCATGCGCGACATCCa                                  >  1:2900813/1‑150 (MQ=255)

                                                                                                                  cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                 <  2:304520/150‑1 (MQ=255)

                                                                                                                  cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                 <  2:1110061/150‑1 (MQ=255)

                                                                                                                  cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                 <  1:1192762/150‑1 (MQ=255)

                                                                                                                  cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                 <  1:2723483/150‑1 (MQ=255)

                                                                                                                    gTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAGCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc                               >  1:95971/1‑150 (MQ=255)

                                                                                                                      gCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCa                             <  1:1098536/150‑1 (MQ=255)

                                                                                                                         ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                          <  2:1570273/150‑1 (MQ=255)

                                                                                                                         ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                          <  2:1570274/150‑1 (MQ=255)

                                                                                                                         ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                          <  2:1377787/150‑1 (MQ=255)

                                                                                                                          gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGTCCGCGGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGGAACACTTTGCGGGTTAATACCGCAGGCATGCCCGACATCCATGCCATCTg                         >  1:2228284/1‑150 (MQ=255)

                                                                                                                          gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                         >  2:2084858/1‑150 (MQ=255)

                                                                                                                           ccccACATCCCCCCCCGTACGGTAAACGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACAACGCGCCCCGCCCCCAATCAACGTAACACATTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGc                        >  2:1247573/1‑150 (MQ=255)

                                                                                                                             ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCCTCTGCCa                      >  1:16712/1‑150 (MQ=255)

                                                                                                                             ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGccc                      >  2:2300908/1‑149 (MQ=255)

                                                                                                                             ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                      <  2:2814478/150‑1 (MQ=255)

                                                                                                                             ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                      <  1:1306560/150‑1 (MQ=255)

                                                                                                                             ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                      <  1:1478731/150‑1 (MQ=255)

                                                                                                                                  tCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCa                 <  1:199710/150‑1 (MQ=255)

                                                                                                                                  tCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATACGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCa                 >  1:1114565/1‑150 (MQ=255)

                                                                                                                                      cccccGTACGGTAATCTAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTAAATACGAGCACCGCCCCCAATCAACGTAACACTTTGCGTTTTAATACCGCAGGCATGCACGACATCCATTCCATCTGACAGCGCATAAc             >  2:441461/1‑150 (MQ=255)

                                                                                                                                      cccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAc             >  2:441466/1‑150 (MQ=255)

                                                                                                                                          cGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc         <  1:1874019/150‑1 (MQ=255)

                                                                                                                                          cGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc         <  1:2359664/150‑1 (MQ=255)

                                                                                                                                            tACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGc       >  1:1720479/1‑150 (MQ=255)

                                                                                                                                               ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCtt    >  1:40590/1‑150 (MQ=255)

AGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727980‑3728274

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