BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,566	A→G	G515G (GGT→GGC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,566	0	A	G	99.0%	291.5 / NA	95	G515G (GGT→GGC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (46/48); minor base C (0/1); total (46/49)
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.99e-01

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAG  >  NC_000913/1594476‑1594655

aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                                                                         >  2:43319/1‑93 (MQ=255)

aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                                                                         >  1:191015/1‑93 (MQ=255)

 ttCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCa                                                                                        <  2:1746635/93‑1 (MQ=255)

   ccAGAGCCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAtt                                                                                      <  2:392494/93‑1 (MQ=255)

    cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                                                                     <  1:2260723/93‑1 (MQ=255)

    cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                                                                     <  1:2680223/93‑1 (MQ=255)

     agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTGGAATTAAGGGTGGCGCCATTat                                                                                    >  1:45645/1‑92 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAGATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 <  1:2044224/93‑1 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 >  1:2974936/1‑93 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                 >  2:4676393/1‑93 (MQ=255)

         aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGAc                                                                                <  1:3962506/93‑1 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                               <  1:4067425/93‑1 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                               >  2:3566373/1‑93 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                               >  1:1313277/1‑93 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                               >  2:489082/1‑93 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                               <  2:3253978/93‑1 (MQ=255)

           tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                                                              >  2:3319087/1‑93 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                             <  1:4373634/93‑1 (MQ=255)

             tCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAttt                                                                            >  1:4084283/1‑93 (MQ=255)

                gTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                                                                         >  1:4600901/1‑93 (MQ=255)

                gTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                                                                         >  1:1913978/1‑93 (MQ=255)

                 ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                        >  1:3736504/1‑93 (MQ=255)

                   gAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATa                                                                      >  1:1138966/1‑93 (MQ=255)

                    aaaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAtcc                                                                     <  2:417440/93‑3 (MQ=255)

                    aaaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAc                                                                     <  1:1664774/93‑1 (MQ=255)

                     aaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAcc                                                                    <  1:2497056/93‑1 (MQ=255)

                      atatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCg                                                                   >  2:3221804/1‑93 (MQ=255)

                       tatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                                                                  <  1:206226/93‑1 (MQ=255)

                       tatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                                                                  <  2:3833627/93‑1 (MQ=255)

                             cTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaa                                                            >  2:688162/1‑93 (MQ=255)

                               aTCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaaa                                                          <  2:2603015/93‑1 (MQ=255)

                                tCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt                                                         >  1:3877447/1‑93 (MQ=255)

                                 ccccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAtt                                                        >  1:4328961/1‑93 (MQ=255)

                                 ccccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAtt                                                        <  1:4631823/93‑1 (MQ=255)

                                  cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa                                                       >  2:3107321/1‑93 (MQ=255)

                                  cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa                                                       >  2:1067318/1‑93 (MQ=255)

                                   ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                      >  1:3933511/1‑93 (MQ=255)

                                    ctgaCCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTTACATTTACGATACCGCTACCAAAATTAcc                                                     <  2:3785434/90‑1 (MQ=255)

                                    cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                     >  2:3957813/1‑93 (MQ=255)

                                    cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                     <  1:958856/93‑1 (MQ=255)

                                    cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCAATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                     >  1:1988760/1‑93 (MQ=255)

                                       aCCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc                                                  >  1:639633/1‑93 (MQ=255)

                                        ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                 <  1:1095757/93‑1 (MQ=255)

                                        ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                 <  2:2321563/93‑1 (MQ=255)

                                         cGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAt                                                <  1:2565018/93‑1 (MQ=255)

                                          gggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                               <  2:1119666/90‑1 (MQ=255)

                                             ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                            <  1:4556486/93‑1 (MQ=255)

                                             ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                            >  2:3483715/1‑93 (MQ=255)

                                             ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                            <  2:1017143/93‑1 (MQ=255)

                                               cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAc                                          <  1:4046319/93‑1 (MQ=255)

                                               cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAc                                          <  1:1753061/93‑1 (MQ=255)

                                                 tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                        >  2:598576/1‑93 (MQ=255)

                                                  ttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCca                                       <  2:268770/93‑1 (MQ=255)

                                                   tcctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                      <  2:2105263/93‑1 (MQ=255)

                                                    cctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                     >  2:927881/1‑93 (MQ=255)

                                                    cctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                     >  1:1229646/1‑93 (MQ=255)

                                                     ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAg                                    >  2:1793568/1‑93 (MQ=255)

                                                     ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAg                                    >  2:1998934/1‑93 (MQ=255)

                                                      tccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGt                                   <  2:1318360/93‑1 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                 >  1:2366005/1‑93 (MQ=255)

                                                         gATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                                >  1:4087215/1‑93 (MQ=255)

                                                           tAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGa                              <  1:3592074/93‑1 (MQ=255)

                                                            aaaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGaa                             <  2:2277707/93‑1 (MQ=255)

                                                            aaaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGaa                             <  2:2364394/93‑1 (MQ=255)

                                                            aaaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGaa                             >  1:3139756/1‑93 (MQ=255)

                                                             aaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAc                            <  2:2918262/93‑1 (MQ=255)

                                                             aaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAc                            >  1:2631110/1‑93 (MQ=255)

                                                             aaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAc                            <  2:3406900/93‑1 (MQ=255)

                                                              aTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACACACAGTATCGAACg                           <  1:4351120/93‑1 (MQ=255)

                                                               tCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACGGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                          <  1:2506337/93‑1 (MQ=255)

                                                               tCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                          <  1:1191892/93‑1 (MQ=255)

                                                                  acaGCCCGTTGAATTAAGTGTGGCCCCGTTACTGACATTTACGATACCGCTACCAAAATTCCCGACATTGATACCCACAGTATCGAACgatga                       <  2:1367469/91‑1 (MQ=255)

                                                                   tAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGAGATTGATACCCACAGTATCGAACgatgat                      >  1:4509468/1‑93 (MQ=255)

                                                                   tAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatgat                      >  1:2209695/1‑93 (MQ=255)

                                                                   tAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatgat                      <  2:1975145/93‑1 (MQ=255)

                                                                   tAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatgat                      <  1:250587/93‑1 (MQ=255)

                                                                     gCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGgtgt                    >  2:2193600/1‑93 (MQ=255)

                                                                     gCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGt                    <  1:1252687/93‑1 (MQ=255)

                                                                     gCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGt                    >  1:3460003/1‑93 (MQ=255)

                                                                      cccGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTc                   >  2:1834734/1‑93 (MQ=255)

                                                                        cGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGc                 <  1:1742323/93‑1 (MQ=255)

                                                                         gTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGca                >  2:3116384/1‑93 (MQ=255)

                                                                         gTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGca                <  2:292436/93‑1 (MQ=255)

                                                                          ttGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcat               >  2:4675979/1‑93 (MQ=255)

                                                                            gAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACACACAGTATCGAACGATGATGTCGcatca             >  1:3390675/1‑92 (MQ=255)

                                                                             aaTTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcca            >  2:989231/1‑93 (MQ=255)

                                                                             aaTTAAGGGGGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcca            <  1:1637065/93‑1 (MQ=255)

                                                                              aTTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccat           >  2:3645196/1‑93 (MQ=255)

                                                                              aTTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccat           <  1:1324793/93‑1 (MQ=255)

                                                                              aTTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccat           >  2:2704485/1‑93 (MQ=255)

                                                                              aTTAAGGGTGGCGCCATAACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATCCCCACAGTATCGAACGATGATGTCGcatccat           <  2:2404593/93‑1 (MQ=255)

                                                                                tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa         <  2:2853546/93‑1 (MQ=255)

                                                                                 aaGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAg        >  2:3534778/1‑93 (MQ=255)

                                                                                  aGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAgg       <  1:1760486/93‑1 (MQ=255)

                                                                                   gggTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGt      >  1:713740/1‑93 (MQ=255)

                                                                                   gggTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGt      >  1:4140230/1‑93 (MQ=255)

                                                                                    ggTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGtt     <  2:852844/93‑1 (MQ=255)

                                                                                      tGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCa   <  1:95748/93‑1 (MQ=255)

                                                                                      tGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCa   <  1:255900/93‑1 (MQ=255)

                                                                                       ggCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAg  <  1:2844279/93‑1 (MQ=255)

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAG  >  NC_000913/1594476‑1594655

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