BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,154,920	C→T	P192S (CCG→TCG)	tmk →	thymidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,154,920	0	C	T	99.1%	810.9 / ‑3.9	234	P192S (CCG→TCG)	tmk	thymidylate kinase
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (123/108); total (125/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 245 total aligned reads displayed.

GAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACGA  >  NC_000913/1154830‑1155008

gAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCg                                                                                        >  2:6519758/1‑93 (MQ=255)

gAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCg                                                                                        >  2:4354260/1‑93 (MQ=255)

 aaTCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGc                                                                                       <  1:4600091/93‑1 (MQ=255)

 aaTCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGc                                                                                       >  2:7338271/1‑93 (MQ=255)

  aTCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCt                                                                                      <  1:4059301/93‑1 (MQ=255)

  aTCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCt                                                                                      <  1:3987129/93‑1 (MQ=255)

   tCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTg                                                                                     >  1:2878189/1‑93 (MQ=255)

   tCTTTCCATTTCTTTAATTGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTg                                                                                     <  2:7169817/93‑1 (MQ=255)

     tttCGATTTCTTTAATCGCCCCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGa                                                                                   <  1:6525481/93‑1 (MQ=255)

     tttCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGa                                                                                   <  2:3626251/93‑1 (MQ=255)

     tttCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGa                                                                                   <  1:3296493/93‑1 (MQ=255)

      ttCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAg                                                                                  <  1:919578/93‑1 (MQ=255)

        cGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGc                                                                                >  1:621898/1‑93 (MQ=255)

        cGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCACTCGCTGGAGGc                                                                                >  2:2541011/1‑93 (MQ=255)

         gATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGcc                                                                               <  2:3792760/93‑1 (MQ=255)

         gATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGcc                                                                               >  2:8155980/1‑93 (MQ=255)

          aTTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCg                                                                              >  1:6300351/1‑93 (MQ=255)

           tttCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGt                                                                             <  2:3051051/93‑1 (MQ=255)

            ttCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTg                                                                            <  1:8580026/93‑1 (MQ=255)

             tCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTga                                                                           <  1:4127978/93‑1 (MQ=255)

               tttAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTgatg                                                                         >  1:5867293/1‑93 (MQ=255)

                 tAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTgatgga                                                                       >  2:4449083/1‑93 (MQ=255)

                 tAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTgatgga                                                                       >  1:748535/1‑93 (MQ=255)

                 tAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTgatgga                                                                       >  1:1838910/1‑93 (MQ=255)

                   aTCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTgatggatg                                                                     >  2:3596518/1‑93 (MQ=255)

                      gCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCaa                                                                  <  2:2649227/93‑1 (MQ=255)

                       cACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAAt                                                                 <  1:2938481/93‑1 (MQ=255)

                          ccGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCg                                                              >  1:4241564/1‑93 (MQ=255)

                           cgcgCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGc                                                             >  1:1031205/1‑93 (MQ=255)

                            gcgcCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCa                                                            >  1:4693922/1‑93 (MQ=255)

                            gcgcCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGAAATCCGCa                                                            >  2:2985430/1‑93 (MQ=255)

                             cgcCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCAc                                                           <  2:8658462/93‑1 (MQ=255)

                               cccGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTa                                                         <  1:6036351/93‑1 (MQ=255)

                                ccGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCTCCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAc                                                        <  2:1201996/93‑1 (MQ=255)

                                ccGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAc                                                        >  1:83162/1‑93 (MQ=255)

                                ccGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAc                                                        >  1:5713255/1‑93 (MQ=255)

                                 cGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAcc                                                       >  1:8157812/1‑93 (MQ=255)

                                 cGCTATCTGGAACTGGCAGCACAAGATAAAAGCATGCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAcc                                                       <  2:8496019/93‑1 (MQ=255)

                                  gCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCg                                                      >  2:8766670/1‑93 (MQ=255)

                                  gCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCg                                                      >  1:2609974/1‑93 (MQ=255)

                                    tATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGta                                                    >  2:790973/1‑92 (MQ=255)

                                    tATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTg                                                    >  1:5380975/1‑93 (MQ=255)

                                      tCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGAc                                                  <  2:4205335/93‑1 (MQ=255)

                                      tCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGAc                                                  >  1:5140904/1‑93 (MQ=255)

                                      tCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGAc                                                  >  1:3828651/1‑93 (MQ=255)

                                         ggAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCa                                               <  1:5266867/93‑1 (MQ=255)

                                          gAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCAc                                              <  1:4814537/93‑1 (MQ=255)

                                           aaCTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACt                                             >  1:2714241/1‑93 (MQ=255)

                                             cTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTgg                                           >  2:137681/1‑93 (MQ=255)

                                             cTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTgg                                           <  1:7082601/93‑1 (MQ=255)

                                              tGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTggg                                          <  2:554581/93‑1 (MQ=255)

                                                 cagcaCAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGa                                       <  1:3559049/93‑1 (MQ=255)

                                                  agcaCAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGGACTACCGTGACCCACTGGGTGaa                                      >  2:4766755/1‑93 (MQ=255)

                                                  agcaCAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGaa                                      <  2:7136135/93‑1 (MQ=255)

                                                  agcaCAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGaa                                      <  1:51865/93‑1 (MQ=255)

                                                   gcaCAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAg                                     <  2:6151975/93‑1 (MQ=255)

                                                    cacaAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATACGCACTACCGTGACCCACTGGGGGAAgg                                    >  2:5461358/1‑93 (MQ=255)

                                                     acaAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGa                                   >  2:1941995/1‑93 (MQ=255)

                                                     acaAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGa                                   >  1:7077852/1‑93 (MQ=255)

                                                     acaAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGa                                   >  1:1436646/1‑93 (MQ=255)

                                                      caAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAg                                  >  1:2517104/1‑93 (MQ=255)

                                                        agagAAAATCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGtt                                <  2:3960892/93‑1 (MQ=255)

                                                        aGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGtt                                >  1:4525081/1‑93 (MQ=255)

                                                         gATAAAAGCATTCATACCATTGATGCCACCCCGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGTAGGAGTTg                               >  2:5237984/1‑93 (MQ=255)

                                                         gATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTg                               >  1:3685289/1‑93 (MQ=255)

                                                           tAAAAGCATTCATACCATTGATGCCTCCCAGTCGCTGGAGACCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGa                             <  2:1446329/93‑1 (MQ=255)

                                                            aaaaGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGAc                            <  2:1793346/93‑1 (MQ=255)

                                                            aaaaGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGAc                            <  1:1878710/93‑1 (MQ=255)

                                                             aaaGCATTCATACCATTGATGCCACCCAGTCGCTGTAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACg                           >  2:6996652/1‑93 (MQ=255)

                                                             aaaGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACg                           >  1:2041638/1‑93 (MQ=255)

                                                             aaaGCATTAATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACg                           >  2:5358973/1‑93 (MQ=255)

                                                                gCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCAt                        >  2:7263597/1‑93 (MQ=255)

                                                                acATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGCGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGGGTTGGACGCAt                        <  2:2093327/92‑1 (MQ=255)

                                                                  aTTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATga                      <  1:1172244/93‑1 (MQ=255)

                                                                   ttCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATgag                     >  2:3165298/1‑93 (MQ=255)

                                                                    tCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATgaga                    <  1:8071456/93‑1 (MQ=255)

                                                                     cATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGAt                   <  1:4185073/93‑1 (MQ=255)

                                                                      aTACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATg                  <  1:6615201/93‑1 (MQ=255)

                                                                       tACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATgg                 <  2:8248300/93‑1 (MQ=255)

                                                                       tACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGAACAACCTTGACCCACTGGGGGAAGGAGTTGGACGCATGGGATgg                 >  2:6385767/1‑93 (MQ=255)

                                                                        aCCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGt                <  2:1513615/93‑1 (MQ=255)

                                                                         ccATTGATGCCACCCAGTCGCTGGAGGCGGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTa               <  1:3388693/93‑1 (MQ=255)

                                                                         ccATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTa               >  2:2325676/1‑93 (MQ=255)

                                                                         ccATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTa               <  2:5048056/93‑1 (MQ=255)

                                                                         ccATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTa               >  2:6701378/1‑93 (MQ=255)

                                                                           aTTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATc             >  2:2340160/1‑93 (MQ=255)

                                                                           aTTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATc             >  1:221452/1‑93 (MQ=255)

                                                                           aTTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGCGCGAAGGAGGTAGAGGCATGAGATGGGATc             >  2:6808308/1‑93 (MQ=255)

                                                                             tgctGCCGCCCAGTCGCTGGAGGTCGTTATGGATGCAATCCGCCCTACCGTGACCCCCTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCa           <  2:5975803/90‑1 (MQ=255)

                                                                             tGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCa           <  2:1702339/93‑1 (MQ=255)

                                                                             tGAAGCCACCCAGACGCTGGAGGCCGTGATGGATGCAATCCGCTCTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCa           <  1:4956261/93‑1 (MQ=255)

                                                                              gATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCAt          <  1:3896977/93‑1 (MQ=255)

                                                                               aTGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATg         <  1:4893434/93‑1 (MQ=255)

                                                                               aTGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATg         <  1:7887992/93‑1 (MQ=255)

                                                                                tGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATgg        <  1:6798276/93‑1 (MQ=255)

                                                                                  ccACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGtt      >  1:5763661/1‑93 (MQ=255)

                                                                                    aCCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTAc    >  1:5603776/1‑93 (MQ=255)

                                                                                     cccAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACg   >  1:580375/1‑93 (MQ=255)

                                                                                      ccAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACGa  >  1:8714857/1‑93 (MQ=255)

                                                                                      ccAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACGa  >  1:4702540/1‑93 (MQ=255)

GAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACGA  >  NC_000913/1154830‑1155008

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