BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,362,769	G→A	E66K (GAA→AAA)	yfaO →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,362,769	0	G	A	99.5%	713.3 / ‑5.3	201	E66K (GAA→AAA)	yfaO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (117/83); total (118/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 209 total aligned reads displayed.

CGGGTGGCGGCGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAGAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGATT  >  NC_000913/2362680‑2362858

cGGGAGGCGGCGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGCACTGGAAGAACAGCTCCTTTTGACAaaaa                                                                                        <  2:924596/93‑1 (MQ=255)

 gggTGGCGGCGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAAt                                                                                       >  1:3930897/1‑93 (MQ=255)

 gggTGGCGGCGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAAt                                                                                       >  1:3193005/1‑93 (MQ=255)

 gggTGGCGGCGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAAt                                                                                       <  2:4632591/93‑1 (MQ=255)

 gggTGGCGGCGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAAt                                                                                       >  1:3378584/1‑93 (MQ=255)

 gggTGGCGGCGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAAt                                                                                       <  2:2140/93‑1 (MQ=255)

    tGGCGGCGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCAc                                                                                    <  2:2380289/93‑1 (MQ=255)

     ggcggcGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACg                                                                                   <  1:1044082/93‑1 (MQ=255)

      gcggcgTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGc                                                                                  <  1:1277592/93‑1 (MQ=255)

       cggcgTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAAAAGCTGCTTTTGACAAAAATCACGcc                                                                                 >  2:2185068/1‑93 (MQ=255)

        ggcgTGGAGCGTGGCGAACGAATTGAAGAGTCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCg                                                                                <  2:1762745/93‑1 (MQ=255)

        ggcgTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCg                                                                                >  2:297781/1‑93 (MQ=255)

        ggcgTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCg                                                                                <  1:2876600/93‑1 (MQ=255)

         gcgTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGt                                                                               >  1:3438836/1‑93 (MQ=255)

           gTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTgg                                                                             >  2:4276470/1‑93 (MQ=255)

            tGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGa                                                                            <  2:3333466/93‑1 (MQ=255)

            tGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGa                                                                            >  1:646228/1‑93 (MQ=255)

               aGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCt                                                                         >  1:3517736/1‑93 (MQ=255)

               aGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCt                                                                         >  1:2439552/1‑93 (MQ=255)

                gCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCtt                                                                        >  2:1930589/1‑93 (MQ=255)

                  cTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCa                                                                      <  1:56088/93‑1 (MQ=255)

                  cTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCa                                                                      <  1:3228514/93‑1 (MQ=255)

                  cTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCa                                                                      >  2:2887019/1‑93 (MQ=255)

                    ggCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGc                                                                    <  1:3039829/93‑1 (MQ=255)

                    ggCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGc                                                                    >  2:1257149/1‑93 (MQ=255)

                     gCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAGAAATCACGCCGTGGACCTTCAGCg                                                                   >  1:4138112/1‑93 (MQ=255)

                     gCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCg                                                                   >  1:90349/1‑93 (MQ=255)

                     gCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCg                                                                   >  2:2687243/1‑93 (MQ=255)

                      cgaacggaTTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCga                                                                  >  1:3163036/1‑93 (MQ=255)

                        aacgaaTTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCgatg                                                                >  1:388647/1‑93 (MQ=255)

                          cgaaTTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCgatgat                                                              >  2:4376903/1‑93 (MQ=255)

                           gaaTTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGTTGata                                                             >  2:1629968/1‑93 (MQ=255)

                           gaaTTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGata                                                             >  1:750454/1‑93 (MQ=255)

                            aaTTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGGGGACCTTCAGCGATGatat                                                            >  1:4615807/1‑93 (MQ=255)

                             aTTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATAtt                                                           <  2:4680514/93‑1 (MQ=255)

                              ttGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTc                                                          >  2:2938889/1‑93 (MQ=255)

                                     ggCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCaa                                                   >  1:364786/1‑93 (MQ=255)

                                       cACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGa                                                 <  1:1379606/93‑1 (MQ=255)

                                       cACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGa                                                 >  2:3637823/1‑93 (MQ=255)

                                       cACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGa                                                 <  1:2034830/93‑1 (MQ=255)

                                       cACTACGCAGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCCCGACGTGGACCTTCAGCGATGATATTCGCACCAAGa                                                 <  1:1201437/93‑1 (MQ=255)

                                        aCTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGAc                                                <  1:3799966/93‑1 (MQ=255)

                                          tACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGt                                              >  1:3784089/1‑93 (MQ=255)

                                          tACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGt                                              >  2:3079465/1‑93 (MQ=255)

                                           aCGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTa                                             <  2:2542622/93‑1 (MQ=255)

                                           aCGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTa                                             >  2:2744036/1‑93 (MQ=255)

                                            cgccgcGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTAt                                            <  1:1250065/93‑1 (MQ=255)

                                            cgccgcGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTAt                                            <  2:1145025/93‑1 (MQ=255)

                                             gccgcGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATg                                           >  2:2242878/1‑93 (MQ=255)

                                             gccgcGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATCCGCACCAAGACGTATg                                           >  2:1565948/1‑93 (MQ=255)

                                              ccgcGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGc                                          >  2:4230796/1‑93 (MQ=255)

                                                 cgAAATTCGCGAAGAACTGGGAGAACAGTGACTGTGGACAAAGATTACGCCGTGGACTGTCAGCGATGATATTCGATTCAGGATGGATGCAGa                                       <  2:2717647/93‑1 (MQ=255)

                                                 cgAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGa                                       <  1:1589767/93‑1 (MQ=255)

                                                  gAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCTCACCAAGACGTATGCAGAt                                      >  2:605989/1‑93 (MQ=255)

                                                  gAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGAt                                      <  1:2793765/93‑1 (MQ=255)

                                                  gAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGAt                                      <  2:3885029/93‑1 (MQ=255)

                                                  gAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGAt                                      <  2:4719756/93‑1 (MQ=255)

                                                  gAAATTCGAGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGAt                                      >  1:488579/1‑93 (MQ=255)

                                                  gAAAATGGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGAt                                      <  1:2134938/93‑1 (MQ=255)

                                                   aaaTTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATg                                     >  2:1215626/1‑93 (MQ=255)

                                                    aaTTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATgg                                    >  2:4615647/1‑93 (MQ=255)

                                                     aTTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGt                                   >  1:4542641/1‑93 (MQ=255)

                                                      ttCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTc                                  <  2:475193/93‑1 (MQ=255)

                                                       tCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCg                                 >  1:2674054/1‑93 (MQ=255)

                                                          cgAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCaa                              >  1:2715821/1‑93 (MQ=255)

                                                          cgAAGAACTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCaa                              >  1:2808630/1‑93 (MQ=255)

                                                           gaagaaCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATCTTCGCACCAAGACGTATGCAGATGGTCGCAAg                             >  2:1644407/1‑93 (MQ=255)

                                                           gaagaaCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAg                             <  1:3391266/93‑1 (MQ=255)

                                                            aagaaCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAgg                            >  2:3493814/1‑93 (MQ=255)

                                                            aagaaCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAgg                            >  1:560009/1‑93 (MQ=255)

                                                             agaaCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGa                           >  2:3422492/1‑93 (MQ=255)

                                                              gaaCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGaa                          >  2:1284740/1‑93 (MQ=255)

                                                              gaaCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGaa                          >  1:2407863/1‑93 (MQ=255)

                                                              gaaCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGaa                          >  2:1783940/1‑93 (MQ=255)

                                                               aaCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAag                         >  1:398798/1‑93 (MQ=255)

                                                                aCTGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAaga                        <  1:2395607/93‑1 (MQ=255)

                                                                  tGGGAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAgaga                      >  2:1929484/1‑93 (MQ=255)

                                                                   gggAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGAt                     >  2:1464017/1‑93 (MQ=255)

                                                                   gggAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGAt                     >  2:3595965/1‑93 (MQ=255)

                                                                   gggAGAACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGAt                     >  1:3916317/1‑93 (MQ=255)

                                                                      agaACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTta                  <  2:950049/93‑1 (MQ=255)

                                                                       gaACAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTtat                 >  1:2957595/1‑93 (MQ=255)

                                                                        aaCAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTtata                <  2:2830695/93‑1 (MQ=255)

                                                                        aaCAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTtata                >  1:1644677/1‑93 (MQ=255)

                                                                        aaCAGCTGCTTTTGACAAAAATCACACCGTGGACCTTCAGCGATGATATTCGCACCCAGACGTATGCAGATGGTCGCCAGGAAGAGATTtaga                >  2:53490/1‑91 (MQ=255)

                                                                        aaCAGCTGCTTTTGACAAAAATCACACCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTtata                >  2:53492/1‑93 (MQ=255)

                                                                         aCAGCTGCTTTTGACAAAAATCATGCCGTGGACCTTTAGCGATGATATTCGCACCATGACGTCTGCAGAAGGACGCTAGGAAGAGACatatat               >  2:2352895/1‑93 (MQ=255)

                                                                         aCAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTatat               >  1:4121205/1‑93 (MQ=255)

                                                                         aCAGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTatat               >  2:768013/1‑93 (MQ=255)

                                                                           aGCTGCTTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGa             <  1:1903303/93‑1 (MQ=255)

                                                                            gctgcgTTTGACAAAAATCACGCCGTGGACCTGCAGCGATGATATTCGCACCAAGACGTATGCAGATGGGCGCAAGGAAGAGATTTATATGAt            >  1:4203968/1‑93 (MQ=255)

                                                                             ctgctTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGAtt           <  1:2296341/93‑1 (MQ=255)

                                                                              tgctTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCGAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGAttt          >  1:1766394/1‑93 (MQ=255)

                                                                              tgctTTTGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGAttt          <  1:3262088/93‑1 (MQ=255)

                                                                                 ttttGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTAcc       <  2:4183479/93‑1 (MQ=255)

                                                                                  tgtGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCt      <  2:3343168/91‑1 (MQ=255)

                                                                                   ttGACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTg     <  1:804469/93‑1 (MQ=255)

                                                                                     gACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGAt   <  1:117570/93‑1 (MQ=255)

                                                                                     gACAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGAt   >  1:4400971/1‑93 (MQ=255)

                                                                                      aCAAAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGAtt  >  2:3061045/1‑93 (MQ=255)

CGGGTGGCGGCGTGGAGCCTGGCGAACGAATTGAAGAGGCACTACGCCGCGAAATTCGCGAAGAACTGGGAGAACAGCTGCTTTTGACAGAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGATT  >  NC_000913/2362680‑2362858

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