BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,129,641	C→T	V261V (GTG→GTA)	b2983 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,129,641	0	C	T	98.1%	529.9 / ‑2.4	157	V261V (GTG→GTA)	b2983	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (2/1); new base T (75/78); total (77/80)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 161 total aligned reads displayed.

CAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGC  >  NC_000913/3129551‑3129730

cAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAatac                                                                                         >  1:2470573/1‑93 (MQ=255)

  accccGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACtg                                                                                       <  2:903143/91‑1 (MQ=255)

  aaCCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACtg                                                                                       >  1:1813804/1‑93 (MQ=255)

   aCCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGCTTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACtgt                                                                                      >  1:2046980/1‑93 (MQ=255)

    cccGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACtgtg                                                                                     <  2:867938/93‑1 (MQ=255)

     ccGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgc                                                                                    <  2:3269492/93‑1 (MQ=255)

     ccGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgc                                                                                    <  1:552693/93‑1 (MQ=255)

     ccGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgc                                                                                    >  1:3367444/1‑93 (MQ=255)

      cGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgct                                                                                   >  1:723704/1‑93 (MQ=255)

       gCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgctg                                                                                  >  1:3135267/1‑93 (MQ=255)

        cTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgctgc                                                                                 >  2:2504156/1‑93 (MQ=255)

        cTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGtgctgc                                                                                 >  1:3488084/1‑93 (MQ=255)

          ggCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCAc                                                                               >  2:2555349/1‑93 (MQ=255)

          ggCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCAc                                                                               <  1:2956299/93‑1 (MQ=255)

                ccGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTa                                                                         <  1:4425610/93‑1 (MQ=255)

                 cggcggCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTAc                                                                        <  1:1011862/93‑1 (MQ=255)

                 cggcggCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTAc                                                                        >  1:2227653/1‑93 (MQ=255)

                   gcggagTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTAcgc                                                                      >  1:843095/1‑93 (MQ=255)

                      gcgTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGcc                                                                   <  2:3750871/93‑1 (MQ=255)

                       cgTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCg                                                                  <  2:2077999/93‑1 (MQ=255)

                         tCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGAc                                                                >  2:1431322/1‑93 (MQ=255)

                           gaagaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCa                                                              <  2:2703596/93‑1 (MQ=255)

                           gaagaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCa                                                              <  2:541901/93‑1 (MQ=255)

                             agaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGa                                                            <  2:2889138/93‑1 (MQ=255)

                             agaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGa                                                            <  2:2889139/93‑1 (MQ=255)

                              gaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACCCGCCGACCAGag                                                           >  1:3582212/1‑92 (MQ=255)

                              gaaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGAt                                                           >  1:2122944/1‑93 (MQ=255)

                               aaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATg                                                          <  2:1576933/93‑1 (MQ=255)

                               aaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATg                                                          >  1:2435514/1‑93 (MQ=255)

                               aaTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATg                                                          >  1:3995486/1‑93 (MQ=255)

                                aTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATgg                                                         >  2:2619278/1‑93 (MQ=255)

                                 tGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGa                                                        <  2:3496489/93‑1 (MQ=255)

                                 tGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGa                                                        <  1:1942758/93‑1 (MQ=255)

                                  gTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAg                                                       >  1:1534231/1‑93 (MQ=255)

                                  gTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAg                                                       <  1:1618863/93‑1 (MQ=255)

                                  gTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAg                                                       <  1:1715358/93‑1 (MQ=255)

                                    cATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTg                                                     >  1:2877028/1‑93 (MQ=255)

                                     aTGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTgg                                                    <  2:3330185/93‑1 (MQ=255)

                                     aTGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTgg                                                    >  1:1720471/1‑93 (MQ=255)

                                      tGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTggg                                                   <  2:4301583/93‑1 (MQ=255)

                                       ggCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGc                                                  >  2:3691428/1‑93 (MQ=255)

                                          gATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAAt                                               <  1:1890059/93‑1 (MQ=255)

                                          gATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAAt                                               <  1:754809/93‑1 (MQ=255)

                                               tAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGGGGGCAATgttgg                                          >  1:1589614/1‑93 (MQ=255)

                                                 aaTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATgttggtt                                        >  1:3657708/1‑93 (MQ=255)

                                                  aTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATgttggttg                                       >  1:3441896/1‑93 (MQ=255)

                                                  aTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATgttggttg                                       <  1:2311431/93‑1 (MQ=255)

                                                   tAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGa                                      <  2:2690541/93‑1 (MQ=255)

                                                    aaCCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGAc                                     <  2:4501292/93‑1 (MQ=255)

                                                    aaCCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGAc                                     >  1:3860952/1‑93 (MQ=255)

                                                       ccGGCGGCAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGCTGGTTGACCAg                                  >  1:4101722/1‑93 (MQ=255)

                                                        cggcggCAGCGTGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGa                                 <  1:890335/93‑1 (MQ=255)

                                                          gcggcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCCGACa                               >  1:2927060/1‑93 (MQ=255)

                                                          gcggcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACa                               <  2:1755489/93‑1 (MQ=255)

                                                             gcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACaaaa                            >  1:1646740/1‑93 (MQ=255)

                                                             gcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACaaaa                            >  1:2196479/1‑93 (MQ=255)

                                                             gcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACaaaa                            >  2:2862789/1‑93 (MQ=255)

                                                             gcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACaaaa                            <  2:110023/93‑1 (MQ=255)

                                                             gcAGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACaaaa                            <  2:4166110/93‑1 (MQ=255)

                                                               aGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAAcc                          >  1:1206672/1‑92 (MQ=255)

                                                               aGCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACt                          >  1:4329223/1‑93 (MQ=255)

                                                                gCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTc                         >  1:3003269/1‑93 (MQ=255)

                                                                gCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTc                         >  1:895735/1‑93 (MQ=255)

                                                                gCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTc                         >  2:3416955/1‑93 (MQ=255)

                                                                gCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTc                         >  2:3570606/1‑93 (MQ=255)

                                                                gCGGGTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTc                         >  2:4616266/1‑93 (MQ=255)

                                                                  gggTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCa                       <  1:3164540/93‑1 (MQ=255)

                                                                  gggTCCTAACACGATCCCCAACAATACTGTGCTGCACGAGTTCCGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCa                       >  2:582467/1‑93 (MQ=255)

                                                                   ggTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAc                      <  1:3715605/93‑1 (MQ=255)

                                                                    gTCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAcg                     <  1:545486/93‑1 (MQ=255)

                                                                     tCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAcgc                    >  2:1891592/1‑93 (MQ=255)

                                                                     tCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAcgc                    <  1:4507049/93‑1 (MQ=255)

                                                                     tCCTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAcgc                    >  1:1140166/1‑93 (MQ=255)

                                                                       cTAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACgcgc                  >  1:1001641/1‑93 (MQ=255)

                                                                        tAACACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGcc                 >  2:3338887/1‑93 (MQ=255)

                                                                         aaCACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCt                >  2:863074/1‑93 (MQ=255)

                                                                         aaCACGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCt                <  1:4292182/93‑1 (MQ=255)

                                                                          acacGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCtt               <  2:1273142/93‑1 (MQ=255)

                                                                          acacGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCtt               <  1:4621506/93‑1 (MQ=255)

                                                                            acGATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTa             >  2:1024889/1‑93 (MQ=255)

                                                                             cGACCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTaa            <  2:4380126/93‑1 (MQ=255)

                                                                              gATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAAt           >  2:699042/1‑93 (MQ=255)

                                                                              gATCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAAt           >  1:2752281/1‑93 (MQ=255)

                                                                              gATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAAt           <  2:4697797/93‑1 (MQ=255)

                                                                               aTCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATg          <  1:262672/93‑1 (MQ=255)

                                                                               aTCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATg          >  2:2152319/1‑93 (MQ=255)

                                                                                tCCCCACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGt         >  1:1472920/1‑93 (MQ=255)

                                                                                 ccccACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTa        <  1:675979/93‑1 (MQ=255)

                                                                                  cccACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTAt       <  1:632122/93‑1 (MQ=255)

                                                                                  cccACCAATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTAt       <  1:3720642/93‑1 (MQ=255)

                                                                                   ccaccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATc      <  1:3315146/93‑1 (MQ=255)

                                                                                   ccaccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATc      <  1:1298901/93‑1 (MQ=255)

                                                                                   ccaccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATc      >  2:257340/1‑93 (MQ=255)

                                                                                    caccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCg     >  2:1718277/1‑93 (MQ=255)

                                                                                    caccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCg     >  1:3059338/1‑93 (MQ=255)

                                                                                    caccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCg     >  2:2665163/1‑93 (MQ=255)

                                                                                     cccaATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCgg    <  1:2748261/92‑1 (MQ=255)

                                                                                       caATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGc  >  2:4085797/1‑93 (MQ=255)

                                                                                       caATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGc  >  2:4227768/1‑93 (MQ=255)

                                                                                       caATACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGc  >  1:798728/1‑93 (MQ=255)

CAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGC  >  NC_000913/3129551‑3129730

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