BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,472,339	(C)_6→7	coding (95/861 nt)	ydbC →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,472,333	1	.	C	98.6%	440.8 / ‑1.0	140	coding (89/861 nt)	ydbC	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (65/73); total (65/75)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 167 total aligned reads displayed.

TGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTA  >  NC_000913/1472241‑1472420

   gATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑ccc                                                                                      >  2:1002567/1‑93 (MQ=255)

   gATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑ccc                                                                                      >  2:4677888/1‑93 (MQ=255)

    aTGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑cccc                                                                                     >  2:6115903/1‑93 (MQ=255)

    aTGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTCG‑cccc                                                                                     >  2:806119/1‑93 (MQ=255)

     tGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑ccccc                                                                                    >  2:2314851/1‑93 (MQ=255)

       agcagcAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑ccccccc                                                                                  <  1:6245057/93‑2 (MQ=255)

       agcagcAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑ccccccc                                                                                  >  2:3975697/1‑92 (MQ=255)

        gcagcaATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTTGAGTTTTTGG‑ccccccca                                                                                 >  1:6055877/1‑91 (MQ=255)

        gcagcaATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑ccccccca                                                                                 >  2:319828/1‑91 (MQ=255)

        gcagcaATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑ccccccca                                                                                 >  1:2878488/1‑91 (MQ=255)

         cagcaATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑cccccccac                                                                                <  2:3373431/93‑4 (MQ=255)

          agcaATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑cccccccacg                                                                               <  1:234514/93‑5 (MQ=255)

          agcaATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑cccccccacg                                                                               >  1:4838492/1‑89 (MQ=255)

            caATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACgag                                                                              >  1:3521531/1‑93 (MQ=255)

             aaTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACgaga                                                                             >  1:2686579/1‑93 (MQ=255)

              aTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGAt                                                                            >  1:1005899/1‑93 (MQ=255)

              aTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGAt                                                                            >  2:148983/1‑93 (MQ=255)

              aTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGAt                                                                            <  2:3026842/93‑1 (MQ=255)

              aTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGAt                                                                            <  2:2588171/93‑1 (MQ=255)

              aTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGAt                                                                            <  2:6069680/93‑1 (MQ=255)

                 cATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGc                                                                         <  1:2004417/93‑1 (MQ=255)

                  aTTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGcc                                                                        >  1:1703225/1‑93 (MQ=255)

                  aTTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGcc                                                                        >  1:5595603/1‑93 (MQ=255)

                   tttACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGccg                                                                       >  2:485147/1‑92 (MQ=255)

                   tttACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGccg                                                                       >  1:2568119/1‑92 (MQ=255)

                   tttACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGccg                                                                       <  1:4543175/93‑2 (MQ=255)

                      aCTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCcgcgc                                                                    >  2:4302160/1‑92 (MQ=255)

                       ctctCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTc                                                                   >  2:667798/1‑93 (MQ=255)

                        tctcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCg                                                                  >  2:422128/1‑93 (MQ=255)

                        tctcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCg                                                                  >  2:1472714/1‑93 (MQ=255)

                          tcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCt                                                                >  2:3842556/1‑93 (MQ=255)

                          tcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCt                                                                >  2:5866379/1‑93 (MQ=255)

                           cGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCta                                                               <  1:2752817/93‑1 (MQ=255)

                           cGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCta                                                               <  1:2476766/93‑1 (MQ=255)

                            ggTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCtat                                                              >  2:4069966/1‑93 (MQ=255)

                            ggTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCtat                                                              >  1:5308121/1‑93 (MQ=255)

                               aCAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAAc                                                           >  2:506064/1‑93 (MQ=255)

                               aCAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAAc                                                           <  1:1380651/93‑1 (MQ=255)

                                cAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAAcc                                                          <  2:616934/93‑1 (MQ=255)

                                 aaaaTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCg                                                         >  1:680375/1‑93 (MQ=255)

                                 aaaaTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCg                                                         <  1:6113534/93‑1 (MQ=255)

                                 aaaaTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCg                                                         >  2:2318360/1‑93 (MQ=255)

                                   aaTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGtg                                                       >  1:4732056/1‑93 (MQ=255)

                                     tCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGtgct                                                     <  1:3446510/93‑1 (MQ=255)

                                     tCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGtgct                                                     <  1:3921124/93‑1 (MQ=255)

                                      ccGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGtgctg                                                    <  2:1046529/93‑1 (MQ=255)

                                       cGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGtgctgc                                                   <  1:2490490/93‑1 (MQ=255)

                                        gTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCg                                                  <  2:5626054/93‑1 (MQ=255)

                                        gTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCg                                                  <  1:102742/93‑1 (MQ=255)

                                           aaCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGa                                               >  1:648837/1‑93 (MQ=255)

                                           aaCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGa                                               <  2:4578497/93‑1 (MQ=255)

                                           aaCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGa                                               <  2:4205829/93‑1 (MQ=255)

                                           aaCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCCGCGTGa                                               <  2:5237743/93‑1 (MQ=255)

                                            aCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAg                                              >  1:2827167/1‑93 (MQ=255)

                                            aCCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAg                                              <  1:3565707/93‑1 (MQ=255)

                                             ccGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAgg                                             >  1:3258831/1‑93 (MQ=255)

                                              cGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgc                                            >  1:694770/1‑93 (MQ=255)

                                              cGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgc                                            <  1:1582647/93‑1 (MQ=255)

                                              cGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgc                                            <  1:3949871/93‑1 (MQ=255)

                                                 cTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCt                                         <  1:3508016/93‑1 (MQ=255)

                                                  ttGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTg                                        >  2:6256063/1‑93 (MQ=255)

                                                  ttGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTg                                        >  2:5441137/1‑93 (MQ=255)

                                                  ttGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTg                                        >  2:3268797/1‑93 (MQ=255)

                                                    ggTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGc                                      >  1:5276245/1‑93 (MQ=255)

                                                    ggTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGc                                      >  2:1210452/1‑93 (MQ=255)

                                                       tATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCAtt                                   <  2:1863786/93‑1 (MQ=255)

                                                       tATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCAtt                                   <  2:1844908/93‑1 (MQ=255)

                                                         tGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTgg                                 >  2:84512/1‑93 (MQ=255)

                                                           gcgcGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGc                               <  1:1548713/93‑1 (MQ=255)

                                                           gcgcGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGc                               >  1:4696298/1‑93 (MQ=255)

                                                                aTGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCaa                          <  1:6294046/93‑1 (MQ=255)

                                                                aTGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCaa                          >  1:3770376/1‑93 (MQ=255)

                                                                aTGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCGCGCATTGGTCGTCaa                          >  1:2402475/1‑93 (MQ=255)

                                                                 tGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAAt                         >  2:5200916/1‑93 (MQ=255)

                                                                   cAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCa                       <  2:1151514/93‑1 (MQ=255)

                                                                   cAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCa                       <  2:5941582/93‑1 (MQ=255)

                                                                   cAACTGGCAGGCCCTGGAGTTTTTGGCCCCCCCACGAGGTCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCa                       >  2:238057/1‑93 (MQ=255)

                                                                    aaCTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCat                      >  1:2016942/1‑93 (MQ=255)

                                                                     aCTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCata                     <  1:4828229/93‑1 (MQ=255)

                                                                     aCTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCata                     <  1:4904031/93‑1 (MQ=255)

                                                                      cTGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCatat                    <  2:4646331/93‑1 (MQ=255)

                                                                       tGGCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATAtt                   <  2:6146569/93‑1 (MQ=255)

                                                                         gCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGa                 >  1:5823889/1‑93 (MQ=255)

                                                                         gCAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGa                 <  1:3281489/93‑1 (MQ=255)

                                                                          cAGGTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGAt                <  1:941992/93‑1 (MQ=255)

                                                                             gTCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATAcc             <  2:4688888/93‑1 (MQ=255)

                                                                              tCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa            >  1:4492920/1‑93 (MQ=255)

                                                                              tCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa            <  2:5815575/93‑1 (MQ=255)

                                                                              tCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa            >  2:2628924/1‑93 (MQ=255)

                                                                              tCCTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa            <  2:680375/93‑1 (MQ=255)

                                                                                cTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGc          <  2:5313065/93‑1 (MQ=255)

                                                                                cTGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCAGTGGGCGTCAATCATATTGATACCAGc          >  1:6304687/1‑93 (MQ=255)

                                                                                 tGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCATTCATATTGATACCAGCg         >  2:6058526/1‑93 (MQ=255)

                                                                                 tGGAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCg         >  2:5718016/1‑93 (MQ=255)

                                                                                  ggAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGa        >  2:1503227/1‑93 (MQ=255)

                                                                                  ggAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGa        >  2:1693068/1‑93 (MQ=255)

                                                                                   gAGTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGAc       <  1:3351764/93‑1 (MQ=255)

                                                                                     gTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACtt     <  1:3081316/93‑1 (MQ=255)

                                                                                     gTTTTTGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACtt     >  1:3815186/1‑93 (MQ=255)

                                                                                        tttGGCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTa  <  2:2227021/93‑1 (MQ=255)

TGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑CCCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTA  >  NC_000913/1472241‑1472420

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