BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,097,357	G→C	5.5%	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.5%	403.3 / 12.5	128	T31R (ACG→AGG)	rhaB	rhamnulokinase
Reads supporting (aligned to +/- strand): ref base G (54/67); new base C (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 = 4.75e-01

Only 100 of 131 total aligned reads displayed.

TATCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAC  >  NC_000913/4097208‑4097505

tATCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCg                                                                                                                                                      >  2:1300054/1‑150 (MQ=255)

   cAATCCCAATGCTATCGATACGAATGCCTTACTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTAAACACATCCCAGGAGACATAGCCCTTCTGAATATGCAGCCAAATGTTAAAAAGATGCATTTCGAGCAGCgaaa                                                                                                                                                   >  2:939079/1‑147 (MQ=255)

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

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

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

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

    aaTCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAg                                                                                                                                                  <  1:1364851/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:277587/150‑3 (MQ=255)

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

               tATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCAtt                                                                                                                                       <  1:1887795/150‑1 (MQ=255)

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

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

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

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

                      aCGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTc                                                                                                                                <  1:1333001/150‑1 (MQ=255)

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

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

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

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

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

                                  cTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAg                                                                                                                    >  1:1692776/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)

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

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

                                                     ccAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGAt                                                                                                 >  2:863829/1‑150 (MQ=255)

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

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

                                                       aaGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATgc                                                                                               >  1:1802695/1‑150 (MQ=255)

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

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

                                                            gAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCga                                                                                          <  2:1692776/150‑1 (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)

                                                                                tCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCCCGCGCCCACTGGATGCGCCGAGATCGACGGCGCCACAATTg                                                                      >  1:1420531/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                         cAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTc                                                             >  2:286509/1‑150 (MQ=255)

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

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

                                                                                            gTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATa                                                          <  1:639606/150‑1 (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)

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

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

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

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

                                                                                                                        cATTGTTAAAACGTTGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAATGGTCATAATGTGATCCTGCTTAATTTCATTACGAc                              <  1:1542246/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                            tttAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGt                          <  1:1979543/149‑1 (MQ=255)

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

                                                                                                                              tAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCt                        <  2:1601172/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)

                                                                                                                                       ggATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACTGCGACACAATTGCGAAAGGTCATAATTTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcgc               >  2:1124264/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                          aatgccgcAGCCTGACGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATACTGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTg            >  1:968608/6‑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        >  2:943985/1‑150 (MQ=255)

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

                                                                                                                                                gcAGCGGCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGTATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcg      <  1:719280/150‑1 (MQ=255)

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

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

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

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

TATCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAC  >  NC_000913/4097208‑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: `‑`