BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,866,583	A→G	S119S (AGT→AGC)	yidR ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,866,583	0	A	G	84.1%	729.0 / 125.0	320	S119S (AGT→AGC)	yidR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (32/19); new base G (128/141); total (160/160)
Fisher's exact test for biased strand distribution p-value = 6.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

Only 100 of 322 total aligned reads displayed.

CTTTCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAA  >  NC_000913/3866491‑3866674

cTTTCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTa                                                                                             >  1:4046983/1‑93 (MQ=255)

   tCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTc                                                                                          <  2:4459814/93‑1 (MQ=255)

     aCCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCAt                                                                                        >  2:1992361/1‑93 (MQ=255)

       cGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCt                                                                                      <  1:4299165/93‑1 (MQ=255)

       cGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCt                                                                                      <  1:3788877/93‑1 (MQ=255)

       cGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCt                                                                                      <  2:5249801/93‑1 (MQ=255)

       cGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCt                                                                                      <  2:6935852/93‑1 (MQ=255)

       cGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCt                                                                                      <  1:1007853/93‑1 (MQ=255)

        gTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCtt                                                                                     <  2:3815239/93‑1 (MQ=255)

        gTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCtt                                                                                     >  1:2938002/1‑93 (MQ=255)

        gTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCtt                                                                                     >  1:5888106/1‑93 (MQ=255)

         ttCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTg                                                                                    >  2:4715722/1‑93 (MQ=255)

          tCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTgc                                                                                   <  1:7113158/93‑1 (MQ=255)

          tCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTgc                                                                                   <  2:5897378/93‑1 (MQ=255)

             ggCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTgccgc                                                                                >  2:13190/1‑93 (MQ=255)

               cTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCt                                                                              >  1:5115325/1‑93 (MQ=255)

                  aaGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTca                                                                           <  1:134291/93‑1 (MQ=255)

                   aGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTcaa                                                                          >  2:336103/1‑93 (MQ=255)

                    gACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTcaac                                                                         <  1:6256473/93‑1 (MQ=255)

                     aCATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTcaaca                                                                        <  2:38803/93‑1 (MQ=255)

                      cATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTcaacaa                                                                       <  1:3449484/93‑1 (MQ=255)

                      cATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTGCCGCCTTcaacaa                                                                       <  1:4332053/93‑1 (MQ=255)

                       aTGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAAt                                                                      >  2:7524882/1‑93 (MQ=255)

                        tGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATc                                                                     <  2:6279771/93‑1 (MQ=255)

                             aTGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTcc                                                                <  1:5992605/93‑1 (MQ=255)

                               ggCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCcgc                                                              <  2:4635074/93‑1 (MQ=255)

                                 cTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCgcgc                                                            <  2:1091306/93‑1 (MQ=255)

                                     cgccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgat                                                        <  1:5529703/93‑1 (MQ=255)

                                      gccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgatg                                                       >  1:6367046/1‑93 (MQ=255)

                                      gccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgatg                                                       >  2:2786351/1‑93 (MQ=255)

                                       ccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgatga                                                      <  2:5082619/93‑1 (MQ=255)

                                       ccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgatga                                                      <  1:5179878/93‑1 (MQ=255)

                                       ccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgatga                                                      <  2:6184690/93‑1 (MQ=255)

                                       ccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgatga                                                      <  2:2322501/93‑1 (MQ=255)

                                         gcgcAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGatgatg                                                    >  2:2100187/1‑93 (MQ=255)

                                            cagcaCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGaa                                                 >  1:4141145/1‑93 (MQ=255)

                                             agcaCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAg                                                <  1:1642536/93‑1 (MQ=255)

                                             agcaCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAg                                                <  1:5734320/93‑1 (MQ=255)

                                              gcaCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGt                                               <  1:6516255/93‑1 (MQ=255)

                                               caCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTc                                              <  2:498737/93‑1 (MQ=255)

                                                aCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCg                                             >  2:6588864/1‑93 (MQ=255)

                                                 cTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGt                                            >  1:3572292/1‑93 (MQ=255)

                                                  tCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTa                                           <  2:3922705/93‑1 (MQ=255)

                                                   ccTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTaa                                          <  1:2239007/93‑1 (MQ=255)

                                                     tGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATg                                        >  1:4168793/1‑93 (MQ=255)

                                                      ggCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGc                                       >  1:250872/1‑93 (MQ=255)

                                                       gCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGcc                                      >  1:6123687/1‑93 (MQ=255)

                                                         gTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCAt                                    <  2:4108128/93‑1 (MQ=255)

                                                         gTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCAt                                    >  2:5538927/1‑93 (MQ=255)

                                                         gTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCAt                                    <  1:2835309/93‑1 (MQ=255)

                                                         gTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCAt                                    >  1:3073890/1‑93 (MQ=255)

                                                          tataCGGCGCGGTAATATCCATCGCATCAAGGTTGTTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATg                                   >  1:4803376/1‑93 (MQ=255)

                                                          tataCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATg                                   >  1:2190543/1‑93 (MQ=255)

                                                          tataCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATg                                   >  2:1687013/1‑93 (MQ=255)

                                                            taCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGtt                                 >  2:3539340/1‑93 (MQ=255)

                                                            taCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGtt                                 <  1:4446386/93‑1 (MQ=255)

                                                             aCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCTTCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGttt                                >  2:6778632/1‑93 (MQ=255)

                                                             aCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGttt                                <  2:5718163/93‑1 (MQ=255)

                                                              cGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCTTAATGCCATGTTtc                               <  2:5183893/93‑1 (MQ=255)

                                                              cGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTtc                               <  2:1519259/93‑1 (MQ=255)

                                                              cGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTtc                               <  2:7233479/93‑1 (MQ=255)

                                                                gcgcGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTtcat                             <  2:3659504/93‑1 (MQ=255)

                                                                gcgcGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTtcat                             >  2:3121559/1‑93 (MQ=255)

                                                                 cgcgGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTtcatc                            <  1:4739305/93‑1 (MQ=255)

                                                                  gcgGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTtcatca                           <  1:1277176/93‑1 (MQ=255)

                                                                   cgGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAg                          <  2:104543/93‑1 (MQ=255)

                                                                     gTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGa                        <  1:56943/93‑1 (MQ=255)

                                                                     gTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGa                        <  1:2034617/93‑1 (MQ=255)

                                                                      tAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAt                       <  2:181536/93‑1 (MQ=255)

                                                                      tAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAt                       >  2:4292656/1‑93 (MQ=255)

                                                                       aaTATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAtt                      >  2:4863995/1‑93 (MQ=255)

                                                                       aaTATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAtt                      >  2:5146720/1‑93 (MQ=255)

                                                                        atatCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAttt                     <  1:2654157/93‑1 (MQ=255)

                                                                        atatCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAttt                     >  1:382198/1‑93 (MQ=255)

                                                                         tatCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAtttt                    >  2:1130961/1‑93 (MQ=255)

                                                                          atCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTc                   >  2:4303102/1‑93 (MQ=255)

                                                                           tCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTcc                  <  1:813687/93‑1 (MQ=255)

                                                                           tCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTcc                  <  1:7297138/93‑1 (MQ=255)

                                                                             cATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCgg                <  1:3034644/93‑1 (MQ=255)

                                                                               tCGCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGc              >  1:6799868/1‑93 (MQ=255)

                                                                                 gCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCg            >  1:4582198/1‑93 (MQ=255)

                                                                                 gCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCg            >  1:6589718/1‑93 (MQ=255)

                                                                                 gCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCg            >  2:3423526/1‑93 (MQ=255)

                                                                                  cATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGt           <  2:851653/93‑1 (MQ=255)

                                                                                  cATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGt           >  2:3236298/1‑93 (MQ=255)

                                                                                  cATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGt           <  1:1514421/93‑1 (MQ=255)

                                                                                   aTCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTg          <  1:3280119/93‑1 (MQ=255)

                                                                                    tCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGtnn         >  2:1896725/1‑91 (MQ=255)

                                                                                    tCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTgg         >  1:3941964/1‑93 (MQ=255)

                                                                                    tCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGTCATGTTTCATCAGGATTTTCCGGGCCGTgg         >  2:4358437/1‑93 (MQ=255)

                                                                                     cAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGa        <  2:648387/93‑1 (MQ=255)

                                                                                      aaGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGAt       <  2:4230243/93‑1 (MQ=255)

                                                                                      aaGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGAt       <  2:2908273/93‑1 (MQ=255)

                                                                                      aaGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGAt       <  2:2274448/93‑1 (MQ=255)

                                                                                       aGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATg      >  1:1932449/1‑93 (MQ=255)

                                                                                        ggTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGa     <  1:1166066/93‑1 (MQ=255)

                                                                                        ggTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGa     >  2:7495841/1‑93 (MQ=255)

                                                                                         gTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGaa    >  1:6871873/1‑93 (MQ=255)

                                                                                         gTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGaa    <  1:5490251/93‑1 (MQ=255)

                                                                                           tACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGaaaa  >  1:4950389/1‑93 (MQ=255)

CTTTCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAA  >  NC_000913/3866491‑3866674

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