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.6%	868.2 / ‑5.6	238	P163P (CCG→CCA)	b2881	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (99/138); minor base C (1/0); total (100/138)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 247 total aligned reads displayed.

CCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGA  >  NC_000913/3019735‑3019914

 cgcgATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACt                                                                                        <  1:1990509/93‑1 (MQ=255)

   cgATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACttt                                                                                      <  2:3000531/93‑1 (MQ=255)

   cgATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACttt                                                                                      >  2:5580205/1‑93 (MQ=255)

    gATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTc                                                                                     >  2:2175987/1‑93 (MQ=255)

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

    gATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACCTTc                                                                                     >  1:3520667/1‑93 (MQ=255)

     aTGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTcc                                                                                    <  2:4156691/93‑1 (MQ=255)

     aTGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTcc                                                                                    >  2:600242/1‑93 (MQ=255)

     aTGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTcc                                                                                    <  1:4532274/93‑1 (MQ=255)

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

       gctggctgGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGt                                                                                  <  1:1188389/93‑1 (MQ=255)

       gctggctgGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGt                                                                                  <  2:3370407/93‑1 (MQ=255)

       gctggctgGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGt                                                                                  <  1:428018/93‑1 (MQ=255)

       gctgactgGCAGCAATACTATCAGGTCATTGAACAGGCAGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGt                                                                                  <  2:6297432/93‑1 (MQ=255)

         tggctggcAGCAATACTATCAGGTCATTTAACTGTCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTga                                                                                <  2:2327915/93‑1 (MQ=255)

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

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

             tggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacgac                                                                            <  2:1341879/93‑1 (MQ=255)

             tggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAAGCCACCATTGATATCGCTTCAACTTTCCGTgacgac                                                                            >  2:4919863/1‑93 (MQ=255)

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

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

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

                cagcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTa                                                                         <  2:4530374/93‑1 (MQ=255)

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

                  gcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGa                                                                       >  2:1155629/1‑93 (MQ=255)

                    aaTACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAg                                                                     >  1:2348292/1‑93 (MQ=255)

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

                      tACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTc                                                                   >  2:4716603/1‑93 (MQ=255)

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

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

                        cTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAt                                                                 >  1:4109800/1‑93 (MQ=255)

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

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

                             aGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTa                                                            <  2:5800621/93‑1 (MQ=255)

                             aGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTa                                                            >  1:2949875/1‑93 (MQ=255)

                             aGGTCATTGAACTGGCGGTTGCAAGTAAAAATAATCCGCAGGCCACGATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTa                                                            >  2:5075065/1‑93 (MQ=255)

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

                              ggTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTaa                                                           >  1:573099/1‑93 (MQ=255)

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

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

                                 cATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCa                                                        <  1:2659408/93‑1 (MQ=255)

                                   ttGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAtt                                                      >  2:3564918/1‑93 (MQ=255)

                                     gAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTac                                                    >  1:1087824/1‑92 (MQ=255)

                                      aaCTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATc                                                   <  2:2391169/93‑1 (MQ=255)

                                      aaCTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATc                                                   >  2:3922129/1‑93 (MQ=255)

                                        cTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCt                                                 <  1:1793700/93‑1 (MQ=255)

                                         tGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTa                                                <  1:342436/93‑1 (MQ=255)

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

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

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

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

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

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

                                            cGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTaaaa                                             >  2:1106233/1‑93 (MQ=255)

                                             ggTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAAc                                            >  2:3798924/1‑93 (MQ=255)

                                             ggTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAAc                                            >  2:6149980/1‑93 (MQ=255)

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

                                              gTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACt                                           >  1:5631214/1‑93 (MQ=255)

                                              gTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACt                                           >  1:5182225/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                    cGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGcc                                     <  1:3032295/93‑1 (MQ=255)

                                                      tAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcgc                                   >  2:2800806/1‑93 (MQ=255)

                                                      tAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcgc                                   >  2:1243086/1‑93 (MQ=255)

                                                      tAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcgc                                   <  2:1205829/93‑1 (MQ=255)

                                                      tAAAAATAACCCGCAGGCCACCATTGATATCGCTCCAACTTCCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcgc                                   <  1:673165/93‑1 (MQ=255)

                                                           ataatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGaaaa                              <  1:5346360/93‑1 (MQ=255)

                                                            taatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAAt                             <  1:2370444/93‑1 (MQ=255)

                                                            taatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAAt                             >  2:1761576/1‑93 (MQ=255)

                                                             aatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATg                            <  2:193398/93‑1 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           <  1:2570713/93‑1 (MQ=255)

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

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

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

                                                                ccGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTg                         >  1:6002462/1‑93 (MQ=255)

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

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

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

                                                                   cAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGTTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAg                      <  2:4066004/93‑1 (MQ=255)

                                                                    aGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAgg                     >  1:1309948/1‑93 (MQ=255)

                                                                    aGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAgg                     <  2:2748772/93‑1 (MQ=255)

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

                                                                      gCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCg                   >  2:877403/1‑93 (MQ=255)

                                                                        caccaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaa                 >  2:4858880/1‑93 (MQ=255)

                                                                         accaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGGCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaaa                >  2:5403138/1‑93 (MQ=255)

                                                                          ccaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAc               >  1:4502975/1‑93 (MQ=255)

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

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

                                                                              tGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTg           >  2:2224677/1‑93 (MQ=255)

                                                                               aatatCGCTCAAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAGCCCTGc          <  2:5197322/92‑1 (MQ=255)

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

                                                                                  atCGCTCGAACTGTCCGTGACGACCTAGAAGTCATTGGTAAGCATTAGCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTAt       <  1:3975557/93‑1 (MQ=255)

                                                                                   tCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATg      <  2:4960926/93‑1 (MQ=255)

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

                                                                                       tCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGa  >  1:2280238/1‑93 (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: `‑`