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%	871.5 / NA	240	P163P (CCG→CCA)	b2881	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (96/142); total (97/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

CCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGA  >  NC_000913/3019735‑3019914

ccGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGGCACCATTGATATCGCTCCCAc                                                                                         >  1:536133/1‑93 (MQ=255)

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

 cgcgATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGAAGGCCACCATTGATATCGCTCCAACt                                                                                        <  2:6263572/93‑1 (MQ=255)

  gcgATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACtt                                                                                       <  1:1257580/93‑1 (MQ=255)

  gcgATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACtt                                                                                       >  2:4419456/1‑93 (MQ=255)

  gcgATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACGGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACtt                                                                                       <  1:1537014/93‑1 (MQ=255)

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

   cgATGCTGGCTGGCAGCAATACTATCAAGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACttt                                                                                      >  2:5453131/1‑93 (MQ=255)

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

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

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

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

        ctggctggCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTg                                                                                 >  2:5338235/1‑93 (MQ=255)

         tggctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCGCCATTGATATCGCTCCAACTTTCCGTga                                                                                <  2:1321321/93‑1 (MQ=255)

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

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

            ctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGcgacga                                                                             <  2:1546250/93‑1 (MQ=255)

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

            ctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAAGCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacga                                                                             <  1:4448173/93‑1 (MQ=255)

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

              ggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGAcc                                                                           >  1:6229581/1‑93 (MQ=255)

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

              ggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGACACCATTGATATCGCTCCAACTTTCCGTGACGAcc                                                                           >  2:6570420/1‑93 (MQ=255)

              ggcAGCAATACTATCAGGTCATTGAACTGGCGGATGCACGTAAAAATAATCCGCAGGCCAACATTGATATCGCTCCAACTTTCCGTGACGAcc                                                                           >  2:4106711/1‑93 (MQ=255)

              ggcACCAATACTATCAGGTCATTGAACTGGCGGTTGCAAGTAAAAATAATCCGCAGGCCACCATTGAAATCGCTCCAACTTTCCGTGACAAcc                                                                           <  1:4637077/93‑1 (MQ=255)

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

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

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

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

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

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

                    aaTACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTGCCGTGACGACCTAGAAg                                                                     <  1:4903320/93‑1 (MQ=255)

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

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

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

                         tATCAGGTCATTGAACTGGCGGTTGAACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAtt                                                                <  2:6480018/93‑1 (MQ=255)

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

                          aTCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTg                                                               <  1:4563932/93‑1 (MQ=255)

                           tCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTgg                                                              <  1:4708345/93‑1 (MQ=255)

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

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

                               gTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAg                                                          >  2:193216/1‑93 (MQ=255)

                                tCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGc                                                         >  2:4309162/1‑93 (MQ=255)

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

                                  aTTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAt                                                       <  1:832282/93‑1 (MQ=255)

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

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

                                   ttGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAtt                                                      <  2:4209064/93‑1 (MQ=255)

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

                                       aCTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATcc                                                  <  1:2387243/93‑1 (MQ=255)

                                        cTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCt                                                 >  2:71277/1‑93 (MQ=255)

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

                                               ttGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTg                                          >  1:1485637/1‑93 (MQ=255)

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

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

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

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

                                                  cACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATg                                       >  2:6726617/1‑93 (MQ=255)

                                                   aCGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGc                                      <  2:611224/93‑1 (MQ=255)

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

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

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

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

                                                            taatCCGCAGGCCGCCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAAt                             <  1:3303112/93‑1 (MQ=255)

                                                             actCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATg                            <  1:646698/91‑1 (MQ=255)

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

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

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

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

                                                                 cGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTga                        >  2:1375802/1‑92 (MQ=255)

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

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

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

                                                                   cAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAg                      >  2:2921748/1‑93 (MQ=255)

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

                                                                     ggCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGc                    <  2:4644731/93‑1 (MQ=255)

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

                                                                       ccaccaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGa                  >  1:5462628/1‑93 (MQ=255)

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

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

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

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

                                                                         accaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaaa                >  2:1019216/1‑93 (MQ=255)

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

                                                                          ccaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAc               <  2:5120183/93‑1 (MQ=255)

                                                                           caTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAcc              >  2:3118346/1‑93 (MQ=255)

                                                                           caTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAcc              >  2:4055083/1‑93 (MQ=255)

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

                                                                            aTTGATATCGCTCCAACCTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAccc             <  1:1425380/93‑1 (MQ=255)

                                                                             ttGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCt            <  2:584535/93‑1 (MQ=255)

                                                                             ttGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCt            <  2:5272256/93‑1 (MQ=255)

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

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

                                                                                atatCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCt         <  1:3281731/93‑1 (MQ=255)

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

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

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

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

                                                                                       tCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGa  >  2:754577/1‑93 (MQ=255)

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

CCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGA  >  NC_000913/3019735‑3019914

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