BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	238,990	C→T	pseudogene (113/357 nt)	yafU ←	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	238,990	0	C	T	98.0%	519.4 / ‑1.3	147	pseudogene (113/357 nt)	yafU	pseudogene
Reads supporting (aligned to +/- strand): ref base C (2/1); new base T (67/77); total (69/78)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

TAAAACAACTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTCCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAA  >  NC_000913/238893‑239086

 aaaaCAACTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCaa                                                                                               >  2:1201419/1‑100 (MQ=255)

 aaaaCAACTGTACACACCATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCaa                                                                                               >  2:866005/1‑100 (MQ=255)

  aaaCAACTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAAt                                                                                              <  2:79972/100‑1 (MQ=255)

    acaacTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAg                                                                                            >  1:2157369/1‑100 (MQ=255)

    acaacTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAg                                                                                            >  1:1919941/1‑100 (MQ=255)

     caacTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGc                                                                                           <  1:558747/100‑1 (MQ=255)

     caacTGTACACACCATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGc                                                                                           <  1:866005/100‑1 (MQ=255)

      aacTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGc                                                                                           <  1:1004257/99‑1 (MQ=255)

       acTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCgg                                                                                         >  2:3364478/1‑100 (MQ=255)

       acTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCgg                                                                                         <  1:1268032/100‑1 (MQ=255)

       acTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCgg                                                                                         >  2:2497421/1‑100 (MQ=255)

           tACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGc                                                                                     <  2:489640/100‑1 (MQ=255)

           tACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGc                                                                                     <  2:3005868/100‑1 (MQ=255)

            acacacTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACGCAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCa                                                                                    <  1:2405917/100‑1 (MQ=255)

             cacacTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCaa                                                                                   >  2:1249753/1‑100 (MQ=255)

             cacacTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCaa                                                                                   <  2:1469122/100‑1 (MQ=255)

              acacTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAac                                                                                  <  2:1403209/100‑1 (MQ=255)

              acacTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAac                                                                                  <  1:92495/100‑1 (MQ=255)

              acacTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAac                                                                                  <  2:1647841/100‑1 (MQ=255)

              acacTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAac                                                                                  <  1:3032652/100‑1 (MQ=255)

               cacTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTGCTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAaca                                                                                 >  2:2825918/1‑100 (MQ=255)

               cacTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAaca                                                                                 >  2:90052/1‑100 (MQ=255)

                acTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAacac                                                                                >  1:872196/1‑100 (MQ=255)

                acTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAacac                                                                                >  1:3288679/1‑100 (MQ=255)

                acTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAacac                                                                                >  2:2485806/1‑100 (MQ=255)

                       cTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGc                                                                         >  1:123686/1‑100 (MQ=255)

                       cTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGc                                                                         <  1:3364491/100‑1 (MQ=255)

                           cAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCtata                                                                     >  1:1967177/1‑100 (MQ=255)

                           cAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCtata                                                                     <  1:278496/100‑1 (MQ=255)

                           cAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCtata                                                                     <  2:222098/100‑1 (MQ=255)

                             gTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCCGACTTTCAATAGCGGTTGCAACACCACCGGCTATAAc                                                                   >  2:2286567/1‑100 (MQ=255)

                             gTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGCCTTTCAATAGCGGTTGCAACACCACCGGCTATAAc                                                                   >  1:1542513/1‑100 (MQ=255)

                             gTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAAc                                                                   <  1:1222447/100‑1 (MQ=255)

                             gTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAAc                                                                   >  2:2929098/1‑100 (MQ=255)

                               aaaaTTGCTCCTCCGACAAGAATTACTGGACAGTCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCTGCTATAACTg                                                                 >  2:2242651/1‑100 (MQ=255)

                               aaaaTTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCACCACCACCGGCTATAACTg                                                                 <  2:3143547/100‑1 (MQ=255)

                                 aaTTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCt                                                               <  1:2056214/100‑1 (MQ=255)

                                  aTTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCtt                                                              >  1:2606714/1‑100 (MQ=255)

                                      ctcctcCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGc                                                          <  1:3308106/100‑1 (MQ=255)

                                      ctcctcCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGc                                                          <  1:2937924/100‑1 (MQ=255)

                                      ctcctcCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGc                                                          <  2:1444153/100‑1 (MQ=255)

                                      ctcctcCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTCCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGc                                                          >  2:1653746/1‑100 (MQ=255)

                                       tcctccGACAAGAATTACTGTAACGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGcc                                                         <  1:1749837/100‑1 (MQ=255)

                                       tcctccGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGcc                                                         <  1:3290778/100‑1 (MQ=255)

                                       tcctccGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGcc                                                         <  1:777216/100‑1 (MQ=255)

                                        cctccGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCca                                                        >  1:1000553/1‑100 (MQ=255)

                                        cctccGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCca                                                        >  1:82269/1‑100 (MQ=255)

                                         ctccGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACTCCACCGGCTATAACTGCTTTGGCcac                                                       <  1:2071337/100‑1 (MQ=255)

                                         ctccGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCcaa                                                       >  1:2214642/1‑99 (MQ=255)

                                         ctccGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCcaa                                                       >  2:837283/1‑99 (MQ=255)

                                            cGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTCCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATc                                                    <  2:2346416/100‑1 (MQ=255)

                                                aaGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCAtt                                                <  1:1165121/100‑1 (MQ=255)

                                                aaGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCAtt                                                <  2:2559927/100‑1 (MQ=255)

                                                 aGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTg                                               >  1:2757246/1‑100 (MQ=255)

                                                    aTTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGaaa                                            <  1:948072/100‑1 (MQ=255)

                                                     ttACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGaaaa                                           <  2:805622/100‑1 (MQ=255)

                                                        cTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATc                                        <  1:2807830/100‑1 (MQ=255)

                                                         tGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCg                                       <  2:1505277/100‑1 (MQ=255)

                                                          ggAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGc                                      >  2:1810811/1‑100 (MQ=255)

                                                          ggAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGc                                      >  2:1357681/1‑100 (MQ=255)

                                                          ggAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGc                                      >  1:2514983/1‑100 (MQ=255)

                                                           gAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACTGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGcc                                     <  2:3288003/100‑1 (MQ=255)

                                                           gAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGcc                                     <  2:2027754/100‑1 (MQ=255)

                                                            aaaGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCa                                    <  1:591251/100‑1 (MQ=255)

                                                              aGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAg                                  >  2:3016794/1‑100 (MQ=255)

                                                              aGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAg                                  >  1:217446/1‑100 (MQ=255)

                                                              aGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAg                                  <  1:3280426/100‑1 (MQ=255)

                                                               gCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGa                                 <  2:176891/100‑1 (MQ=255)

                                                               gCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGa                                 >  1:1017248/1‑100 (MQ=255)

                                                                  aaaGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGaaaa                              >  1:3212935/1‑100 (MQ=255)

                                                                  aaaGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGaaaa                              <  1:32756/100‑1 (MQ=255)

                                                                   aaGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAAt                             <  1:2929109/100‑1 (MQ=255)

                                                                   aaGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAAt                             <  1:3112243/100‑1 (MQ=255)

                                                                   aaGCTAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAAt                             <  1:36152/100‑1 (MQ=255)

                                                                       tAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAAc                         >  1:1425404/1‑100 (MQ=255)

                                                                       tAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAAc                         >  1:1425333/1‑100 (MQ=255)

                                                                       tAACACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAAc                         >  1:2633544/1‑100 (MQ=255)

                                                                        aaCACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAAcc                        >  2:2848293/1‑100 (MQ=255)

                                                                        aaCACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGACAAAAAAAAAAAAAAATCGCCAAGAAAATTAAcc                        <  2:1517358/100‑1 (MQ=255)

                                                                        aaCACAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTAGCCACATCCATTGAAAAATCGCCAAGAAAATTAAcc                        >  1:1074081/1‑100 (MQ=255)

                                                                         acacAGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCa                       <  1:2344794/100‑1 (MQ=255)

                                                                          cacaGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAg                      <  2:2757235/100‑1 (MQ=255)

                                                                          cacaGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAg                      <  2:1579778/100‑1 (MQ=255)

                                                                          cacaGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAg                      >  2:836789/1‑100 (MQ=255)

                                                                             aGGCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGAtc                   <  2:2157364/100‑1 (MQ=255)

                                                                              ggCAAAAGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGAtct                  >  2:3136195/1‑100 (MQ=255)

                                                                                 aaaaGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGc               <  2:1188785/100‑1 (MQ=255)

                                                                                 aaaaGAAGCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGc               <  1:735153/100‑1 (MQ=255)

                                                                                    agaagCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGctct            <  1:90051/100‑1 (MQ=255)

                                                                                     gaagCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTg           <  2:2153590/100‑1 (MQ=255)

                                                                                     gaagCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTg           <  1:2502048/100‑1 (MQ=255)

                                                                                     gaagCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTg           <  1:3312548/100‑1 (MQ=255)

                                                                                       agCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGaa         >  2:1694497/1‑100 (MQ=255)

                                                                                        gCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAAc        <  1:2242657/100‑1 (MQ=255)

                                                                                        gCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAAc        >  2:662682/1‑100 (MQ=255)

                                                                                        gCCAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAAc        >  1:1329636/1‑100 (MQ=255)

                                                                                          cAGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTc      >  1:2411230/1‑100 (MQ=255)

                                                                                           aGACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCa     <  1:3016806/100‑1 (MQ=255)

                                                                                            gACTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTTGCTCTGAACTCat    <  1:2286573/100‑1 (MQ=255)

                                                                                              cTTTCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCataa  >  1:2014451/1‑100 (MQ=255)

TAAAACAACTGTACACACTATCCCTGTCAGTAAAATTGCTCCTCCGACAAGAATTACTGGAAAGCCAAAGCTAACACAGGCAAAAGAAGCCAGACTTCCAATAGCGGTTGCAACACCACCGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAA  >  NC_000913/238893‑239086

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