BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,109	0	A	C	8.2%	539.4 / 4.8	199	D421E (GAT→GAG)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base A (67/109); new base C (17/3); total (87/112)
Fisher's exact test for biased strand distribution p-value = 6.73e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

GTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGC  >  NC_000913/3727960‑3728258

gTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCcaca                                                                                                                                                        >  2:2801186/1‑150 (MQ=255)

 tttCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCcacag                                                                                                                                                       >  1:2967045/1‑149 (MQ=255)

 tttCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCcacag                                                                                                                                                       >  1:3719222/1‑149 (MQ=255)

 tttCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCcacag                                                                                                                                                       >  1:1943567/1‑149 (MQ=255)

  ttCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATc                                                                                                                                                      >  2:16286/1‑150 (MQ=255)

      cgtcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTccccc                                                                                                                                                  >  1:3272186/1‑150 (MQ=255)

      cgtcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccc                                                                                                                                                  <  2:3551406/150‑1 (MQ=255)

      cgtcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccc                                                                                                                                                  <  1:702179/150‑1 (MQ=255)

       gtcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCCGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccc                                                                                                                                                 <  2:2494521/150‑1 (MQ=255)

        tcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccc                                                                                                                                                <  1:3313922/150‑1 (MQ=255)

         cgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccc                                                                                                                                               >  2:1234850/1‑150 (MQ=255)

         cgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccc                                                                                                                                               <  1:3227863/150‑1 (MQ=255)

          gtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCGCGCCCg                                                                                                                                              >  1:2819103/1‑150 (MQ=255)

           tGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGt                                                                                                                                             <  2:211935/150‑1 (MQ=255)

           tGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGt                                                                                                                                             >  2:4086743/1‑150 (MQ=255)

               tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCggccgg                                                                                                                                         >  2:3833366/1‑145 (MQ=255)

               tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCggccgg                                                                                                                                         >  2:226793/1‑145 (MQ=255)

               tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACgg                                                                                                                                         >  2:2181353/1‑150 (MQ=255)

                gATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCggggggg                                                                                                                                        >  1:1371592/1‑144 (MQ=255)

                   aaGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAt                                                                                                                                     <  1:1431512/150‑1 (MQ=255)

                    aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGTGCCCGCCCGTGCggggggg                                                                                                                                    >  2:960737/1‑145 (MQ=255)

                    aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACggggatc                                                                                                                                    >  2:2701616/1‑145 (MQ=255)

                     gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCATCCCCCCCCGTACGGtgcgcg                                                                                                                                   >  2:3771168/1‑145 (MQ=255)

                     gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTGCGGtgctcg                                                                                                                                   >  2:2279945/1‑145 (MQ=255)

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

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

                            aaCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgctgct                                                                                                                            <  1:1311166/150‑1 (MQ=255)

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

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

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

                                      aTCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGat                                                                                                                  <  2:2882345/150‑1 (MQ=255)

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

                                                   gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                     >  2:3275163/1‑134 (MQ=255)

                                                   gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGAGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCGGACCGCTGATATCCGCCAGCATc                                                                                                     <  2:3372074/150‑1 (MQ=255)

                                                    aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGTCCGCTGGGATGCGCCGGcggcg                                                                                                    >  2:698480/1‑146 (MQ=255)

                                                     cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTg                                                                                                   <  2:422084/150‑1 (MQ=255)

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

                                                          cTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTTCGGTATTCGAGCTGCTGTCCGCTGCTCTCCGCCGGCCTCTGTCggc                                                                                              >  2:63351/1‑148 (MQ=255)

                                                           tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTGCGGTGATCGGGCTGCTGACCGCTGGTAGCCGCCGGCCTCTGACggcc                                                                                             >  2:2173680/1‑147 (MQ=255)

                                                           aaaCGGTACTAGCCGCAACAATTCAAACAGCGATTTCTCTGGATTCGCCGCGAACTGCGCCACCCTAGCTGCGCCCAGTCCTGGCCCCACATCCCCCCCCGTACGCTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                             <  2:3926031/149‑1 (MQ=255)

                                                             aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGGCGGGCATCGCGCTGCTGTCGGCTGAGCTCCGCCAGCATCTGGCGCCGGt                                                                                           >  2:36030/1‑150 (MQ=255)

                                                             aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGGTGGTGACCGCTGATATCCGCCAGCTTCTGACGCCAgg                                                                                           >  2:3977164/1‑149 (MQ=255)

                                                                  aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                                      <  1:987560/150‑1 (MQ=255)

                                                                   gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGGACGGTTGTCGGGGTGCTGTCCGCTGGGAGCCGCCGGCGGCTGGCGCGAGTACTctc                                                                                     >  1:2171045/1‑148 (MQ=255)

                                                                   gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGAGTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                     <  2:1946936/150‑1 (MQ=255)

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

                                                                          caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAGCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACgcgc                                                                              <  1:2971520/150‑1 (MQ=255)

                                                                            acaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc                                                                            <  1:2950148/150‑1 (MQ=255)

                                                                             caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgcccc                                                                           <  1:2872296/150‑1 (MQ=255)

                                                                             caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgcccc                                                                           <  1:231580/150‑1 (MQ=255)

                                                                               aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTAGTCGAGCTGGTGCCCGCTGAGGGCCGCCAGCAGCTGACGCCGGTTCTCACTCCGcgccccgc                                                                         >  2:1319563/1‑150 (MQ=255)

                                                                                ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgcc                                                                        <  1:1379717/150‑1 (MQ=255)

                                                                                     tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCac                                                                   >  2:2924203/1‑149 (MQ=255)

                                                                                     tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCaa                                                                   <  1:3381491/150‑1 (MQ=255)

                                                                                        gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCa                                                                <  1:3536087/150‑1 (MQ=255)

                                                                                        gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCa                                                                <  2:1033957/150‑1 (MQ=255)

                                                                                         cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGGAATCGGGCTGCGGGCCGCTGATCTCCGCCCGCCTCTGACGCCAGTCCTCACTCCGCGCCCCGcccccccgcaa                                                               >  2:1500717/1‑144 (MQ=255)

                                                                                          gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCGAGTGCTGGCCCCACATCCCCCCCCGTACGGTTATCGAGCTGCTGACCGCTGATATCCGCCTGCATCTGTCGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                              <  2:3240067/150‑1 (MQ=255)

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

                                                                                                  tGGATTCGCCGCGATCATCGCCAGCCTTGCTGCTCCCAGTGCTGGCCCCACATCCCCACCCGTACGGAAATCGAGCAGGTGACCGCGGATATCAGCCAGCAGCGGGCGCCAGTACTCAAAACGCGCCCCGCCCCCATGCAACGTCacacg                                                      >  2:1376471/1‑149 (MQ=255)

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

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

                                                                                                          ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGAGCTGCTGTCCGCTGATATCCGGCAGCGTCTGTCGGGAGTACGGGCTGCGCGCCCCGCCCCCAATCGACGGACCACTTTGCgggt                                              >  2:300592/1‑148 (MQ=255)

                                                                                                          ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCTCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGt‑t                                             <  2:3100352/150‑1 (MQ=255)

                                                                                                          ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGt‑t                                             <  1:2037169/150‑1 (MQ=255)

                                                                                                            gcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGTTATCCGCCAGCGTCGGTCGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTGTGCGGg‑ttt                                           >  2:1724246/1‑149 (MQ=255)

                                                                                                             cgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTaa                                          <  1:2689528/150‑1 (MQ=255)

                                                                                                               aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCTGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATa                                        <  2:3310117/150‑1 (MQ=255)

                                                                                                                tCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATAc                                       <  2:50219/150‑1 (MQ=255)

                                                                                                                tCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATAc                                       <  1:1568423/150‑1 (MQ=255)

                                                                                                                     gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAg                                  <  1:739427/150‑1 (MQ=255)

                                                                                                                      ccAGCGTGGCTGCGCCCAGGGCTGGCCCCACATCCCCCCCCGGACGGGACTCGAGCTGCGGACCGCTGATATCCGCCAGCATCTGACTCCAGTACTCACTACGCGCCCCGCCCCCAATCCACGTAACACTTTGCGGTATTAAT‑CCGCaaa                                 <  2:4005872/150‑3 (MQ=255)

                                                                                                                      ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑ATTAAACCGCaaa                                 <  2:597919/150‑3 (MQ=255)

                                                                                                                         gCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCAt                              <  2:2587408/150‑1 (MQ=255)

                                                                                                                            ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCa                           <  1:1251546/150‑1 (MQ=255)

                                                                                                                               ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGa                        <  2:345176/150‑1 (MQ=255)

                                                                                                                               ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGa                        <  2:720406/150‑1 (MQ=255)

                                                                                                                               ctgcGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTCCGGTAAGCGAGCTGCTGCCCGCTGATATCCGCCCGCATCTGGCGCCAGGACTCCCGACGCGCCCCGCCCCCCAGCAACGTAACACGGTGCGGG‑TGACTACCGCAGGCCTGCCCgcc                       >  1:635003/1‑148 (MQ=255)

                                                                                                                                   gcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACGTTGCGGT‑TTAATACCGCAGGCATGCACGACATc                    >  2:1039183/1‑150 (MQ=255)

                                                                                                                                    cccAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGAGCTGCTGACCGCGGATATCCGCCAGCATCTGACGCCAGTACTCACGACGCGCCCCGCCCCCAAGCAACGGAACACTTTGCGGT‑TTAATACCGCAGGCCTGCACGACATcc                   >  1:1105659/1‑150 (MQ=255)

                                                                                                                                     ccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCa                  >  2:3600288/1‑150 (MQ=255)

                                                                                                                                     ccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCa                  >  2:763757/1‑150 (MQ=255)

                                                                                                                                      cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTCCGCACTACGCGCCCCGCCCCCACTCAACGTAACACGTTGCGGT‑GTAATACCGCAGGCATGCCCGACATCCAt                 >  2:867680/1‑150 (MQ=255)

                                                                                                                                      cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCAt                 <  1:4049426/150‑1 (MQ=255)

                                                                                                                                      cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCAt                 <  1:3452463/150‑1 (MQ=255)

                                                                                                                                       aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATg                <  2:208626/150‑1 (MQ=255)

                                                                                                                                        gTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGc               <  2:514895/150‑1 (MQ=255)

                                                                                                                                           cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACGTTGCGGT‑TTAATACCGCAGGCATGAACGACATCCAAGccct            >  2:2355513/1‑148 (MQ=255)

                                                                                                                                            tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTACCACTGTGCGGT‑TTAATACCGCCGGCATGCCCGACATCCATGCCagc           >  2:3064879/1‑148 (MQ=255)

                                                                                                                                             ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATCt          <  1:3076354/150‑1 (MQ=255)

                                                                                                                                              gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATcag         >  1:2469808/1‑148 (MQ=255)

                                                                                                                                              gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGG‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTg         >  2:261143/1‑150 (MQ=255)

                                                                                                                                              gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGG‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTg         >  2:2676833/1‑150 (MQ=255)

                                                                                                                                               ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTGc        >  2:221832/1‑150 (MQ=255)

                                                                                                                                               ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGGAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTGc        >  1:1295936/1‑150 (MQ=255)

                                                                                                                                                cccACATCCCCCCCCCTACGGTAAACGAGCTGCTGACCGCTGATAACAGCCAGCATCTGACGCCAGTCATCACTACGCGCCCCGCCCCCAATCAACCAAACACTATGCGGT‑TTAAAACCGCAGGCATGCACGACAACCATGCCAACAGcc       >  1:122407/1‑150 (MQ=255)

                                                                                                                                                 ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATAACCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa      >  2:3289369/1‑150 (MQ=255)

                                                                                                                                                   acaTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc    >  1:3998115/1‑150 (MQ=255)

                                                                                                                                                    caTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg   <  1:2745361/150‑1 (MQ=255)

                                                                                                                                                    caTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCACCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg   <  1:3279171/150‑1 (MQ=255)

GTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGT‑TTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGC  >  NC_000913/3727960‑3728258

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