BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	409,218	A→G	V9A (GTT→GCT)	yaiD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	409,218	0	A	G	99.4%	540.9 / ‑5.0	165	V9A (GTT→GCT)	yaiD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (74/90); total (75/90)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 173 total aligned reads displayed.

TGTCCTGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAACCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCTTGCA  >  NC_000913/409126‑409308

  tCCTGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGcc                                                                                          >  1:1280962/1‑93 (MQ=255)

    cTGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCAt                                                                                        <  1:157196/93‑1 (MQ=255)

     tGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCAtt                                                                                       <  2:3273832/93‑1 (MQ=255)

     tGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCAtt                                                                                       >  1:28699/1‑93 (MQ=255)

      ggCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTa                                                                                      >  1:3960322/1‑93 (MQ=255)

       gCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTaa                                                                                     >  2:2864697/1‑93 (MQ=255)

         tGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAt                                                                                   <  2:4018185/93‑1 (MQ=255)

            cGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAtttt                                                                                <  2:2131100/93‑1 (MQ=255)

              cATGGGGTAAATGCCATCGCGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCACGGCTAAGAAGGTCAGCCATTAAATTTTTg                                                                              >  2:4524858/1‑93 (MQ=255)

               aTGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGa                                                                             >  2:4463922/1‑93 (MQ=255)

               aTGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGa                                                                             <  1:1582535/93‑1 (MQ=255)

               aTGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGa                                                                             >  2:3136889/1‑93 (MQ=255)

                tGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGaa                                                                            >  1:444906/1‑93 (MQ=255)

                tGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGaa                                                                            <  1:2249941/93‑1 (MQ=255)

                 ggggTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAAc                                                                           <  2:3243083/93‑1 (MQ=255)

                  gggTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAAcc                                                                          <  2:1957084/93‑1 (MQ=255)

                  gggTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAAcc                                                                          >  2:3215640/1‑93 (MQ=255)

                    gTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACcac                                                                        <  1:3720154/93‑1 (MQ=255)

                     tAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACcaca                                                                       >  2:18233/1‑93 (MQ=255)

                     tAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACcaca                                                                       <  2:2381647/93‑1 (MQ=255)

                     tAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACcaca                                                                       <  2:1305333/93‑1 (MQ=255)

                     tAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACcaca                                                                       <  1:105929/93‑1 (MQ=255)

                      aaaTGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAg                                                                      >  2:135891/1‑93 (MQ=255)

                       aaTGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGc                                                                     >  1:571822/1‑93 (MQ=255)

                         tGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCAt                                                                   >  1:1197663/1‑93 (MQ=255)

                           ccATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCATCGGAATCTCGCGGCTAAGACGGTCAGCCATTAAATTTTTGAACCACAGCATGa                                                                 >  2:2597569/1‑93 (MQ=255)

                           ccATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGa                                                                 <  1:1360121/93‑1 (MQ=255)

                            cATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGAt                                                                >  1:367900/1‑93 (MQ=255)

                              tCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATaa                                                              >  1:2685404/1‑93 (MQ=255)

                               cGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAAt                                                             <  2:3219834/93‑1 (MQ=255)

                               cGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAAt                                                             <  1:4651571/93‑1 (MQ=255)

                                gAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAAtt                                                            >  2:436317/1‑93 (MQ=255)

                                gAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAAtt                                                            <  1:3218486/93‑1 (MQ=255)

                                gAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAAtt                                                            >  2:1079873/1‑93 (MQ=255)

                                   gCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTcc                                                         >  1:562307/1‑93 (MQ=255)

                                    ccaGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCa                                                        >  1:2729699/3‑93 (MQ=255)

                                    cTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCa                                                        <  2:656101/93‑1 (MQ=255)

                                    cTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCa                                                        <  1:1491510/93‑1 (MQ=255)

                                     tAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCAc                                                       <  2:2349332/93‑1 (MQ=255)

                                     tAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCAc                                                       >  1:1030083/1‑93 (MQ=255)

                                     tAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCAc                                                       >  1:1019980/1‑93 (MQ=255)

                                     tAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCAc                                                       <  1:1764895/93‑1 (MQ=255)

                                     tAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCAc                                                       <  1:3469246/93‑1 (MQ=255)

                                      aGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACg                                                      >  2:4733686/1‑93 (MQ=255)

                                      aGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACg                                                      >  2:3584066/1‑93 (MQ=255)

                                      aGCTGTTTCTCCATCTCTTCTGCACGCAGCGAAATCTCGGGGCTACGACGTTAAGCCATTAAATTTGTGAACCACAGCATGAGAATTTCCACg                                                      >  1:153166/1‑93 (MQ=255)

                                       gCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACgg                                                     <  2:3881188/93‑1 (MQ=255)

                                          ttttttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCt                                                  <  2:1956323/92‑1 (MQ=255)

                                          gTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCt                                                  <  2:2501077/93‑1 (MQ=255)

                                           tttttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCtt                                                 <  2:468423/93‑1 (MQ=255)

                                            ttttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTg                                                <  1:279713/93‑1 (MQ=255)

                                             tttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGt                                               <  2:1799638/93‑1 (MQ=255)

                                             tttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGt                                               <  2:3028485/93‑1 (MQ=255)

                                             tttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGt                                               <  1:4413141/93‑1 (MQ=255)

                                             tttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGt                                               <  1:476466/93‑1 (MQ=255)

                                              ttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTc                                              <  2:1859818/93‑1 (MQ=255)

                                               tCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCg                                             <  1:2914727/93‑1 (MQ=255)

                                               tCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCg                                             >  2:518718/1‑93 (MQ=255)

                                                ccATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGt                                            <  1:1436415/93‑1 (MQ=255)

                                                   tctcTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTaa                                         >  2:2537654/1‑93 (MQ=255)

                                                   tctcTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTaa                                         <  1:2117150/93‑1 (MQ=255)

                                                     tcttctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGTTAATTTCCACGGCCTTGTCGTTAAAt                                       >  2:2230788/1‑93 (MQ=255)

                                                     tcttctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACCGCCTTGTCGTTAAAt                                       >  2:3120132/1‑93 (MQ=255)

                                                     tcttctGCACGCAGCGAAATCTCGCGGCTAAGAAGGTAAGCTATTAAATTTTTGAACCACAGAATGATAATTTCCACGGTCTTGTCGTTAAAt                                       >  2:2965645/1‑93 (MQ=255)

                                                       ttctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAAttt                                     >  1:357256/1‑93 (MQ=255)

                                                       ttctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAAttt                                     <  2:4540897/93‑1 (MQ=255)

                                                        tctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTa                                    <  2:3099723/93‑1 (MQ=255)

                                                          tGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGc                                  >  2:4157987/1‑93 (MQ=255)

                                                          tGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGc                                  >  1:4513869/1‑93 (MQ=255)

                                                          tGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCAATTAAATTTTTGAACCACAGCATGATAATTTACACGGCCTTGTCGTTAAATTTAGc                                  >  2:4636652/1‑93 (MQ=255)

                                                           gCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCg                                 <  1:403915/93‑1 (MQ=255)

                                                               gCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCa                             <  2:1941041/93‑1 (MQ=255)

                                                               gCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCa                             >  1:1420065/1‑93 (MQ=255)

                                                               gCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCa                             <  1:765584/93‑1 (MQ=255)

                                                                   cGAAATCTCGCGGCTAAGACGGTTAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGAt                         <  2:3387705/93‑1 (MQ=255)

                                                                   cGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGAt                         <  1:1387751/93‑1 (MQ=255)

                                                                    gAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATa                        <  2:4653443/93‑1 (MQ=255)

                                                                    gAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATa                        >  1:2941338/1‑93 (MQ=255)

                                                                      aaTCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAAc                      >  1:1936703/1‑93 (MQ=255)

                                                                       aTCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACg                     >  1:1626149/1‑93 (MQ=255)

                                                                       aTCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACg                     >  1:3432198/1‑93 (MQ=255)

                                                                        tctcGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGa                    <  1:2661270/93‑1 (MQ=255)

                                                                         ctcGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGaa                   <  2:1566061/93‑1 (MQ=255)

                                                                         ctcGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGaa                   <  1:4605524/93‑1 (MQ=255)

                                                                          tcGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTGCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAAt                  >  2:3503147/1‑93 (MQ=255)

                                                                          tcGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAAt                  >  2:4016060/1‑93 (MQ=255)

                                                                          tcGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAAt                  >  1:910945/1‑93 (MQ=255)

                                                                           cgcgGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAAtt                 <  1:1678496/93‑1 (MQ=255)

                                                                           cgcgGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAAtt                 <  1:2828796/93‑1 (MQ=255)

                                                                           cgcgGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAAtt                 <  2:551337/93‑1 (MQ=255)

                                                                             cgGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGt               <  1:199331/93‑1 (MQ=255)

                                                                              ggCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTc              >  1:3907053/1‑93 (MQ=255)

                                                                               gCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTcg             <  1:3253628/93‑1 (MQ=255)

                                                                                cTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTcgg            >  2:801471/1‑93 (MQ=255)

                                                                                  aaGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTcggcg          >  1:1909576/1‑93 (MQ=255)

                                                                                   aGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTcgccgg         >  1:3563086/1‑93 (MQ=255)

                                                                                     aCGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGcc       <  1:2864697/93‑1 (MQ=255)

                                                                                      cGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCt      <  1:3794716/93‑1 (MQ=255)

                                                                                       ggTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCtt     <  1:4444831/93‑1 (MQ=255)

                                                                                          aaaCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCTTGCa  <  1:1976033/93‑1 (MQ=255)

TGTCCTGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAACCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCTTGCA  >  NC_000913/409126‑409308

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