BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,097,359	C→G	5.5%	L30L (CTG→CTC) ‡	rhaB ←	rhamnulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,097,359	0	C	G	5.5%	413.8 / 12.2	128	L30L (CTG→CTC) ‡	rhaB	rhamnulokinase
Reads supporting (aligned to +/- strand): ref base C (54/67); new base G (3/4); total (57/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

Only 100 of 129 total aligned reads displayed.

TCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAC  >  NC_000913/4097210‑4097505

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

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

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

 cAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATAGTTAAAACGATGGATTTCGCGCAGCGTCa                                                                                                                                                   >  2:1098719/1‑150 (MQ=255)

  aaTCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAg                                                                                                                                                  <  1:1336000/150‑1 (MQ=255)

    tCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGc                                                                                                                                                <  1:1761470/150‑1 (MQ=255)

     cccAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGAGTCAGGCt                                                                                                                                               >  2:1367057/1‑150 (MQ=255)

      ccAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGCCATAGCCGTTCTGACTATGCCGCCCATTGTTAAAACGATGGATTTCTCGCAGCGTCAGGCTg                                                                                                                                              >  1:201848/1‑150 (MQ=255)

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

          tGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGc                                                                                                                                          <  2:828162/150‑1 (MQ=255)

          tGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCCTGACGCTGCgcg                                                                                                                                          <  1:1624287/150‑3 (MQ=255)

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

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

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

              aTCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCCTGACGCTGCgcgaaat                                                                                                                                      <  2:848417/150‑7 (MQ=255)

              aTCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATCCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCCTGACGCTGCgcgaaat                                                                                                                                      <  2:471220/150‑7 (MQ=255)

                 gATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACg                                                                                                                                   >  2:1019685/1‑150 (MQ=255)

                    aCGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTc                                                                                                                                >  2:1945085/1‑150 (MQ=255)

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

                      gAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGt                                                                                                                              <  2:649895/150‑1 (MQ=255)

                           ccTTCCTCGCACCCCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGAATTCGCGCAGCGTCAGGCTGCGGCATTAACGCTCGtcccgc                                                                                                                         >  2:72132/1‑145 (MQ=255)

                           ccTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACCTAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAAcgc                                                                                                                         <  2:1606276/150‑1 (MQ=255)

                            cTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAAcgcg                                                                                                                        >  1:1224256/1‑150 (MQ=255)

                                  cGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCa                                                                                                                  <  1:2000134/150‑1 (MQ=255)

                                     acacATTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCAGTTCTGACTATGCAGACGATTGTTAAAACGATGGATTACGCGCAGCGTCAGGCTGCGGCATACACGCTCGTAACGAGAAAGCATCa                                                                                                               >  1:1393421/1‑150 (MQ=255)

                                      cacCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCAc                                                                                                              >  2:536689/1‑150 (MQ=255)

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

                                              tttaTCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGAAATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCAc                                                                                                      <  2:327759/147‑1 (MQ=255)

                                                aaTCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGGGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTg                                                                                                    <  2:1868945/150‑1 (MQ=255)

                                                 acccAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTgg                                                                                                   >  2:1367183/3‑150 (MQ=255)

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

                                                   ccAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGAt                                                                                                 <  1:1505927/150‑1 (MQ=255)

                                                   ccAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGAt                                                                                                 <  1:663047/150‑1 (MQ=255)

                                                    cAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATg                                                                                                >  2:590225/1‑150 (MQ=255)

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

                                                       ggCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATgcgc                                                                                             >  2:673453/1‑150 (MQ=255)

                                                                 aCTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGAAGGATTTCGCGCAGCGTCAGGCTGCGGCCTTCACGCTCCTACCGCGCCAGCATCCCGCGCCCACTGGATGCGCCGGCATCGAc                                                                                   >  1:492471/1‑150 (MQ=255)

                                                                 aCTTTCAAGGCAATCCACATCCCAGGTGCCATAGCCGTTCTGACTATGCAGCCCATTGTAAAAACGATGGAATTCGCTCAGCGTCAGCCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGAc                                                                                   <  2:1941449/150‑1 (MQ=255)

                                                                    ttCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGc                                                                                <  1:956614/150‑1 (MQ=255)

                                                                    ttCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGc                                                                                <  1:1691639/150‑1 (MQ=255)

                                                                     tCAAGGCTATCCACATCCCAGTTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCg                                                                               <  2:1609343/150‑1 (MQ=255)

                                                                     tCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCg                                                                               <  1:602635/150‑1 (MQ=255)

                                                                        aGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGaca                                                                            <  1:1945085/150‑1 (MQ=255)

                                                                        aGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTAAAAACGAGGGATTTCGCGCCGCGTCAGGCTGCGGACTTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGcca                                                                            >  1:289868/1‑150 (MQ=255)

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

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

                                                                              tCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCAGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTg                                                                      <  1:181546/150‑1 (MQ=255)

                                                                               ccACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGc                                                                     <  1:982032/150‑1 (MQ=255)

                                                                               ccACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGc                                                                     >  1:833005/1‑150 (MQ=255)

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

                                                                                cacaTCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCg                                                                    <  2:1097944/150‑1 (MQ=255)

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

                                                                                   aTCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGaaa                                                                 >  1:822080/1‑150 (MQ=255)

                                                                                    tCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGTCGACACAATTGCGAAAg                                                                <  1:1367183/150‑1 (MQ=255)

                                                                                    tCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAg                                                                <  1:1757612/150‑1 (MQ=255)

                                                                                      ccAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGt                                                              >  2:1045809/1‑150 (MQ=255)

                                                                                       cAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCCTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTc                                                             <  2:370627/150‑1 (MQ=255)

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

                                                                                          gTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATa                                                          <  1:695155/150‑1 (MQ=255)

                                                                                          gTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATa                                                          <  1:639606/150‑1 (MQ=255)

                                                                                              cATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCATGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGAGACACAATTGCGAAAGGTCATAAtgt                                                      >  2:77830/1‑150 (MQ=255)

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

                                                                                                  gCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATc                                                  >  2:58427/1‑150 (MQ=255)

                                                                                                  gCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGGCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATc                                                  <  2:1783081/150‑1 (MQ=255)

                                                                                                   ccGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCGGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATcc                                                 <  1:1891279/150‑1 (MQ=255)

                                                                                                   ccGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGTAAGGTCATAATGTGATcc                                                 <  1:1231255/150‑1 (MQ=255)

                                                                                                    cGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCct                                                >  2:1944579/1‑150 (MQ=255)

                                                                                                     gTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCctg                                               <  2:1508181/150‑1 (MQ=255)

                                                                                                          gACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCTAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGaa                                          >  1:1356015/1‑150 (MQ=255)

                                                                                                           aCTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAAt                                         <  2:1948151/150‑1 (MQ=255)

                                                                                                              aTGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTc                                      <  1:1018514/150‑1 (MQ=255)

                                                                                                              aTGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTc                                      >  2:1610257/1‑150 (MQ=255)

                                                                                                               tGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCTCGCTCGTAACGCGCCTGCATGACGCGCCCACTGGATGCGCCGAGATCTACGGCGACACAATTGCGAAAGGTCATAATTTGATCCTGCTGAATTTCa                                     <  2:1902586/150‑1 (MQ=255)

                                                                                                               tGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCa                                     <  1:1573890/150‑1 (MQ=255)

                                                                                                                gCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCAt                                    <  2:1802695/150‑1 (MQ=255)

                                                                                                                    cccATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCTTTATGTGATCCTGCTGAATTTCATTACg                                <  1:1942850/150‑1 (MQ=255)

                                                                                                                      cATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGAc                              <  1:1271584/150‑1 (MQ=255)

                                                                                                                      cATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGAc                              <  2:1015099/150‑1 (MQ=255)

                                                                                                                       aTTGTTAACACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCAAGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGAGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGAcc                             <  2:1095558/150‑1 (MQ=255)

                                                                                                                       aTTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGATTTTCATTACGAcc                             >  1:1173184/1‑150 (MQ=255)

                                                                                                                            tAAAACGATTGATTTCGCGCAGCGTCAGGCTGCGGCTTTCACGCTCGTAACGCGCCAGGATCACGCGCCCACTGGATGCGCCGAGTTCGACTGCGACACAATTGCGAAAGGTCATATTGTGATCCTGCTGAATTTCATTACGACCAGTCt                        <  1:1849807/150‑1 (MQ=255)

                                                                                                                             aaaaCGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTa                       <  2:1472809/150‑1 (MQ=255)

                                                                                                                                  gATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAg                  >  1:27961/1‑150 (MQ=255)

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

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

                                                                                                                                     ggATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGAGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcgc               <  1:1300054/150‑1 (MQ=255)

                                                                                                                                      gATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGcc              >  1:1514928/1‑150 (MQ=255)

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

                                                                                                                                      gATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGcc              >  2:1316712/1‑150 (MQ=255)

                                                                                                                                       aTTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACTGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCt             >  1:1704222/1‑150 (MQ=255)

                                                                                                                                         atgccgcAGCCTGACGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGCCGGCGACACAAGTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGCCCAGTCTAAAAAGCGCCTGa           >  1:619158/5‑150 (MQ=255)

                                                                                                                                         atgccgcAGCCTGACGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGGTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGa           >  1:1179497/5‑150 (MQ=255)

                                                                                                                                           cgcgCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt         >  2:1871440/1‑150 (MQ=255)

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

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

                                                                                                                                              gcAGCGTCAGGCTGCGGCATTCACTCTCGTAACGCGCCTTGATCTCGCGCAGACTTGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTTAATTcg      <  1:143631/150‑1 (MQ=255)

                                                                                                                                              gcAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcg      <  1:1418485/150‑1 (MQ=255)

                                                                                                                                                aGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgcg    <  2:17888/150‑1 (MQ=255)

                                                                                                                                                 gCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGa   <  2:1584106/150‑1 (MQ=255)

                                                                                                                                                  ggTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAc  <  2:415120/149‑1 (MQ=255)

TCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAC  >  NC_000913/4097210‑4097505

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