BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,556,716	C→T	G127D (GGC→GAC)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,556,716	0	C	T	85.4%	867.3 / 93.5	315	G127D (GGC→GAC)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand): ref base C (19/27); new base T (148/121); total (167/148)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

Only 100 of 320 total aligned reads displayed.

ACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGAC  >  NC_000913/4556623‑4556807

acacGATCGATAATCGCCACGTCTTTTTCC‑CGGAACCCGGAATGGGGCGGGGAGGGACCTGATAAGCGCCGGGCCACAATCAGGGGCCGATGc                                                                                             >  1:5422851/1‑93 (MQ=255)

   cgatcgatAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTct                                                                                           >  1:3410593/1‑93 (MQ=255)

      tcgatAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGtcttct                                                                                        >  2:5615645/1‑93 (MQ=255)

      tcgatAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGtcttct                                                                                        >  1:6468536/1‑93 (MQ=255)

         atAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCa                                                                                     >  1:7388950/1‑93 (MQ=255)

          tAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCAt                                                                                    >  2:899619/1‑93 (MQ=255)

          tAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCAt                                                                                    <  2:4430391/93‑1 (MQ=255)

            aTCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTg                                                                                  >  2:1615065/1‑93 (MQ=255)

               gCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAgc                                                                               >  1:4059653/1‑93 (MQ=255)

                 cACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAgcgc                                                                             >  1:7228692/1‑93 (MQ=255)

                  aCGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCa                                                                            >  2:324760/1‑93 (MQ=255)

                  aCGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCa                                                                            <  1:2329568/93‑1 (MQ=255)

                   cGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCAc                                                                           <  2:3177640/93‑1 (MQ=255)

                    gTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACg                                                                          >  2:5459701/1‑93 (MQ=255)

                    gTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACg                                                                          <  1:4839816/93‑1 (MQ=255)

                     tCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACgg                                                                         <  2:4705405/93‑1 (MQ=255)

                      cTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACggg                                                                        >  2:1476949/1‑93 (MQ=255)

                       ttttcCCCCGGACCCCG‑AATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTc                                                                      <  1:7114406/93‑1 (MQ=255)

                         tttCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCt                                                                     <  1:2370019/93‑1 (MQ=255)

                         tttCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCt                                                                     <  2:102331/93‑1 (MQ=255)

                          ttCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCtt                                                                    <  2:4578082/93‑1 (MQ=255)

                           tCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTg                                                                   >  1:4410255/1‑93 (MQ=255)

                           tCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTg                                                                   <  1:1846838/93‑1 (MQ=255)

                             cACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGc                                                                 >  2:2005925/1‑93 (MQ=255)

                              aCGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCg                                                                <  2:4154957/93‑1 (MQ=255)

                                  aaCCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCa                                                            <  1:6051399/93‑1 (MQ=255)

                                  aaCCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCa                                                            <  2:5684901/93‑1 (MQ=255)

                                  aaCCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCa                                                            <  1:1091662/93‑1 (MQ=255)

                                    cccGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAgg                                                          <  2:5096860/93‑1 (MQ=255)

                                      cGTAATGGCGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGa                                                        >  1:6512498/1‑93 (MQ=255)

                                        tAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGAtt                                                      >  1:4330842/1‑93 (MQ=255)

                                        tAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGAtt                                                      >  2:7510855/1‑93 (MQ=255)

                                        tAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGAtt                                                      >  2:718803/1‑93 (MQ=255)

                                        tAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGAtt                                                      >  2:3064762/1‑93 (MQ=255)

                                        tAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGAtt                                                      <  2:627638/93‑1 (MQ=255)

                                        tAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGAtt                                                      >  2:8033/1‑93 (MQ=255)

                                          aTGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTcc                                                    >  1:2990176/1‑93 (MQ=255)

                                           tGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCg                                                   <  1:3937686/93‑1 (MQ=255)

                                            ggTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCgg                                                  <  2:6960766/93‑1 (MQ=255)

                                             gTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCggg                                                 <  1:588273/93‑1 (MQ=255)

                                              tGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGt                                                <  1:2917269/93‑1 (MQ=255)

                                              tGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGt                                                <  2:3414875/93‑1 (MQ=255)

                                               gCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTg                                               <  1:6305528/93‑1 (MQ=255)

                                                 gggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGc                                             >  1:1010215/1‑93 (MQ=255)

                                                  ggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCg                                            >  1:2581902/1‑93 (MQ=255)

                                                   gAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCga                                           <  1:3759526/93‑1 (MQ=255)

                                                   gAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCga                                           >  1:7481420/1‑93 (MQ=255)

                                                   gAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCga                                           <  1:1942833/93‑1 (MQ=255)

                                                     aGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCgaga                                         >  1:2099487/1‑93 (MQ=255)

                                                     aGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCgaga                                         >  2:4608958/1‑93 (MQ=255)

                                                      gggACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCgagag                                        >  1:3061961/1‑93 (MQ=255)

                                                       ggACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCgagaga                                       <  2:1143091/93‑1 (MQ=255)

                                                       ggACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCgagaga                                       >  2:4850919/1‑93 (MQ=255)

                                                        gACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGAt                                      <  1:2153286/93‑1 (MQ=255)

                                                        gACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGAt                                      <  1:5043337/93‑1 (MQ=255)

                                                           aTGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATaga                                   <  2:4985318/93‑1 (MQ=255)

                                                           aTGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATaga                                   <  1:7013368/93‑1 (MQ=255)

                                                           aTGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATaga                                   >  2:4364847/1‑93 (MQ=255)

                                                           aTGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATaga                                   <  2:5735064/93‑1 (MQ=255)

                                                            tGATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATagag                                  >  2:586929/1‑93 (MQ=255)

                                                             gATAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGt                                 <  1:6909022/93‑1 (MQ=255)

                                                              aTAAGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTc                                <  1:1698259/93‑1 (MQ=255)

                                                              aTAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTc                                <  1:310857/93‑1 (MQ=255)

                                                               taggcgcCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCg                               <  2:287552/90‑1 (MQ=255)

                                                                aaGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCgg                              <  2:1342185/93‑1 (MQ=255)

                                                                 aGCGCCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGt                             >  2:2561231/1‑93 (MQ=255)

                                                                  gcgcCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTg                            <  2:227677/93‑1 (MQ=255)

                                                                   cgcCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGc                           <  2:2139331/93‑1 (MQ=255)

                                                                    gcCGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGcc                          >  2:2378801/1‑93 (MQ=255)

                                                                     ccGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGccc                         <  2:3636250/93‑1 (MQ=255)

                                                                     ccGGTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGccc                         <  1:5503959/93‑1 (MQ=255)

                                                                        gTCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCcagc                      >  2:4278868/1‑93 (MQ=255)

                                                                         tCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCcagca                     <  2:2641319/93‑1 (MQ=255)

                                                                         tCAGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCcagca                     >  2:6893666/1‑93 (MQ=255)

                                                                         tCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCcagca                     >  2:5990636/1‑93 (MQ=255)

                                                                           aGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCTCAGCAGa                   >  2:7031902/1‑93 (MQ=255)

                                                                           aGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGa                   >  1:4768897/1‑93 (MQ=255)

                                                                           aGCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGa                   >  2:4744970/1‑93 (MQ=255)

                                                                            gCATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGac                  >  2:1136103/1‑93 (MQ=255)

                                                                            gCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGac                  >  1:4264767/1‑93 (MQ=255)

                                                                             cATCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGacc                 <  1:4655650/93‑1 (MQ=255)

                                                                              aTCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGacca                >  1:6825651/1‑93 (MQ=255)

                                                                              aTCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGacca                >  1:227094/1‑93 (MQ=255)

                                                                              aTCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGacca                >  1:620193/1‑93 (MQ=255)

                                                                              aTCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGacca                >  2:67793/1‑93 (MQ=255)

                                                                               tCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGaccaa               >  1:2704700/1‑93 (MQ=255)

                                                                               tCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGaccaa               >  2:6349069/1‑93 (MQ=255)

                                                                               tCCAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGaccaa               <  1:326158/93‑1 (MQ=255)

                                                                               tCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGaccaa               <  1:4175748/93‑1 (MQ=255)

                                                                                ccAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGaccaac              >  1:2001349/1‑93 (MQ=255)

                                                                                ccAGGCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGaccaac              <  1:5179881/93‑1 (MQ=255)

                                                                                   ggCGCTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCAc           <  2:3777293/93‑1 (MQ=255)

                                                                                      gcTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtga        >  2:6560171/1‑93 (MQ=255)

                                                                                      gcTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtga        <  2:1739882/93‑1 (MQ=255)

                                                                                      gcTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtga        <  1:2469240/93‑1 (MQ=255)

                                                                                       cTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgac       >  1:3520030/1‑93 (MQ=255)

                                                                                       cTGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGGGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgac       >  1:1630420/1‑93 (MQ=255)

                                                                                        tGATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgacg      >  1:5577262/1‑93 (MQ=255)

                                                                                         gATGTCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgacgt     <  2:6667248/93‑1 (MQ=255)

                                                                                         gATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgacgt     <  1:1244047/93‑1 (MQ=255)

ACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGAC  >  NC_000913/4556623‑4556807

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