BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,019,825	G→A	P163P (CCG→CCA)	b2881 →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,019,825	0	G	A	99.0%	673.3 / ‑3.7	194	P163P (CCG→CCA)	b2881	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (70/122); total (70/124)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

Only 100 of 200 total aligned reads displayed.

CCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGA  >  NC_000913/3019735‑3019914

ccGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAAc                                                                                         <  2:3805927/93‑1 (MQ=255)

 cgcgATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACt                                                                                        >  2:4027657/1‑93 (MQ=255)

 cgccATGCTGGCCGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTTCACGTAAACATAACCCGCAGGCCACCATTGATATCGCTCCAACt                                                                                        <  1:805717/93‑1 (MQ=255)

    gATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTc                                                                                     <  2:3107288/93‑1 (MQ=255)

      tGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCg                                                                                   >  2:1854789/1‑93 (MQ=255)

      tGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCg                                                                                   <  1:4406710/93‑1 (MQ=255)

      tGCTGGCCGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCCCCATTGCTATCGCTCCAACTTTCCg                                                                                   >  2:3586641/1‑93 (MQ=255)

         tggctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTga                                                                                <  1:4451862/93‑1 (MQ=255)

         tggctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTga                                                                                >  2:4178705/1‑93 (MQ=255)

         tggctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTga                                                                                <  2:3785449/93‑1 (MQ=255)

          ggctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgac                                                                               >  1:3611830/1‑93 (MQ=255)

           gctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacg                                                                              >  1:2490792/1‑93 (MQ=255)

           gctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacg                                                                              >  2:36622/1‑93 (MQ=255)

           gctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacg                                                                              <  2:808997/93‑1 (MQ=255)

            ctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacga                                                                             <  2:1225982/93‑1 (MQ=255)

            ctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacga                                                                             <  1:306858/93‑1 (MQ=255)

              ggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGAcc                                                                           <  2:4139523/93‑1 (MQ=255)

              ggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGAcc                                                                           <  2:1297638/93‑1 (MQ=255)

               gcagcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCt                                                                          <  1:131700/93‑1 (MQ=255)

                cagcaATACTATCAGGTCATTGAACTGTCGGTTGCACGTAAAAATAAGCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTa                                                                         <  2:3002940/93‑1 (MQ=255)

                 agcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTaa                                                                        >  1:3158494/1‑92 (MQ=255)

                 agcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAg                                                                        <  2:2761832/93‑1 (MQ=255)

                 agcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAg                                                                        >  2:1129856/1‑93 (MQ=255)

                   caATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGaa                                                                      <  1:1387320/93‑1 (MQ=255)

                     aTACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGt                                                                    <  2:80359/93‑1 (MQ=255)

                     aTACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGt                                                                    >  1:4137186/1‑93 (MQ=255)

                       aCTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCa                                                                  <  1:284303/93‑1 (MQ=255)

                       aCTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTGCCGTGACGACCTAGAAGTCa                                                                  <  1:1951665/93‑1 (MQ=255)

                        cTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAt                                                                 <  1:3868587/93‑1 (MQ=255)

                         tATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCAtt                                                                <  1:2084570/93‑1 (MQ=255)

                         tATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAtt                                                                <  2:2268112/93‑1 (MQ=255)

                         tATCAGGTCATTGAACTGGAGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAtt                                                                <  2:2521306/93‑1 (MQ=255)

                          aTCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTg                                                               >  2:955190/1‑93 (MQ=255)

                            cAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGt                                                             <  2:712287/93‑1 (MQ=255)

                             aGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGAAGACCTAGAATTCATTggga                                                            >  2:2469234/1‑91 (MQ=255)

                              ggTCATTGAAGTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTaa                                                           <  1:1911697/93‑1 (MQ=255)

                              ggTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTaa                                                           <  2:1492489/93‑1 (MQ=255)

                              ggTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGAACTAGAAGTCATTGGTaa                                                           <  1:3514686/93‑1 (MQ=255)

                               gTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAg                                                          <  1:4031980/93‑1 (MQ=255)

                                tCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGc                                                         <  2:1986026/93‑1 (MQ=255)

                                 cATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCa                                                        >  2:1366182/1‑93 (MQ=255)

                                  aTTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAt                                                       >  1:1623733/1‑93 (MQ=255)

                                   ttGAGCTGGCGGTTGCCCGTAAAAATAATCCGCAGGCCACCATTGATCTCGCTCCAACTTTCCGTGACGACCTAGAAGTCATAAGTCAGCAtt                                                      <  2:4225100/93‑1 (MQ=255)

                                    tGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTa                                                     <  2:1928166/93‑1 (MQ=255)

                                     gAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTAt                                                    <  2:2077544/93‑1 (MQ=255)

                                     gAACTGGCGGTGGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTAt                                                    <  1:724231/93‑1 (MQ=255)

                                     gAAATGGCGGTTGTACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTAt                                                    <  1:369625/93‑1 (MQ=255)

                                       aCTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATcc                                                  >  1:4109738/1‑93 (MQ=255)

                                          ggCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTaa                                               >  1:1650218/1‑93 (MQ=255)

                                          ggCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTaa                                               <  2:3494286/93‑1 (MQ=255)

                                          ggCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTaa                                               >  1:3843943/1‑93 (MQ=255)

                                            cGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTaaaa                                             <  2:1108366/93‑1 (MQ=255)

                                             ggTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAAc                                            <  1:269430/93‑1 (MQ=255)

                                              gTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACt                                           <  2:271257/93‑1 (MQ=255)

                                                tGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGa                                         <  1:4288223/93‑1 (MQ=255)

                                                tGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGa                                         <  2:3222177/93‑1 (MQ=255)

                                                tGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGa                                         <  1:544294/93‑1 (MQ=255)

                                                 gCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGAt                                        >  1:1523910/1‑93 (MQ=255)

                                                 gCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGAt                                        <  1:2303506/93‑1 (MQ=255)

                                                 gCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGAt                                        <  1:3532495/93‑1 (MQ=255)

                                                   aCGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGc                                      >  2:589678/1‑93 (MQ=255)

                                                   aCGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGc                                      >  1:1192221/1‑93 (MQ=255)

                                                     gTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcg                                    >  1:1465056/1‑93 (MQ=255)

                                                       aaaaaTAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcgcg                                  >  2:2548428/1‑93 (MQ=255)

                                                          aataatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGaaa                               >  1:1722881/1‑93 (MQ=255)

                                                            aaatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAAt                             <  2:3927538/92‑1 (MQ=255)

                                                             aatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATg                            >  2:3982288/1‑93 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           >  1:2285202/1‑93 (MQ=255)

                                                               tCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGt                          <  2:1623345/93‑1 (MQ=255)

                                                               tCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGt                          <  2:1623346/93‑1 (MQ=255)

                                                               tCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGt                          <  1:2575644/93‑1 (MQ=255)

                                                                ccGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTg                         <  2:2182824/93‑1 (MQ=255)

                                                                 cGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGc                        <  2:2640979/93‑1 (MQ=255)

                                                                 cGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGc                        <  2:354408/93‑1 (MQ=255)

                                                                  gCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCa                       >  2:4538837/1‑93 (MQ=255)

                                                                  gCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCa                       <  1:221207/93‑1 (MQ=255)

                                                                   cAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAg                      <  2:4591563/93‑1 (MQ=255)

                                                                   cAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAg                      <  1:599915/93‑1 (MQ=255)

                                                                     ggCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGc                    >  2:4641133/1‑93 (MQ=255)

                                                                      gCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCg                   <  2:260219/93‑1 (MQ=255)

                                                                       ccaccaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGa                  <  1:808660/93‑1 (MQ=255)

                                                                         accaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaac                >  1:652264/1‑92 (MQ=255)

                                                                         accaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaaa                <  2:3839801/93‑1 (MQ=255)

                                                                            aTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTACAACTGATGCTGCGCAAATGGTGCAGGCGAAAccc             >  2:2586252/1‑93 (MQ=255)

                                                                            aTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAccc             <  2:4363626/93‑1 (MQ=255)

                                                                            aTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAccc             <  1:3426172/93‑1 (MQ=255)

                                                                            aTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAccc             >  1:402113/1‑93 (MQ=255)

                                                                             ttGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCt            >  1:1833587/1‑93 (MQ=255)

                                                                               gATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGc          <  2:3341149/93‑1 (MQ=255)

                                                                                atatCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCt         >  1:1569169/1‑93 (MQ=255)

                                                                                atatCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCt         >  2:2660110/1‑93 (MQ=255)

                                                                                 tatCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTa        >  2:1694562/1‑93 (MQ=255)

                                                                                  atCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTAt       >  1:4252815/1‑93 (MQ=255)

                                                                                    cGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGt     <  1:2216980/93‑1 (MQ=255)

                                                                                    cGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGt     <  1:3285599/93‑1 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   >  1:460285/1‑93 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   >  1:2430032/1‑93 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   >  1:4738639/1‑93 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   <  1:309298/93‑1 (MQ=255)

                                                                                       tCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGa  <  2:1537021/93‑1 (MQ=255)

CCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGA  >  NC_000913/3019735‑3019914

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