BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,525,330	T→C	T408T (ACT→ACC)	yrfF →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,525,330	0	T	C	99.2%	1128.4 / ‑2.2	365	T408T (ACT→ACC)	yrfF	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (168/194); total (170/195)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 378 total aligned reads displayed.

AGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATT  >  NC_000913/3525240‑3525419

aGCCACCAGCTTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTc                                                                                         <  1:4231571/93‑1 (MQ=255)

 gCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTcc                                                                                        >  1:1892168/1‑93 (MQ=255)

  ccaccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCg                                                                                       >  1:4859713/1‑93 (MQ=255)

   caccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCgg                                                                                      >  2:5164880/1‑93 (MQ=255)

    accaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGa                                                                                     >  2:1567972/1‑93 (MQ=255)

     ccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGaa                                                                                    <  1:1096867/93‑1 (MQ=255)

     ccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGaa                                                                                    >  1:2806947/1‑93 (MQ=255)

     ccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGaa                                                                                    >  1:2528447/1‑93 (MQ=255)

      caGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAAc                                                                                   >  2:2137706/1‑93 (MQ=255)

        gCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCt                                                                                 <  1:6266449/93‑1 (MQ=255)

         cGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTg                                                                                <  2:6123167/93‑1 (MQ=255)

           tAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGa                                                                              >  2:4575150/1‑93 (MQ=255)

           tAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGa                                                                              <  1:928646/93‑1 (MQ=255)

             aaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgc                                                                            <  1:7839487/93‑1 (MQ=255)

             aaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgc                                                                            <  1:3338362/93‑1 (MQ=255)

             aaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgc                                                                            <  2:4447684/93‑1 (MQ=255)

             aaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgc                                                                            >  2:350812/1‑93 (MQ=255)

             aaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgc                                                                            <  2:2578967/93‑1 (MQ=255)

              aacaacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgcg                                                                           >  1:4726534/1‑93 (MQ=255)

               acaacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgcgc                                                                          >  1:3001378/1‑93 (MQ=255)

               acaacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgcgc                                                                          >  1:168934/1‑93 (MQ=255)

                caacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGcgcg                                                                         >  2:2342650/1‑93 (MQ=255)

                 aacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGa                                                                        <  1:3904960/93‑1 (MQ=255)

                  acTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaa                                                                       <  1:1463234/93‑1 (MQ=255)

                  acTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaa                                                                       >  2:3060113/1‑93 (MQ=255)

                   cTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaaa                                                                      <  1:979979/93‑1 (MQ=255)

                      gCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAAcc                                                                   <  1:3124936/93‑1 (MQ=255)

                      gCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAAcc                                                                   <  2:2513272/93‑1 (MQ=255)

                      gCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAAcc                                                                   <  2:5224202/93‑1 (MQ=255)

                       cTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCa                                                                  <  2:556509/93‑1 (MQ=255)

                        tGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCaa                                                                 >  1:3519857/1‑93 (MQ=255)

                           cGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTc                                                              >  1:6503825/1‑93 (MQ=255)

                            gCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAAttaa                                                             >  2:479590/1‑91 (MQ=255)

                            gCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCa                                                             <  2:7217862/93‑1 (MQ=255)

                              tGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCAcc                                                           >  2:3670822/1‑93 (MQ=255)

                                gCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCtt                                                         <  2:5816718/93‑1 (MQ=255)

                                   tGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCttttt                                                      <  2:6275361/93‑1 (MQ=255)

                                     cGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTa                                                    >  2:5070265/1‑93 (MQ=255)

                                      gggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAc                                                   <  1:4548531/93‑1 (MQ=255)

                                      gggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAc                                                   >  2:5397641/1‑93 (MQ=255)

                                      gggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAc                                                   >  1:1248264/1‑93 (MQ=255)

                                      gggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAc                                                   <  2:1965102/93‑1 (MQ=255)

                                      gggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGCAATCTTCGAACCTCCGGGACCTGGAGCGCGAAAACCAATTCACCTTTTTTAc                                                   >  2:6531047/1‑93 (MQ=255)

                                        gTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTCCCg                                                 <  2:2442883/93‑1 (MQ=255)

                                        gTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCg                                                 >  2:4107235/1‑93 (MQ=255)

                                         tAGTCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGt                                                <  2:7294875/93‑1 (MQ=255)

                                         tAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGt                                                >  1:3667387/1‑93 (MQ=255)

                                          aGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGtt                                               >  2:6886659/1‑93 (MQ=255)

                                           ggCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGttt                                              <  2:4283007/93‑1 (MQ=255)

                                           ggCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGttt                                              <  1:6885700/93‑1 (MQ=255)

                                           ggCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGttt                                              <  1:7536132/93‑1 (MQ=255)

                                           ggCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGttt                                              <  1:7283541/93‑1 (MQ=255)

                                            gCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTg                                             <  1:5485059/93‑1 (MQ=255)

                                              gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           <  2:4364522/93‑1 (MQ=255)

                                               aTACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACt                                          >  1:818204/1‑93 (MQ=255)

                                                tACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTg                                         >  2:5467100/1‑93 (MQ=255)

                                                 aCATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGc                                        <  1:2619644/93‑1 (MQ=255)

                                                 aCATTGCGTATTAGCGATACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGc                                        <  1:7033710/93‑1 (MQ=255)

                                                   aTTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTc                                      <  2:1839593/93‑1 (MQ=255)

                                                    ttGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACGGTTTGACTGCTCg                                     >  1:6384959/1‑93 (MQ=255)

                                                    ttGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCg                                     >  2:3839931/1‑93 (MQ=255)

                                                    ttGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCg                                     >  2:228685/1‑93 (MQ=255)

                                                        gTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGa                                 >  1:5616909/1‑93 (MQ=255)

                                                        gTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGa                                 <  1:5316862/93‑1 (MQ=255)

                                                        gTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGa                                 >  2:1283880/1‑93 (MQ=255)

                                                         tATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGat                                >  1:1032304/1‑93 (MQ=255)

                                                           ttAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatca                              <  1:6654525/93‑1 (MQ=255)

                                                           ttAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatca                              >  2:3534235/1‑93 (MQ=255)

                                                              gCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCt                           <  2:7545774/93‑1 (MQ=255)

                                                               cGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTg                          >  1:2231935/1‑93 (MQ=255)

                                                               cGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTg                          <  2:584613/93‑1 (MQ=255)

                                                                 gTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGa                        <  2:5975000/93‑1 (MQ=255)

                                                                 gTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGa                        >  1:3242922/1‑93 (MQ=255)

                                                                 gTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGa                        >  2:167499/1‑93 (MQ=255)

                                                                  tACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGaa                       >  2:26508/1‑93 (MQ=255)

                                                                  tACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGaa                       <  1:3292768/93‑1 (MQ=255)

                                                                  tACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGaa                       >  1:7450380/1‑93 (MQ=255)

                                                                   aCGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAac                      <  2:4892537/93‑1 (MQ=255)

                                                                     gggAATGTGTAATATTCGAGCTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacga                    <  1:5222068/93‑1 (MQ=255)

                                                                     gggAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacga                    >  2:2774533/1‑93 (MQ=255)

                                                                      ggAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacgcc                   >  2:6798313/1‑91 (MQ=255)

                                                                      ggAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacgac                   >  2:7403142/1‑93 (MQ=255)

                                                                       gAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacgacg                  >  1:4076464/1‑93 (MQ=255)

                                                                         aTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGcc                <  1:5544625/93‑1 (MQ=255)

                                                                          tgtgTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGccc               >  1:3811445/1‑93 (MQ=255)

                                                                          tgtgTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGccc               <  2:3351729/93‑1 (MQ=255)

                                                                          tgtgTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGccc               >  2:7763536/1‑93 (MQ=255)

                                                                          tgtgTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGccc               <  2:7165618/93‑1 (MQ=255)

                                                                          tgtgTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGccc               <  1:6817277/93‑1 (MQ=255)

                                                                            tgtAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGc             <  2:1151009/93‑1 (MQ=255)

                                                                              tAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTc           >  1:4599581/1‑93 (MQ=255)

                                                                               aaTATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCa          <  1:6115431/93‑1 (MQ=255)

                                                                               aaTATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCa          >  1:5947668/1‑93 (MQ=255)

                                                                                atatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAt         <  2:7917501/93‑1 (MQ=255)

                                                                                 tatTCGACCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAtt        <  2:1739016/93‑1 (MQ=255)

                                                                                 tatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAtt        <  1:2074388/93‑1 (MQ=255)

                                                                                 aatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAtt        <  1:5744489/92‑1 (MQ=255)

                                                                                    tCGAGCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGcc     <  2:6160369/93‑1 (MQ=255)

                                                                                     cGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCa    >  1:504197/1‑93 (MQ=255)

                                                                                      gAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAt   >  1:2171122/1‑93 (MQ=255)

AGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATT  >  NC_000913/3525240‑3525419

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