BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,097,361	G→C	5.7%	L30V (CTG→GTG) ‡	rhaB ←	rhamnulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,097,361	0	G	C	5.7%	383.5 / 13.0	122	L30V (CTG→GTG) ‡	rhaB	rhamnulokinase
Reads supporting (aligned to +/- strand): ref base G (49/65); new base C (3/4); total (53/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 125 total aligned reads displayed.

AATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCT  >  NC_000913/4097212‑4097510

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

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

    ccAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGCCATAGCCGTTCTGACTATGCCGCCCATTGTTAAAACGATGGATTTCTCGCAGCGTCAGGCTg                                                                                                                                                   >  1:201848/1‑150 (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                                                                                                                                           <  1:1248411/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                     ggCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATgcgc                                                                                                  >  2:673453/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:1691639/150‑1 (MQ=255)

                                                                   tCAAGGCTATCCACATCCCAGTTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCg                                                                                    <  2:1609343/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)

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

                                                                            tCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTg                                                                           <  2:1080014/150‑1 (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:135228/1‑150 (MQ=255)

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

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

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

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

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

                                                                                       ggTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCAt                                                                <  2:765264/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)

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

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

                                                                                                            aTGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTc                                           <  1:1018514/150‑1 (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)

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

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

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

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

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

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

                                                                                                                           aaaaCGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTa                            <  2:1472809/150‑1 (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                   >  1:1514928/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             >  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)

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

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

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

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

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

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

                                                                                                                                                     ggCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCAGAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGAGCCCTgct  >  2:301313/1‑150 (MQ=255)

AATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCT  >  NC_000913/4097212‑4097510

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