BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,097,357	G→C	5.4%	T31R (ACG→AGG)	rhaB ←	rhamnulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,097,357	0	G	C	5.4%	507.4 / 14.0	148	T31R (ACG→AGG)	rhaB	rhamnulokinase
Reads supporting (aligned to +/- strand): ref base G (80/60); new base C (3/5); total (83/65)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

Only 100 of 153 total aligned reads displayed.

TATCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACC  >  NC_000913/4097208‑4097506

  tCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTc                                                                                                                                                     <  1:1859174/150‑1 (MQ=255)

   cAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCa                                                                                                                                                    >  2:2114021/1‑150 (MQ=255)

       cccAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCt                                                                                                                                                >  2:162448/1‑150 (MQ=255)

          aaTGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCg                                                                                                                                             <  2:646785/150‑1 (MQ=255)

           aTGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCgg                                                                                                                                            <  2:2329389/150‑1 (MQ=255)

            tGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCCTGACGCTGCgcg                                                                                                                                           <  2:67234/150‑3 (MQ=255)

            tGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAACCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCCTGACGCTGCgcg                                                                                                                                           <  1:1819101/150‑3 (MQ=255)

             gCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCa                                                                                                                                          <  2:2371551/150‑1 (MQ=255)

               tATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCAtt                                                                                                                                        >  1:2416025/1‑150 (MQ=255)

                aTCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTc                                                                                                                                       <  1:1206341/150‑1 (MQ=255)

                 tCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCa                                                                                                                                      >  1:1738126/1‑150 (MQ=255)

                  cGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCAc                                                                                                                                     >  2:2231987/1‑150 (MQ=255)

                   caTACGAATACCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCCCg                                                                                                                                    >  2:175175/2‑150 (MQ=255)

                       cGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCg                                                                                                                                <  1:905678/150‑1 (MQ=255)

                        gAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGt                                                                                                                               >  2:1086689/1‑150 (MQ=255)

                        gAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGt                                                                                                                               >  2:2323581/1‑150 (MQ=255)

                          aTCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTaa                                                                                                                             >  2:852239/1‑150 (MQ=255)

                            cccTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAAcg                                                                                                                           >  2:2055026/1‑150 (MQ=255)

                            cccTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAAcg                                                                                                                           >  1:800626/1‑150 (MQ=255)

                             ccTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGTCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAAcgc                                                                                                                          >  1:1877216/1‑150 (MQ=255)

                              cTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAAcgcg                                                                                                                         <  2:1331645/150‑1 (MQ=255)

                                  cTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAg                                                                                                                     <  1:2128010/150‑1 (MQ=255)

                                  cTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAg                                                                                                                     >  1:2189901/1‑150 (MQ=255)

                                      cacaCCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATc                                                                                                                 >  2:44005/1‑150 (MQ=255)

                                         acCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCAcg                                                                                                              >  2:1023880/1‑150 (MQ=255)

                                         acCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCAcc                                                                                                              >  2:506838/1‑149 (MQ=255)

                                            ttGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGTCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACgcgc                                                                                                           <  2:1877216/150‑1 (MQ=255)

                                              gTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGccc                                                                                                         >  1:1130907/1‑150 (MQ=255)

                                                 tAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGATAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACt                                                                                                      >  1:1577374/1‑150 (MQ=255)

                                                  aaTCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTg                                                                                                     <  1:971615/150‑1 (MQ=255)

                                                   aTCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCTACTgg                                                                                                    <  2:1104899/150‑1 (MQ=255)

                                                   aTCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTgg                                                                                                    >  2:15668/1‑150 (MQ=255)

                                                    tCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGa                                                                                                   >  1:1650410/1‑150 (MQ=255)

                                                      cAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATg                                                                                                 <  2:2376609/150‑1 (MQ=255)

                                                        aGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATgcg                                                                                               <  2:1747181/150‑1 (MQ=255)

                                                           cGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCg                                                                                            >  1:2317853/1‑150 (MQ=255)

                                                                 gCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCg                                                                                      >  1:1093967/1‑150 (MQ=255)

                                                                 gCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCg                                                                                      >  1:1302630/1‑150 (MQ=255)

                                                                  cACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGa                                                                                     >  2:824729/1‑150 (MQ=255)

                                                                        cAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGa                                                                               <  2:1340596/150‑1 (MQ=255)

                                                                        cAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGa                                                                               >  1:577016/1‑150 (MQ=255)

                                                                         aaGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAGCGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGac                                                                              >  2:931033/1‑150 (MQ=255)

                                                                           ggCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGacac                                                                            >  2:2337649/1‑150 (MQ=255)

                                                                             cTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACaa                                                                          <  1:1067028/150‑1 (MQ=255)

                                                                              tATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAAt                                                                         >  2:1432478/1‑150 (MQ=255)

                                                                               aTCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAAtt                                                                        >  1:1316252/1‑150 (MQ=255)

                                                                                  cacaTCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCg                                                                     >  1:1711002/1‑150 (MQ=255)

                                                                                   acaTCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGa                                                                    <  1:1018843/150‑1 (MQ=255)

                                                                                       cccAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACTCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAgg                                                                <  2:1192923/150‑1 (MQ=255)

                                                                                           ggTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCAt                                                            <  1:1312556/150‑1 (MQ=255)

                                                                                           ggTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCAt                                                            >  1:1455791/1‑150 (MQ=255)

                                                                                               aCATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAAtg                                                        <  2:1037673/150‑1 (MQ=255)

                                                                                               aCATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAAtg                                                        >  2:2457252/1‑150 (MQ=255)

                                                                                                cATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAAtgt                                                       <  1:967550/150‑1 (MQ=255)

                                                                                                cATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAAtgt                                                       >  2:2190351/1‑150 (MQ=255)

                                                                                                cATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAAtgt                                                       <  2:13097/150‑1 (MQ=255)

                                                                                                  tAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGa                                                     <  1:488369/150‑1 (MQ=255)

                                                                                                  tAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGCGa                                                     <  1:281267/150‑1 (MQ=255)

                                                                                                   aGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAAtgtgct                                                    >  1:683278/1‑148 (MQ=255)

                                                                                                   aGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGAt                                                    <  2:1262235/150‑1 (MQ=255)

                                                                                                     ccGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGAGACACAATTGCGAAAGGTCATAATGTGATcc                                                  <  2:1619590/150‑1 (MQ=255)

                                                                                                       gTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCctg                                                >  2:1660469/1‑150 (MQ=255)

                                                                                                        ttCTGACTATTCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGAATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATTTGATCctgc                                               <  1:388174/150‑1 (MQ=255)

                                                                                                          cTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCctgctg                                             >  1:1025305/1‑150 (MQ=255)

                                                                                                           tGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGa                                            >  1:1341279/1‑150 (MQ=255)

                                                                                                            gACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGaa                                           >  2:1236382/1‑150 (MQ=255)

                                                                                                            gACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGaa                                           >  1:1945772/1‑150 (MQ=255)

                                                                                                              cTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAAtt                                         >  1:104203/1‑150 (MQ=255)

                                                                                                              cTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAAtt                                         <  1:513169/150‑1 (MQ=255)

                                                                                                              cTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAAtt                                         <  1:1279811/150‑1 (MQ=255)

                                                                                                                 tGCAGCCCATTGTTAAAAAGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACTGCGACACAATTGCTAAAGGTCATAATGTGATCCTGCTGAATTTCa                                      <  2:545457/150‑1 (MQ=255)

                                                                                                                   cAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCAtt                                    >  1:37559/1‑150 (MQ=255)

                                                                                                                   cAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCAtt                                    <  1:162448/150‑1 (MQ=255)

                                                                                                                     gCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTAc                                  >  1:601384/1‑150 (MQ=255)

                                                                                                                       ccATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGa                                >  1:959857/1‑150 (MQ=255)

                                                                                                                        cATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGAc                               >  2:846637/1‑150 (MQ=255)

                                                                                                                            gTTAAAACTATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCAGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGt                           >  2:2090144/1‑150 (MQ=255)

                                                                                                                             ttAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTc                          <  1:1660469/150‑1 (MQ=255)

                                                                                                                             ttAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTc                          <  1:1665630/150‑1 (MQ=255)

                                                                                                                             ttAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTc                          <  2:2416025/150‑1 (MQ=255)

                                                                                                                              tAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCt                         <  1:1727404/150‑1 (MQ=255)

                                                                                                                              tAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCt                         >  1:1187469/1‑150 (MQ=255)

                                                                                                                              tAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCt                         >  2:356973/1‑150 (MQ=255)

                                                                                                                                aaaCGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTaa                       >  1:716899/1‑150 (MQ=255)

                                                                                                                                   cGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaaa                    >  2:314471/1‑150 (MQ=255)

                                                                                                                                      tGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                 <  2:722049/150‑1 (MQ=255)

                                                                                                                                      tGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                 >  2:712043/1‑150 (MQ=255)

                                                                                                                                        gATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGcc               <  2:1683847/150‑1 (MQ=255)

                                                                                                                                         aTTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACAGAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCt              >  2:806253/1‑150 (MQ=255)

                                                                                                                                         aTTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCt              >  1:83662/1‑150 (MQ=255)

                                                                                                                                         aTTTAGCGCAGCGTCAGCATGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCTCCCACTGGATGCGCCGAGATAGAATGCGACACAAATTCGAAAGGTCATAATGTAATACTGCTCAAATACATTACGACCATTATAAAAAGCGCCt              >  2:613083/1‑150 (MQ=255)

                                                                                                                                          tttCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTg             <  1:1708216/150‑1 (MQ=255)

                                                                                                                                            tCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGaa           <  2:1711002/150‑1 (MQ=255)

                                                                                                                                             gccgcAGCCTGACGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt          >  1:2314698/3‑150 (MQ=255)

                                                                                                                                             cgcgCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt          <  2:1957070/150‑1 (MQ=255)

                                                                                                                                              gcgcAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt         >  1:1231647/1‑150 (MQ=255)

                                                                                                                                              ccgcAGCCTGACGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt         >  2:346177/2‑150 (MQ=255)

                                                                                                                                                 cAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgc      <  1:1240475/150‑1 (MQ=255)

                                                                                                                                                 cAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgc      <  1:330783/150‑1 (MQ=255)

                                                                                                                                                     gTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAcc  <  1:314471/150‑1 (MQ=255)

TATCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACC  >  NC_000913/4097208‑4097506

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