BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,339,304	+C	coding (398/987 nt)	yrbH →	putative isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,339,301	1	.	C	99.6%	869.1 / ‑4.1	268	coding (395/987 nt)	yrbH	putative isomerase
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (133/134); total (134/134)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 300 total aligned reads displayed.

TGCTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑CCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGAC  >  NC_000913/3339209‑3339388

 tcTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑c                                                                                        <  1:6460250/92‑1 (MQ=255)

 gCTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑c                                                                                        >  2:2741832/1‑93 (MQ=255)

 gCTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑c                                                                                        <  2:1987061/93‑1 (MQ=255)

  cTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑cc                                                                                       >  2:6456681/1‑93 (MQ=255)

    aTCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑cccc                                                                                     >  2:2812273/1‑92 (MQ=255)

    aTCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑cccc                                                                                     <  1:5509319/93‑2 (MQ=255)

     tctcTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccccg                                                                                    <  1:995884/93‑3 (MQ=255)

      ctctAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccccgg                                                                                   <  1:3839503/93‑4 (MQ=255)

       tctAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccccgga                                                                                  >  2:3644605/1‑89 (MQ=255)

       tctAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑ccccgga                                                                                  >  1:2621427/1‑89 (MQ=255)

         tAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGgaga                                                                                 <  1:1189834/93‑1 (MQ=255)

         tAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGgaga                                                                                 <  2:3596638/93‑1 (MQ=255)

         tAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGgaga                                                                                 <  1:72153/93‑1 (MQ=255)

          aaCTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCTGTCGCCCCGGagag                                                                                >  1:6402329/1‑93 (MQ=255)

            ctctGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGagca                                                                              >  1:4116727/1‑93 (MQ=255)

            ctctGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGagca                                                                              <  2:4407730/93‑1 (MQ=255)

            ctctGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGagca                                                                              <  2:1673327/93‑1 (MQ=255)

               tGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAgcagca                                                                           >  1:297443/1‑93 (MQ=255)

               tGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAgcagca                                                                           <  1:2544175/93‑1 (MQ=255)

               tGGGGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAgcagca                                                                           >  2:5246043/1‑93 (MQ=255)

                 gTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATg                                                                         >  1:3381364/1‑93 (MQ=255)

                 gTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATg                                                                         >  2:2016662/1‑93 (MQ=255)

                  tGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATgg                                                                        <  2:4618918/93‑1 (MQ=255)

                   gAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGgc                                                                       <  1:5639047/93‑1 (MQ=255)

                   gAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGgc                                                                       <  2:5799123/93‑1 (MQ=255)

                     aTCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGgcgc                                                                     >  1:5345424/1‑93 (MQ=255)

                      tCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGgcgcg                                                                    >  2:6557015/1‑93 (MQ=255)

                       ccAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGgcgcgc                                                                   >  1:1408970/1‑93 (MQ=255)

                       ccAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGgcgcgc                                                                   >  1:2893175/1‑93 (MQ=255)

                         aGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGgcgcgcgc                                                                 >  1:3456835/1‑93 (MQ=255)

                         aGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGgcgcgcgc                                                                 <  1:2725776/93‑1 (MQ=255)

                            gAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGcgccgc                                                              >  2:1848540/1‑93 (MQ=255)

                            gAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGcgccgc                                                              >  2:6201316/1‑93 (MQ=255)

                             aaaTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCa                                                             >  1:3596406/1‑93 (MQ=255)

                               aTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGa                                                           <  1:6610188/93‑1 (MQ=255)

                               aTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGa                                                           <  2:6061068/93‑1 (MQ=255)

                               aTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGa                                                           <  1:6176770/93‑1 (MQ=255)

                               aTCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGa                                                           <  1:2175822/93‑1 (MQ=255)

                                tCACGGCCTGAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGAt                                                          <  2:779436/93‑1 (MQ=255)

                                 cACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGAtg                                                         <  2:6360357/93‑1 (MQ=255)

                                  aCGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGAtgt                                                        >  2:2558828/1‑93 (MQ=255)

                                  aCGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGAtgt                                                        <  2:2190474/93‑1 (MQ=255)

                                   cGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGAtgtg                                                       <  1:6064380/93‑1 (MQ=255)

                                    gggcTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGc                                                      <  2:1576375/90‑1 (MQ=255)

                                     gCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCa                                                     <  1:1323321/93‑1 (MQ=255)

                                       cTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATc                                                   <  2:3746346/93‑1 (MQ=255)

                                        ttAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCt                                                  >  2:909763/1‑93 (MQ=255)

                                         tAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtg                                                 <  2:2905106/93‑1 (MQ=255)

                                         tAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtg                                                 >  1:516070/1‑93 (MQ=255)

                                          aaTTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtgt                                                >  1:1858044/1‑93 (MQ=255)

                                          aaTTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtgt                                                <  1:4722570/93‑1 (MQ=255)

                                          aaTTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtgt                                                >  2:3053871/1‑93 (MQ=255)

                                             tCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCtgtgtg                                             >  1:5720245/1‑93 (MQ=255)

                                              ccAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTgtgtgt                                            <  1:1068032/93‑1 (MQ=255)

                                                ggTGCTTATGCGTCTTCACGTGCCGTTAATCTGCATCTCCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCGTCTGTGTGTTa                                          <  1:322637/92‑1 (MQ=255)

                                                aGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTa                                          >  1:481088/1‑93 (MQ=255)

                                                 gTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCGCCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTaa                                         <  1:6650313/93‑1 (MQ=255)

                                                   gCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAg                                       <  2:4845384/93‑1 (MQ=255)

                                                    cTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGt                                      <  1:1768254/93‑1 (MQ=255)

                                                     ttAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                     >  1:6227778/1‑93 (MQ=255)

                                                     ttAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                     >  2:4362301/1‑93 (MQ=255)

                                                     ttAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTa                                     >  1:576528/1‑93 (MQ=255)

                                                      tAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAg                                    <  2:5113704/93‑1 (MQ=255)

                                                      tAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAg                                    <  1:3682541/93‑1 (MQ=255)

                                                      tAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAg                                    >  1:1501923/1‑93 (MQ=255)

                                                       aaGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGGAGc                                   >  1:1947297/1‑93 (MQ=255)

                                                        aGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCg                                  >  1:3517252/1‑93 (MQ=255)

                                                        aGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCg                                  <  1:5811213/93‑1 (MQ=255)

                                                        aGCGTCTTAACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCg                                  >  1:6554565/1‑93 (MQ=255)

                                                         gCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGa                                 >  2:5446518/1‑93 (MQ=255)

                                                         gCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGa                                 >  1:4925534/1‑93 (MQ=255)

                                                         gCGTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGa                                 >  2:2307107/1‑93 (MQ=255)

                                                           gTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGaaa                               <  2:3290371/93‑1 (MQ=255)

                                                           gTCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGaaa                               >  1:4554313/1‑93 (MQ=255)

                                                            tCTTCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaag                              <  2:1349490/93‑1 (MQ=255)

                                                             tttCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaaga                             <  2:5476812/92‑1 (MQ=255)

                                                              ttCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaagaa                            <  2:2649941/93‑1 (MQ=255)

                                                              ttCACGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAaagaa                            <  2:2116077/93‑1 (MQ=255)

                                                                 aCGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGcc                         <  2:816847/93‑1 (MQ=255)

                                                                  cGTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCt                        <  2:4108674/93‑1 (MQ=255)

                                                                   gTACCGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTg                       <  1:5398377/93‑1 (MQ=255)

                                                                      ccGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTcc                    <  2:1065607/93‑1 (MQ=255)

                                                                      ccGTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTcc                    >  2:2361700/1‑93 (MQ=255)

                                                                        gTTAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGt                  <  1:861455/93‑1 (MQ=255)

                                                                         ttAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGtt                 <  2:6306296/93‑1 (MQ=255)

                                                                          tAATCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTa                <  2:6694951/93‑1 (MQ=255)

                                                                           aaTCTGCATCACCTGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAg               <  1:4432396/93‑1 (MQ=255)

                                                                           aaTCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAg               <  2:1361936/93‑1 (MQ=255)

                                                                           aaTCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAg               <  2:474879/93‑1 (MQ=255)

                                                                           aaTCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAg               >  1:2458017/1‑93 (MQ=255)

                                                                             tCTGCATCACCGGTCGCCCCGGAGAGCAGCATGGTGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAggg             >  2:4568967/1‑93 (MQ=255)

                                                                                gCATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTg          <  1:1648348/93‑1 (MQ=255)

                                                                                 cATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTgg         >  2:3700522/1‑93 (MQ=255)

                                                                                 cATCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTgg         <  1:4045784/93‑1 (MQ=255)

                                                                                  aTCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGc        >  2:3927099/1‑93 (MQ=255)

                                                                                  aTCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGc        <  1:4292269/93‑1 (MQ=255)

                                                                                  aTCACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGGTAGGGCTGGc        >  2:1326674/1‑93 (MQ=255)

                                                                                    cACCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCAc      <  1:2851844/93‑1 (MQ=255)

                                                                                     aCCGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCAcc     >  1:5955947/1‑93 (MQ=255)

                                                                                      ccGGTCGCCCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCg    >  1:4347656/1‑93 (MQ=255)

TGCTATCTCTAACTCTGGTGAATCCAGCGAAATCACGGCCTTAATTCCAGTGCTTAAGCGTCTTCACGTACCGTTAATCTGCATCACCGGTCG‑CCCGGAGAGCAGCATGGCGCGCGCCGCAGATGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGAC  >  NC_000913/3339209‑3339388

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