BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,895,750	T→C	C34R (TGT→CGT)	yieJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,895,750	0	T	C	85.0%	737.7 / 121.8	315	C34R (TGT→CGT)	yieJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (22/25); new base C (138/129); total (160/155)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 320 total aligned reads displayed.

AATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGT  >  NC_000913/3895660‑3895839

aaTATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCcgt                                                                                         >  1:208961/1‑93 (MQ=255)

aaTATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCcgt                                                                                         >  1:4085587/1‑93 (MQ=255)

aaTATCAGGCCGTTACCCCAATTCAAATATCACCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCcgt                                                                                         <  1:491861/93‑1 (MQ=255)

  tatCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGa                                                                                       >  1:158218/1‑93 (MQ=255)

   atCAGGCCGTTCCCCAAATTAAAATATCACCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGaa                                                                                      <  1:1636940/93‑1 (MQ=255)

   atCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGaa                                                                                      >  2:6595550/1‑93 (MQ=255)

   atCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGaa                                                                                      <  1:393984/93‑1 (MQ=255)

    tCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGaac                                                                                     <  1:6563818/93‑1 (MQ=255)

        gCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAacaaca                                                                                 <  2:7283990/93‑1 (MQ=255)

         ccGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAg                                                                                <  2:2960639/93‑1 (MQ=255)

          cGTTACCCCAATTCAAATATCACCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGa                                                                               <  1:1434/93‑1 (MQ=255)

           gTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGAc                                                                              >  1:5450806/1‑93 (MQ=255)

           gTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGAc                                                                              >  2:742185/1‑93 (MQ=255)

            ttACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACg                                                                             >  2:1930644/1‑93 (MQ=255)

            ttACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACg                                                                             <  1:6348891/93‑1 (MQ=255)

               ccccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCa                                                                          >  1:2450920/1‑93 (MQ=255)

               ccccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCa                                                                          >  2:3225566/1‑93 (MQ=255)

               ccccAATTAAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCa                                                                          <  1:5649149/93‑1 (MQ=255)

                cccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAg                                                                         >  1:5311404/1‑93 (MQ=255)

                cccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAg                                                                         >  1:6636422/1‑93 (MQ=255)

                cccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAg                                                                         >  1:7179564/1‑93 (MQ=255)

                 ccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGt                                                                        >  2:2724000/1‑93 (MQ=255)

                    aTTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTtta                                                                     <  1:3195669/93‑1 (MQ=255)

                    aTTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTtta                                                                     <  2:2779559/93‑1 (MQ=255)

                    aTTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTtta                                                                     >  1:2157573/1‑93 (MQ=255)

                    aTTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAAC‑ACAGACGTCAGTttat                                                                    <  2:3906180/93‑1 (MQ=255)

                      tCAAATATCATGCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTttatt                                                                   >  2:4653147/1‑93 (MQ=255)

                      tCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTttatt                                                                   <  1:2947419/93‑1 (MQ=255)

                       cAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTttatta                                                                  <  1:6161381/93‑1 (MQ=255)

                       cAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTttatta                                                                  <  1:5939898/93‑1 (MQ=255)

                        aaaTATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTAc                                                                 <  2:7455237/93‑1 (MQ=255)

                         aaTATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACt                                                                >  2:5806091/1‑93 (MQ=255)

                         aaTATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACt                                                                >  2:3792538/1‑93 (MQ=255)

                           tatCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCg                                                              >  2:4720511/1‑93 (MQ=255)

                           tatCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCg                                                              <  2:6746753/93‑1 (MQ=255)

                           tatCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCg                                                              >  2:6873188/1‑93 (MQ=255)

                            atcatcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCgg                                                             <  1:6496219/93‑1 (MQ=255)

                            atcatcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCgg                                                             <  1:447937/93‑1 (MQ=255)

                            atcatcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCgg                                                             <  2:5927054/93‑1 (MQ=255)

                               atcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTc                                                          >  2:2601990/1‑93 (MQ=255)

                                ccccAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTcc                                                         <  1:2791047/92‑1 (MQ=255)

                                   cAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCtt                                                      >  2:4427992/1‑93 (MQ=255)

                                    aaGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCttt                                                     >  2:5508290/1‑93 (MQ=255)

                                    aaGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCttt                                                     <  2:4532177/93‑1 (MQ=255)

                                    aaGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCttt                                                     >  2:2235083/1‑93 (MQ=255)

                                      gCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTa                                                   >  2:4506244/1‑93 (MQ=255)

                                      gCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTa                                                   >  2:125893/1‑93 (MQ=255)

                                       ccACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAt                                                  >  1:3499707/1‑93 (MQ=255)

                                       ccACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAt                                                  <  2:1751985/93‑1 (MQ=255)

                                       ccACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAt                                                  <  2:4927900/93‑1 (MQ=255)

                                        cACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAtt                                                 <  1:1127483/93‑1 (MQ=255)

                                           tGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCg                                              >  2:4004057/1‑93 (MQ=255)

                                           tGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCg                                              >  2:1376449/1‑93 (MQ=255)

                                            ggAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGt                                             >  1:7365912/1‑93 (MQ=255)

                                             gAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGtt                                            <  2:1810310/93‑1 (MQ=255)

                                               aaCAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTtga                                          >  1:2697705/1‑93 (MQ=255)

                                                 cAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTtgatg                                        <  1:5855217/93‑1 (MQ=255)

                                                  aGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTtgatga                                       >  2:5404935/1‑93 (MQ=255)

                                                     cgcATTTGAACAGGATAAAACCGTAGAGTGCGATTGGCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGt                                    <  1:3423019/93‑1 (MQ=255)

                                                      gcATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGtt                                   >  2:2552003/1‑93 (MQ=255)

                                                      gcATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGtt                                   <  1:1990976/93‑1 (MQ=255)

                                                       cATTTGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTg                                  >  2:6644491/1‑93 (MQ=255)

                                                       cATTTGAACAGGATAAAACCGTAGAGTGCAATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTg                                  >  1:6040701/1‑93 (MQ=255)

                                                         tttGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGaa                                >  1:67117/1‑93 (MQ=255)

                                                         tttGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGaa                                <  1:2267292/93‑1 (MQ=255)

                                                         tttGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGaa                                >  1:937107/1‑93 (MQ=255)

                                                         tttGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGaa                                >  2:3747516/1‑93 (MQ=255)

                                                          ttGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAAc                               <  1:3754603/93‑1 (MQ=255)

                                                          ttGAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAAc                               <  2:5331096/93‑1 (MQ=255)

                                                            gAACAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTCCGTTGATGAAGTTGAACAt                             >  1:419827/1‑93 (MQ=255)

                                                             aaCAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACAtc                            >  1:7008115/1‑93 (MQ=255)

                                                              aCAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACAtct                           >  2:896504/1‑93 (MQ=255)

                                                               cAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACAtctc                          <  2:3465810/93‑1 (MQ=255)

                                                               cAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACAtctc                          <  2:6289047/93‑1 (MQ=255)

                                                               cAGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACAtctc                          <  1:5769645/93‑1 (MQ=255)

                                                                aGGATAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATctct                         <  1:5115443/93‑1 (MQ=255)

                                                                   aTAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTc                      <  1:3621028/93‑1 (MQ=255)

                                                                    tAAAACCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTcc                     <  1:978467/93‑1 (MQ=255)

                                                                     aaaaCCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCg                    <  2:6789942/93‑1 (MQ=255)

                                                                      aaaCCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCgt                   >  1:4660743/1‑93 (MQ=255)

                                                                      aaaCCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCgt                   >  2:1571889/1‑93 (MQ=255)

                                                                       aaCCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCgtg                  <  1:4284629/93‑1 (MQ=255)

                                                                        aCCGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCgtgg                 <  2:3398178/93‑1 (MQ=255)

                                                                         ccGTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCgtggt                <  1:5646000/93‑1 (MQ=255)

                                                                          cGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCgtggtg               >  2:5203101/1‑93 (MQ=255)

                                                                           gTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGgtgt              >  2:5991134/1‑93 (MQ=255)

                                                                           gTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGgtgt              >  1:5975665/1‑93 (MQ=255)

                                                                           gTAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGgtgt              <  1:2212895/93‑1 (MQ=255)

                                                                            tAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTa             >  2:7089976/1‑93 (MQ=255)

                                                                            tAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTa             <  2:417418/93‑1 (MQ=255)

                                                                            tAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTa             <  2:2051405/93‑1 (MQ=255)

                                                                            tAGAGTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTa             <  1:4562316/93‑1 (MQ=255)

                                                                             agagTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTAt            <  1:17708/93‑1 (MQ=255)

                                                                              gagTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTAtt           >  2:6123308/1‑93 (MQ=255)

                                                                              gagTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCC‑TGGTGGATTg          >  2:4981643/1‑93 (MQ=255)

                                                                                gTGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGc         >  2:1087047/1‑93 (MQ=255)

                                                                                 tGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCg        >  1:5037448/1‑93 (MQ=255)

                                                                                 tGCGATTGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCg        >  2:5611164/1‑93 (MQ=255)

                                                                                      ttGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACgg   <  1:3679681/93‑1 (MQ=255)

                                                                                       tGCCGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGt  <  2:1334800/93‑1 (MQ=255)

AATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGT  >  NC_000913/3895660‑3895839

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