BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,522,275	T→G	H607P (CAT→CCT)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,522,275	0	T	G	98.2%	756.8 / 1.6	217	H607P (CAT→CCT)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand): ref base T (3/1); new base G (112/100); total (116/101)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 217 total aligned reads displayed.

TAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCT  >  NC_000913/4522179‑4522371

tAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGc                                                                                               >  1:2350430/1‑100 (MQ=255)

tAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGc                                                                                               >  1:188830/1‑100 (MQ=255)

tAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGc                                                                                               >  1:1931981/1‑100 (MQ=255)

   acaCAACCGGTGCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTgg                                                                                            >  2:5092052/3‑100 (MQ=255)

   aTACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTgg                                                                                            <  2:859078/100‑1 (MQ=255)

   aTACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTgg                                                                                            >  2:5054589/1‑100 (MQ=255)

   aTACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTgg                                                                                            >  1:1676826/1‑100 (MQ=255)

      cAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTc                                                                                         <  2:1876201/100‑1 (MQ=255)

       aaCCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCt                                                                                        >  2:2310908/1‑100 (MQ=255)

        aCCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTg                                                                                       <  2:4823068/100‑1 (MQ=255)

         ccGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGTCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTgg                                                                                      <  2:2029810/100‑1 (MQ=255)

           ggTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCa                                                                                    <  1:678723/100‑1 (MQ=255)

           ggTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCa                                                                                    <  2:2538857/100‑1 (MQ=255)

           ggTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCa                                                                                    <  2:3051644/100‑1 (MQ=255)

           ggTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCa                                                                                    >  2:3166537/1‑100 (MQ=255)

           ggTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCa                                                                                    >  2:4136359/1‑100 (MQ=255)

           ggTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCa                                                                                    >  1:1460061/1‑100 (MQ=255)

            gTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCAt                                                                                   <  2:1657025/100‑1 (MQ=255)

              ccctGTCCCACCACTCACGGCCTGAAGCATTGCCCATTGCCGGGTATCATCTGGCAATTCAGGCTCGGGAAGTAAATCGTGAGGGCTGGGTCTGGCAttt                                                                                 <  2:1147025/97‑1 (MQ=255)

              ccATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCAttt                                                                                 <  2:753565/100‑1 (MQ=255)

              ccATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCAttt                                                                                 >  1:936581/1‑100 (MQ=255)

              ccATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCAttt                                                                                 >  2:954516/1‑100 (MQ=255)

              ccATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCAttt                                                                                 <  1:2754404/100‑1 (MQ=255)

                aTGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTaa                                                                               <  2:4900773/100‑1 (MQ=255)

                aTGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTaa                                                                               <  2:2856808/100‑1 (MQ=255)

                    cccACCACTCCCGGAATGAAGCATTGCCCATAGAAGGGTATCCACAGGAAATTCAGGATCGAGCAGTAAATCATGAGGGCTGGGTCTGGCATATAATAtt                                                                           <  2:4931596/100‑1 (MQ=255)

                    cccACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATAtt                                                                           <  2:772756/100‑1 (MQ=255)

                     ccaccaCTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTg                                                                          >  2:704018/1‑100 (MQ=255)

                     ccaccaCTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTg                                                                          >  1:1359886/1‑100 (MQ=255)

                      caccacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGc                                                                         <  1:2062059/100‑1 (MQ=255)

                       accacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCa                                                                        <  2:4323958/100‑1 (MQ=255)

                       accacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCa                                                                        <  2:3514166/100‑1 (MQ=255)

                       accacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCa                                                                        >  1:2298602/1‑100 (MQ=255)

                        ccacTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAg                                                                       <  1:2723475/100‑1 (MQ=255)

                           cgcACGGCCTGAAGCATTTCCCATAGCCGGGTATCACCTGGCAATGCAGGATCGGGGATTAAACCCTGAGGTCTGGGTCTGGCATCGAATGTTGCAGTTg                                                                    >  2:2494855/3‑100 (MQ=255)

                           cTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTg                                                                    >  1:2564273/1‑100 (MQ=255)

                           cTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTg                                                                    >  1:240701/1‑100 (MQ=255)

                              aCGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcct                                                                 >  2:1828089/1‑100 (MQ=255)

                               cGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcctc                                                                >  1:2683877/1‑100 (MQ=255)

                                 gCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcctcct                                                              >  2:4439585/1‑100 (MQ=255)

                                   cTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcctcctcc                                                            <  1:712787/100‑1 (MQ=255)

                                     gAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCtg                                                          <  1:4074990/100‑1 (MQ=255)

                                          aTTGTCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAAGATTGCAGTTGCCTCCTCCTGTGAAg                                                     <  1:4991114/100‑1 (MQ=255)

                                          aTTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAg                                                     >  1:2171023/1‑100 (MQ=255)

                                          aTTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAg                                                     <  2:2434453/100‑1 (MQ=255)

                                            tGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGt                                                   >  2:1486180/1‑100 (MQ=255)

                                             gCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTa                                                  >  1:4927450/1‑100 (MQ=255)

                                             gCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTa                                                  >  1:2044014/1‑100 (MQ=255)

                                             gCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTa                                                  >  1:2934943/1‑100 (MQ=255)

                                             gCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTa                                                  >  1:3880522/1‑100 (MQ=255)

                                             gCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTa                                                  >  1:1169063/1‑100 (MQ=255)

                                             gCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGCAAATCCTGAGGGCTGGGTCTGGCATTTACTAATGCAGTTGCCTCCTCCTGTGAAGGTa                                                  >  2:5229003/1‑100 (MQ=255)

                                              cccATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTaa                                                 >  1:1094185/1‑100 (MQ=255)

                                              cccATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTCAAGGTaa                                                 >  1:3743088/1‑100 (MQ=255)

                                               ccATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAAt                                                <  1:3520235/100‑1 (MQ=255)

                                               ccATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAAt                                                <  1:5050214/100‑1 (MQ=255)

                                                cATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGTTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTCAtc                                               >  2:1737045/1‑99 (MQ=255)

                                                 aTAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTg                                              <  2:3348016/100‑1 (MQ=255)

                                                    gCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTc                                           >  2:849236/1‑100 (MQ=255)

                                                    gCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTc                                           <  1:4476322/100‑1 (MQ=255)

                                                     ccGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCa                                          <  1:3053595/100‑1 (MQ=255)

                                                      cGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTCATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCAt                                         <  2:4221753/100‑1 (MQ=255)

                                                            tcatcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTAc                                   >  2:3248475/1‑100 (MQ=255)

                                                             catcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACg                                  >  2:2174102/1‑100 (MQ=255)

                                                             catcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACg                                  >  1:1626842/1‑100 (MQ=255)

                                                                 tGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGtt                              >  2:4851908/1‑100 (MQ=255)

                                                                 tGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGtt                              <  1:3771927/100‑1 (MQ=255)

                                                                   gCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTcc                            <  1:1215564/100‑1 (MQ=255)

                                                                    cAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCg                           <  1:4700376/100‑1 (MQ=255)

                                                                     aaTTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGa                          <  1:3196093/100‑1 (MQ=255)

                                                                      aTTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAg                         >  2:2788948/1‑100 (MQ=255)

                                                                      aTTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAg                         >  2:4473246/1‑100 (MQ=255)

                                                                        tCAGGATCGGGAAGGAAATCCGGAGGGCTGGGTCTGGCATTTAATATTGTAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTTCGAGGa                       <  2:4394094/100‑1 (MQ=255)

                                                                         cAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGaa                      >  2:4049896/1‑100 (MQ=255)

                                                                            gATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAAtt                   <  1:2992949/100‑1 (MQ=255)

                                                                             aTCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTg                  <  2:391038/100‑1 (MQ=255)

                                                                             aTCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTg                  <  1:4207853/100‑1 (MQ=255)

                                                                             aTCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTg                  <  1:3347573/100‑1 (MQ=255)

                                                                              tCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTGAATATGGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCGGAGGAAATTGa                 >  2:2354422/1‑100 (MQ=255)

                                                                                gggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACt               >  2:2556840/1‑100 (MQ=255)

                                                                                gggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACt               >  1:3217026/1‑100 (MQ=255)

                                                                                gggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACt               <  2:4973511/100‑1 (MQ=255)

                                                                                 ggAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACtt              >  1:3709972/1‑100 (MQ=255)

                                                                                 ggAAGTAAATCCTGACGGCCGGGTCCGGCATTTAATATTCCAGCTGCCTCCCCCTGTGAAGGTAATTATTCATCGCTACGGTTTCCGCGGCAAttccctt              >  2:399827/1‑95 (MQ=255)

                                                                                     gTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGcc          <  2:547138/100‑1 (MQ=255)

                                                                                     gTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGcc          >  2:2963093/1‑100 (MQ=255)

                                                                                     gTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGcc          >  1:5108333/1‑100 (MQ=255)

                                                                                     gTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGcc          >  1:5402470/1‑100 (MQ=255)

                                                                                     gTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGcc          <  1:1424469/100‑1 (MQ=255)

                                                                                     gTAAAGCCGGAGGGCGGGGTCTGGCATGTAAGATCGCAGTGCCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGcc          <  2:5243102/100‑1 (MQ=255)

                                                                                      tAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTCAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCg         <  2:3708709/100‑1 (MQ=255)

                                                                                       aaaTCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGt        >  1:1714480/1‑100 (MQ=255)

                                                                                        aaTCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTg       <  1:3338160/100‑1 (MQ=255)

                                                                                         aTCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGa      >  2:4699523/1‑100 (MQ=255)

                                                                                         aTCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGa      >  1:5549480/1‑100 (MQ=255)

                                                                                          tCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGaa     <  2:3300052/100‑1 (MQ=255)

                                                                                            cTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGCAGc   >  2:4604328/1‑100 (MQ=255)

                                                                                            cTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGc   >  1:3246584/1‑100 (MQ=255)

                                                                                            cTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGc   >  2:1012966/1‑100 (MQ=255)

                                                                                             tGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCt  >  2:2076247/1‑100 (MQ=255)

TAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCT  >  NC_000913/4522179‑4522371

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