BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,648,429	G→A	A110T (GCC→ACC)	ydfE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,648,429	0	G	A	95.6%	355.1 / 1.5	114	A110T (GCC→ACC)	ydfE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (3/2); new base A (50/58); total (53/61)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 118 total aligned reads displayed.

ATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAATTGTTG  >  NC_000913/1648339‑1648518

atGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGAcc                                                                                         >  1:1397467/1‑93 (MQ=255)

  gCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGAcccc                                                                                       <  1:1337402/93‑1 (MQ=255)

  gCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGAcccc                                                                                       >  1:4219522/1‑93 (MQ=255)

   cAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCt                                                                                      <  1:4173480/93‑1 (MQ=255)

    aGGAAAACCGCTGGGATGAGGTGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGCACGTCACTTTTGACAAGCTGCCGCCTg                                                                                     >  2:145735/1‑93 (MQ=255)

     ggAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGa                                                                                    <  1:3453643/93‑1 (MQ=255)

     ggAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGa                                                                                    >  1:3188714/1‑93 (MQ=255)

       aaaaCCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAg                                                                                  >  2:1962152/1‑93 (MQ=255)

         aaCCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTa                                                                                >  2:2702962/1‑93 (MQ=255)

         aaCCGCTGGGAGGAGGAGTCTGCTTCGTTTTTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTa                                                                                <  1:3596788/93‑1 (MQ=255)

           ccGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACa                                                                              >  1:4136295/1‑93 (MQ=255)

            cGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAg                                                                             >  2:932994/1‑93 (MQ=255)

             gCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGa                                                                            >  1:3684842/1‑93 (MQ=255)

              cTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGaa                                                                           >  1:295295/1‑93 (MQ=255)

              cTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGaa                                                                           <  1:4291064/93‑1 (MQ=255)

               tGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAAc                                                                          <  1:2681755/93‑1 (MQ=255)

                 ggATGAGGAGTCCGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGc                                                                        >  1:2533033/1‑93 (MQ=255)

                  gATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCt                                                                       >  2:485540/1‑93 (MQ=255)

                   aTGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTg                                                                      <  1:1320247/93‑1 (MQ=255)

                    tGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGt                                                                     <  1:1590657/93‑1 (MQ=255)

                    tGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGt                                                                     <  2:1428969/93‑1 (MQ=255)

                      aggagTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCa                                                                   <  2:424186/93‑1 (MQ=255)

                       ggagTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCAt                                                                  <  1:3652975/93‑1 (MQ=255)

                       ggagTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCAt                                                                  >  2:732659/1‑93 (MQ=255)

                       ggagTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCAt                                                                  <  1:1462527/93‑1 (MQ=255)

                        gagTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATg                                                                 <  1:3744034/93‑1 (MQ=255)

                        gagTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATg                                                                 <  1:1996825/93‑1 (MQ=255)

                          gTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGt                                                               >  1:4001977/1‑93 (MQ=255)

                          gTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGt                                                               <  1:3569925/93‑1 (MQ=255)

                           tCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGtt                                                              >  2:1260549/1‑93 (MQ=255)

                            cTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTa                                                             <  2:3411762/93‑1 (MQ=255)

                             tGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTaa                                                            <  1:3044983/93‑1 (MQ=255)

                             tGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTaa                                                            >  1:2314663/1‑93 (MQ=255)

                             tGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTaa                                                            <  2:2668560/93‑1 (MQ=255)

                              gCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAg                                                           <  1:689142/93‑1 (MQ=255)

                              gCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAg                                                           >  1:1718990/1‑93 (MQ=255)

                               cttcttTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGt                                                          <  1:3598177/93‑1 (MQ=255)

                                ttcttTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGtt                                                         <  1:672759/93‑1 (MQ=255)

                                 tcttTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTc                                                        >  2:4441523/1‑93 (MQ=255)

                                 tcttTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTc                                                        >  1:3013227/1‑93 (MQ=255)

                                 tcttTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTc                                                        >  2:132467/1‑93 (MQ=255)

                                   ttttGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGa                                                      >  2:2858538/1‑93 (MQ=255)

                                   ttttGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGa                                                      >  1:2047474/1‑93 (MQ=255)

                                    tttGTCCGGGTTGAGACTAAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGac                                                     <  2:345974/93‑1 (MQ=255)

                                    tttGTCCCGCTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGac                                                     <  2:2953835/93‑1 (MQ=255)

                                     ttGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGaca                                                    <  1:4351930/93‑1 (MQ=255)

                                       gTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGAcaca                                                  >  1:381789/1‑93 (MQ=255)

                                       gTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGAcaca                                                  <  2:251152/93‑1 (MQ=255)

                                        tCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACAt                                                 <  2:2089454/93‑1 (MQ=255)

                                         cccGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACAtg                                                <  2:4098660/93‑1 (MQ=255)

                                         cccGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACAtg                                                <  1:127928/93‑1 (MQ=255)

                                           cGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACAtgtg                                              >  1:273207/1‑93 (MQ=255)

                                           cGGTTGAGACTGAAGCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACAtgtg                                              <  1:3141239/93‑1 (MQ=255)

                                            ggTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGa                                             >  1:3709378/1‑93 (MQ=255)

                                            ggTTGAGACAGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGa                                             <  2:548035/93‑1 (MQ=255)

                                             gTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGaa                                            >  2:3486613/1‑93 (MQ=255)

                                               tGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGaaaa                                          >  2:669490/1‑93 (MQ=255)

                                                gagaCTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAAc                                         <  1:2242733/93‑1 (MQ=255)

                                                 agaCTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACa                                        <  2:1017657/93‑1 (MQ=255)

                                                 agaCTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACa                                        <  2:451876/93‑1 (MQ=255)

                                                 agaCTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACa                                        <  2:1641156/93‑1 (MQ=255)

                                                  gaCTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACAt                                       >  2:2235761/1‑93 (MQ=255)

                                                   aCTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATc                                      >  1:222670/1‑93 (MQ=255)

                                                   aCTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATc                                      <  1:2310180/93‑1 (MQ=255)

                                                   aCTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATc                                      <  1:1681698/93‑1 (MQ=255)

                                                    cTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCa                                     <  2:3008088/93‑1 (MQ=255)

                                                    cTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCa                                     <  2:1331820/93‑1 (MQ=255)

                                                    cTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCa                                     <  2:697073/93‑1 (MQ=255)

                                                    cTGAATCCAATCCGATGAACGTCGCTTTTGACAAGCTGCCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCa                                     <  1:1953772/93‑1 (MQ=255)

                                                     tGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCAc                                    >  1:2279241/1‑93 (MQ=255)

                                                     tGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCAc                                    <  2:2720560/93‑1 (MQ=255)

                                                     tGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCAc                                    <  1:4088807/93‑1 (MQ=255)

                                                       aaTCCGATCCGATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCg                                  <  1:1389655/93‑1 (MQ=255)

                                                        atccgatccgATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGt                                 <  2:4028691/93‑1 (MQ=255)

                                                        atccgatccgATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGt                                 >  2:3092405/1‑93 (MQ=255)

                                                        atccgatccgATGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGt                                 <  1:3934425/93‑1 (MQ=255)

                                                          ccgatccgatCAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTCAGTTCGACACATGTGAAAACATCACCGTTg                               >  2:94734/1‑93 (MQ=255)

                                                            gatccgatGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGat                             >  1:2183269/1‑93 (MQ=255)

                                                             atccgatGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGata                            >  1:1527883/1‑93 (MQ=255)

                                                               ccgatGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATg                          <  2:320467/93‑1 (MQ=255)

                                                                cgatGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATgg                         >  1:2461983/1‑93 (MQ=255)

                                                                   tGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGtta                      >  2:1207764/1‑93 (MQ=255)

                                                                   tGAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGtta                      <  1:256571/93‑1 (MQ=255)

                                                                    gAACGTCCCTTTTGACAGGCTGCCGCCTGAAGTACAGACCGCTGTCATGGTTAATTTCGACACATGTGAAAACATCTCCGTTGATATGGttat                     <  2:312818/93‑1 (MQ=255)

                                                                    gAACGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGttat                     <  2:3495906/93‑1 (MQ=255)

                                                                      aCGTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGttatta                   >  2:2627045/1‑93 (MQ=255)

                                                                        gTCACTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAgc                 <  2:1427248/93‑1 (MQ=255)

                                                                            cTTTTGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGcac             >  1:2383494/1‑93 (MQ=255)

                                                                             ttttGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGcaca            <  2:4423491/93‑1 (MQ=255)

                                                                             ttttGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGcaca            >  2:112010/1‑93 (MQ=255)

                                                                             ttttGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGcaca            <  2:4299066/93‑1 (MQ=255)

                                                                              tttGACAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAg           <  2:4550194/93‑1 (MQ=255)

                                                                                  aCAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAAt       >  1:1796577/1‑93 (MQ=255)

                                                                                  aCAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAAt       <  2:3929264/93‑1 (MQ=255)

                                                                                  aCAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAAt       <  1:1145566/93‑1 (MQ=255)

                                                                                  aCAAGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAAt       <  2:3091858/93‑1 (MQ=255)

                                                                                   cAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAAtt      >  1:2140837/1‑93 (MQ=255)

                                                                                    aaGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAAttg     >  2:1640575/1‑93 (MQ=255)

                                                                                     aGCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAAttgt    >  1:1625009/1‑93 (MQ=255)

                                                                                      gCTGACCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAAttgtt   >  2:2572217/1‑93 (MQ=255)

ATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTCGACACATGTGAAAACATCACCGTTGATATGGTTATTAGCGCACAGGAATTGTTG  >  NC_000913/1648339‑1648518

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