BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,616,532	G→A	T134I (ACC→ATC)	ydeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,532	0	G	A	85.7%	711.3 / 73.3	252	T134I (ACC→ATC)	ydeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (19/17); new base A (111/105); total (130/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-01

Only 100 of 256 total aligned reads displayed.

AAAAGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCC  >  NC_000913/1616440‑1616624

aaaaGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATgg                                                                                              >  2:5840023/1‑93 (MQ=255)

  aaGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGtt                                                                                            >  2:4523166/1‑93 (MQ=255)

  aaGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGAtt                                                                                            >  1:1602069/1‑93 (MQ=255)

  aaGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGAtt                                                                                            >  2:2240689/1‑93 (MQ=255)

   aGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTc                                                                                           >  1:1995491/1‑93 (MQ=255)

    gATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTcc                                                                                          <  2:1901740/93‑1 (MQ=255)

    gATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTcc                                                                                          <  2:7470123/93‑1 (MQ=255)

    gATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTcc                                                                                          >  1:2984326/1‑93 (MQ=255)

     aTCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTccc                                                                                         <  1:3499453/93‑1 (MQ=255)

     aTCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTccc                                                                                         <  2:6697709/93‑1 (MQ=255)

      tCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCg                                                                                        >  1:1487016/1‑93 (MQ=255)

      tCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCg                                                                                        >  1:6512282/1‑93 (MQ=255)

      tCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCg                                                                                        >  1:3123454/1‑93 (MQ=255)

       cAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCgg                                                                                       >  2:4573209/1‑93 (MQ=255)

         gcggcggcgCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGAc                                                                                     <  2:598544/93‑1 (MQ=255)

          cggcggcgCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGAcc                                                                                    <  2:1556817/93‑1 (MQ=255)

            gcggcgCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCgg                                                                                  >  2:267386/1‑93 (MQ=255)

             cggcgCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGc                                                                                 <  2:5673063/93‑1 (MQ=255)

              ggcgCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCa                                                                                >  1:364079/1‑93 (MQ=255)

              ggcgCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCa                                                                                >  1:5124845/1‑93 (MQ=255)

               gcgcAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAg                                                                               <  1:1008518/93‑1 (MQ=255)

                  cAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTAc                                                                            >  2:6441291/1‑93 (MQ=255)

                  cAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTAc                                                                            >  1:4096017/1‑93 (MQ=255)

                    aCCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACtt                                                                          >  2:2826702/1‑93 (MQ=255)

                    aCCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACtt                                                                          <  2:5065344/93‑1 (MQ=255)

                    aCCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACtt                                                                          <  2:3651958/93‑1 (MQ=255)

                    aCCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACtt                                                                          <  2:3383193/93‑1 (MQ=255)

                    aCCATCAGTACCAAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACtt                                                                          >  1:5769032/1‑93 (MQ=255)

                       aTCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGc                                                                       <  2:7070564/93‑1 (MQ=255)

                        tCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCa                                                                      >  1:5527340/1‑93 (MQ=255)

                          aGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCAtt                                                                    <  1:3312222/93‑1 (MQ=255)

                           gTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTg                                                                   >  2:3470664/1‑93 (MQ=255)

                           gTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTg                                                                   >  2:1644264/1‑93 (MQ=255)

                           gTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTg                                                                   >  2:3336027/1‑93 (MQ=255)

                            tACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGc                                                                  <  1:7465310/93‑1 (MQ=255)

                             aCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGcc                                                                 <  1:1931328/93‑1 (MQ=255)

                              cccAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCa                                                                <  2:2732021/93‑1 (MQ=255)

                               ccAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAg                                                               <  1:6965640/93‑1 (MQ=255)

                                 aaTCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCg                                                             >  1:4324221/1‑93 (MQ=255)

                                  aTCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCgg                                                            >  1:7355126/1‑93 (MQ=255)

                                   tCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGc                                                           <  1:1336617/93‑1 (MQ=255)

                                   tCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGc                                                           >  2:1374921/1‑93 (MQ=255)

                                    cAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGca                                                          <  2:4913017/93‑1 (MQ=255)

                                    cAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGca                                                          <  2:2604339/93‑1 (MQ=255)

                                      gCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGcaca                                                        >  1:3259849/1‑93 (MQ=255)

                                        gAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACaaa                                                      >  1:2291882/1‑93 (MQ=255)

                                        gAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACaaa                                                      >  2:1474010/1‑93 (MQ=255)

                                         aaGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAg                                                     <  2:2435431/93‑1 (MQ=255)

                                         aaGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAg                                                     >  2:589480/1‑93 (MQ=255)

                                         aaGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAg                                                     >  2:3871507/1‑93 (MQ=255)

                                            gTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCa                                                  >  2:2222819/1‑93 (MQ=255)

                                             ttGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCaa                                                 <  1:6559422/93‑1 (MQ=255)

                                                 aCTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATAtt                                             <  1:241306/93‑1 (MQ=255)

                                                  ctctGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTg                                            <  2:5811405/93‑1 (MQ=255)

                                                  ctctGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTg                                            >  2:5975634/1‑93 (MQ=255)

                                                   tctGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTgg                                           >  1:6051614/1‑93 (MQ=255)

                                                   tctGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTgg                                           <  1:2187880/93‑1 (MQ=255)

                                                    ctGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGc                                          <  1:6224169/93‑1 (MQ=255)

                                                     tGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCa                                         <  2:937600/93‑1 (MQ=255)

                                                     tGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCa                                         <  2:3943847/93‑1 (MQ=255)

                                                     tGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCa                                         <  1:1137259/93‑1 (MQ=255)

                                                         gCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCAGGACCGGCAGTACTTGGCATTGCCAGCGGCGCAAAGGCAATATTGGCCCTgg                                     <  1:1155262/93‑1 (MQ=255)

                                                         gCACTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTgg                                     <  2:3056321/93‑1 (MQ=255)

                                                           aCTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGt                                   >  1:3882438/1‑93 (MQ=255)

                                                           aCTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGt                                   <  2:6188670/93‑1 (MQ=255)

                                                            cTGTGGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGtt                                  >  1:3955651/1‑93 (MQ=255)

                                                             tgtgGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTc                                 >  2:2882776/1‑93 (MQ=255)

                                                              gtgGACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCa                                <  2:1222477/93‑1 (MQ=255)

                                                                 gACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCt                             <  1:4113459/93‑1 (MQ=255)

                                                                 gACGCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCt                             <  1:6253932/93‑1 (MQ=255)

                                                                  aCGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCtt                            >  2:4404739/1‑93 (MQ=255)

                                                                    gCCGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTcc                          <  2:83633/93‑1 (MQ=255)

                                                                     ccGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCa                         <  2:3157227/93‑1 (MQ=255)

                                                                      cGAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAg                        >  1:3726698/1‑93 (MQ=255)

                                                                       gAGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGc                       <  1:755505/93‑1 (MQ=255)

                                                                        aGCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGct                      <  1:976850/93‑1 (MQ=255)

                                                                         gCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGctc                     <  2:1582735/93‑1 (MQ=255)

                                                                         gCTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGctc                     >  2:6764942/1‑93 (MQ=255)

                                                                          cTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGctct                    <  2:342397/93‑1 (MQ=255)

                                                                          cTAATAATCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGctct                    <  2:4803300/93‑1 (MQ=255)

                                                                           taataaTCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCtt                   <  1:6704150/93‑1 (MQ=255)

                                                                           taataaTCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCtt                   >  1:5882809/1‑93 (MQ=255)

                                                                            aataatCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTc                  >  1:282711/1‑93 (MQ=255)

                                                                             ataatCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTcc                 >  2:6261115/1‑93 (MQ=255)

                                                                             ataatCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTcc                 <  1:6106710/93‑1 (MQ=255)

                                                                              taatCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCg                >  1:1067028/1‑93 (MQ=255)

                                                                              taatCATCGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCg                >  1:457481/1‑93 (MQ=255)

                                                                                 tcatcGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACt             >  1:1710522/1‑93 (MQ=255)

                                                                                   atcGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTg           >  1:2789258/1‑93 (MQ=255)

                                                                                   atcGCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTg           >  2:7161026/1‑93 (MQ=255)

                                                                                    tcGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGc          >  1:2736952/1‑93 (MQ=255)

                                                                                      gCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCtt        >  1:6770851/1‑93 (MQ=255)

                                                                                      gCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCtt        >  1:5289046/1‑93 (MQ=255)

                                                                                      gCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCtt        >  1:713623/1‑93 (MQ=255)

                                                                                      gCTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCtt        >  1:1792611/1‑93 (MQ=255)

                                                                                       cTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCttt       >  2:3707522/1‑93 (MQ=255)

                                                                                       cTATGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCCCTTCCGACTTGCctt       >  1:3634277/1‑93 (MQ=255)

                                                                                         aTGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTc     <  2:5481868/93‑1 (MQ=255)

                                                                                         aTGATTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTc     <  2:1989086/93‑1 (MQ=255)

                                                                                            gTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGcc  >  2:2112984/1‑93 (MQ=255)

AAAAGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCC  >  NC_000913/1616440‑1616624

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