BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,127	A→C	7.1%	D415E (GAT→GAG)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,127	0	A	C	7.1%	1375.7 / 4.6	482	D415E (GAT→GAG)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base A (157/268); new base C (41/3); total (211/271)
Fisher's exact test for biased strand distribution p-value = 1.08e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 497 total aligned reads displayed.

AAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAG  >  NC_000913/3727979‑3728276

aaGCGGCATCACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCCACAATTCCATGCGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCACCCCCGTAAGGTAAt                                                                                                                                                      <  2:3200653/150‑1 (MQ=255)

 aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATc                                                                                                                                                     <  1:8392393/150‑1 (MQ=255)

  gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCg                                                                                                                                                    <  2:7590229/150‑1 (MQ=255)

   cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                                   >  2:9868722/1‑150 (MQ=255)

   cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                                   <  1:9963237/150‑1 (MQ=255)

    ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAg                                                                                                                                                  <  1:8162554/150‑1 (MQ=255)

     gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCCATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgc                                                                                                                                                 <  2:5029873/150‑1 (MQ=255)

     gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGGCGGGACTCGAgg                                                                                                                                                 >  2:4023927/1‑149 (MQ=255)

        tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgctgc                                                                                                                                              <  2:2450700/150‑1 (MQ=255)

        tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAAGTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgctgc                                                                                                                                              <  2:5759458/150‑1 (MQ=255)

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

           cGCTGCGCATCTGTTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGACCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGCTAATCGAGCTGCTGa                                                                                                                                           <  2:4260679/150‑1 (MQ=255)

           cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCATCCCCCCCCGGCCGGGGATCGCGctgctgg                                                                                                                                           >  2:2338686/1‑149 (MQ=255)

           cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGGCCCCCCCCGTACGGTACTCGGGCTGCTGa                                                                                                                                           >  2:3846744/1‑150 (MQ=255)

           cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                                                                           <  2:2654489/150‑1 (MQ=255)

           cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                                                                           <  2:8368793/150‑1 (MQ=255)

            gCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGAGCTGCTGAc                                                                                                                                          >  2:4805648/1‑150 (MQ=255)

            gCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGAc                                                                                                                                          >  2:1492647/1‑150 (MQ=255)

             cTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGAcc                                                                                                                                         <  2:3580603/150‑1 (MQ=255)

              tGCGCATCTGGTAGATGCGACTGTTCTAACGGGAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCCGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCg                                                                                                                                        <  2:6439306/150‑1 (MQ=255)

                cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACACTTCAATGAGCGAGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                      <  2:6333120/150‑1 (MQ=255)

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

                 gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGTGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGACGCGAGCTGCGGAgcgcgg                                                                                                                                     >  2:5452126/1‑148 (MQ=255)

                 gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                                     <  1:1982606/150‑1 (MQ=255)

                 gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                                     <  2:3807803/150‑1 (MQ=255)

                       ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCTGCTGTCCGCTggtgtcg                                                                                                                               >  2:5230670/1‑144 (MQ=255)

                       ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                               <  1:10088496/150‑1 (MQ=255)

                         tAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCCGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGTCCGGGCCGCGAGCGGCTGCCCGCGGAGCTCCGc                                                                                                                             >  2:9581165/1‑150 (MQ=255)

                           gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                                           <  1:736418/150‑1 (MQ=255)

                           gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                                           <  1:1821232/150‑1 (MQ=255)

                               cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTAATCGGGCTGCTGACCGCTGATATccgcccgccg                                                                                                                       >  2:9206819/1‑148 (MQ=255)

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

                                   tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGACGGGGCGCGGGCTGCTGACCGCTGAGATCCGCCGGCCTCTgc                                                                                                                   >  2:10309211/1‑149 (MQ=255)

                                     ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGTCCGCTGATATCCGCCCGCAGCTGTCg                                                                                                                 >  2:274984/1‑150 (MQ=255)

                                      tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                                <  1:6505720/150‑1 (MQ=255)

                                      tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                                <  2:4428624/150‑1 (MQ=255)

                                        tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTAATCGAGCTGCTGAGCGCTGATGGGCGGCGGCTTCTGACggcg                                                                                                              >  2:10999490/1‑147 (MQ=255)

                                          aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGt                                                                                                            >  1:5198147/1‑150 (MQ=255)

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

                                                            aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGTCGGTTGGCGGGCGGCTGTCCGGGGGTGGGCGCCAGCGTCGGGGGCCAGTTCTCACTGcgcgccgcgc                                                                                          >  1:2746604/1‑150 (MQ=255)

                                                                aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCGGGTGTCCGCTGATATCCGCCAGCCTCTGGCGCCAGTCCTCACTGcgcgccccgcgccc                                                                                      >  2:894246/1‑150 (MQ=255)

                                                                 aTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTCCGGGGCGGGTGCGGCGGGCCGCGGATATCCGCCCGCATCTGGCGCCGGTACTCACTCcgcgccccgcgcccc                                                                                     >  1:3624708/1‑149 (MQ=255)

                                                                 aTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCa                                                                                     <  1:933480/150‑1 (MQ=255)

                                                                   gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                                                                   <  1:1100545/150‑1 (MQ=255)

                                                                    aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCCGCCCCCCCCGGCCGGGCCGCGAGCTGCGGCCCGCGGGGCTCCGCCAGCCGCGGACGCCCGGCCTCACGCCGCGCCCCGcccccgcgc                                                                                  >  1:10811206/1‑146 (MQ=255)

                                                                    aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATc                                                                                  <  1:2338650/150‑1 (MQ=255)

                                                                     gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTACTCGAGCTGCTGACCGCTGATATCCGCCCGCGGCTGACGCCCGGGCTCACTGCGCGCCCCGccccccgtca                                                                                 >  2:9374198/1‑145 (MQ=255)

                                                                     gCGATTTCTCTGGATTCGACGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGGAATCGAGCTGCGGGCCGCTGAGATCCGCCCGCATCTGGCGCCGGTACTCACTACGCGCCCCGCCCCCCATCa                                                                                 >  1:317971/1‑150 (MQ=255)

                                                                       gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCTATCAAc                                                                               <  2:1755203/150‑1 (MQ=255)

                                                                             tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAaca                                                                         <  1:4182191/150‑1 (MQ=255)

                                                                             tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAaca                                                                         <  1:7467566/150‑1 (MQ=255)

                                                                             tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTACGGTAATCGGGCTGCTGACCGCGGATGTCCGCCCGCATCTGACGCCCGTACTCACTACGCGCCCCGCCCCCAATCGACGTAaca                                                                         >  2:1374347/1‑150 (MQ=255)

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

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

                                                                                       ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTAATCGGGCTGCTGTCCGCTGGTATCCGCCAGCTTCTGACGCCCGTACTCGCTGCGCGCCCCGCCCGCGCGCATCGTAACACTGTGCGGtt                                                               >  1:4013613/1‑150 (MQ=255)

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

                                                                                            aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATa                                                          <  1:9111672/150‑1 (MQ=255)

                                                                                               tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                       <  1:2605247/150‑1 (MQ=255)

                                                                                                gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGc                                                      <  1:1327842/150‑1 (MQ=255)

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

                                                                                                   ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTACACGAGCTGCTGACCGCGGATCTCCGCCAGCCTCTGACGCCCGTCCTCAGTACGCGCCCCGCCCCCCAGCCACGGAACACTTTGCGGGTGAAGACCGCAgg                                                   >  2:9497028/1‑150 (MQ=255)

                                                                                                    cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGc                                                  <  2:7738407/150‑1 (MQ=255)

                                                                                                     aGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGGGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCCGTACTCGCTACGCGCCCCGCCCCCAGTCAACGTAACACTTTGCGGGTTGATACCGCCGGcg                                                 >  2:5918330/1‑149 (MQ=255)

                                                                                                      gCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt                                                <  2:10509612/150‑1 (MQ=255)

                                                                                                       ccTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATg                                               <  2:4650094/150‑1 (MQ=255)

                                                                                                         ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                             <  2:3729830/150‑1 (MQ=255)

                                                                                                         ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                             <  1:8886953/150‑1 (MQ=255)

                                                                                                         ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                             <  2:2975379/150‑1 (MQ=255)

                                                                                                             tgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGTGCTGCTGACCGCTGATATCCGCGAGCTTCTGACGGCAGTAGTCAGTTCGCGCCCCGCGCCCTATCAAGGTTTCTCTTTGCGGTTTAATAGCGCAGGCTTGCTCGAc                                         >  2:6828082/1‑150 (MQ=255)

                                                                                                              gcgcCCAGTGCGGGCCCCACATCCCCCCCCGGACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACa                                        <  2:7191410/150‑1 (MQ=255)

                                                                                                                  ccAGTGCTGGCCCCACATCCCCCCCCGTACGGGCATCGAGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCa                                    >  2:3329520/1‑150 (MQ=255)

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

                                                                                                                   cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGGACTCACTACGCGCCCCGCCCCCCCTCACCGGAACACTTTGCGGTTGAAGACCGCAGGCCTGCACGACAGCCag                                   >  2:8272495/1‑149 (MQ=255)

                                                                                                                   cAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                  >  2:10599501/1‑150 (MQ=255)

                                                                                                                     gTGCTGGCCCCACATCCCCCCCCTTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc                                 <  1:6972608/150‑1 (MQ=255)

                                                                                                                     gTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc                                 <  2:7055903/150‑1 (MQ=255)

                                                                                                                      tGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGcc                                <  2:1952144/150‑1 (MQ=255)

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

                                                                                                                        cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGTCGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTGTAATACCGCAGGCCTGCACGACATCCCTGCCAt                              >  1:10361766/1‑150 (MQ=255)

                                                                                                                        cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAGTAAGCGCCCCGCCCCCAATCAACGTAACACTTTTCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAt                              >  2:7960654/1‑150 (MQ=255)

                                                                                                                        cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAt                              <  2:4126293/150‑1 (MQ=255)

                                                                                                                        cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCAGTCCGCGGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATAAACGGAACACTTTGCGGTTTAATACCGCAGGCCTGCACGACATCCATGCCAt                              >  1:6240087/1‑150 (MQ=255)

                                                                                                                         tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCGGGCCGCGGATATCCGCCCGCATCTGACGCCCGGACTCACTACGCGCCCCCCCCCCAATCAACGGAACACTTTGCGGTTGAATAACGCAGGCATGCACGACAACCCTGCCaac                             >  1:2500296/1‑148 (MQ=255)

                                                                                                                          ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCTACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                            <  1:6725035/150‑1 (MQ=255)

                                                                                                                          ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                            >  2:2174586/1‑150 (MQ=255)

                                                                                                                           gCCCCACATCCCCCCCCGTACCGTAATCGAGCTGCTGACCCCTGATATCCGCCAGCATCTGACGCCAGTAATCACTACGCGCCACGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                           >  2:5660088/1‑150 (MQ=255)

                                                                                                                            ccccACATCCCCCCCCGTACGGTAATCGCGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGc                          <  1:9643118/150‑1 (MQ=255)

                                                                                                                            ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCCGGCGTGCACGACATCCATGCCCTCGGc                          >  1:6609960/1‑150 (MQ=255)

                                                                                                                              ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCCGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                        >  1:556637/1‑150 (MQ=255)

                                                                                                                              ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTGACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                        <  1:7836712/150‑1 (MQ=255)

                                                                                                                                acaTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc                      >  2:1010444/1‑150 (MQ=255)

                                                                                                                                 caTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg                     <  1:104797/150‑1 (MQ=255)

                                                                                                                                    ccccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCAt                  >  1:4604435/1‑150 (MQ=255)

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

                                                                                                                                       cccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAc               <  2:3039939/150‑1 (MQ=255)

                                                                                                                                          ccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCa            >  1:3076393/1‑150 (MQ=255)

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

                                                                                                                                            gTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACg          >  2:8556132/1‑150 (MQ=255)

                                                                                                                                                ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCta      >  2:6127647/1‑149 (MQ=255)

                                                                                                                                                  tAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTcc    >  1:3486010/1‑150 (MQ=255)

AAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAG  >  NC_000913/3727979‑3728276

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