BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,194,729	A→G	G83G (GGA→GGG)	yjaD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,194,729	0	A	G	86.5%	481.4 / 61.7	185	G83G (GGA→GGG)	yjaD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (16/9); new base G (90/70); total (106/79)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 192 total aligned reads displayed.

TGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGAT  >  NC_000913/4194639‑4194821

tGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCg                                                                                            <  1:4303625/93‑1 (MQ=255)

tGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACg                                                                                            >  2:5726609/1‑93 (MQ=255)

tGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGAGGTTGGGCTGTTTCAACTGGACGGGCg                                                                                            >  1:1034535/1‑93 (MQ=255)

  tttGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAg                                                                                          <  2:4492921/93‑1 (MQ=255)

   ttGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAgg                                                                                         <  1:5190126/93‑1 (MQ=255)

    tGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGc                                                                                        <  2:573986/93‑1 (MQ=255)

    tGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGc                                                                                        >  2:2011511/1‑93 (MQ=255)

     ggTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCg                                                                                       >  1:4644250/1‑93 (MQ=255)

       ttAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTa                                                                                     <  1:3082172/93‑1 (MQ=255)

        tAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTAc                                                                                    <  1:476643/93‑1 (MQ=255)

         aGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACa                                                                                   >  2:5213287/1‑93 (MQ=255)

          gTACAACAGCAGCGGCGTCACGATATGGTGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACaa                                                                                  <  1:6564776/93‑1 (MQ=255)

            acaacaGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACt                                                                                <  2:2895676/93‑1 (MQ=255)

            acaacaGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACt                                                                                <  2:1657026/93‑1 (MQ=255)

            acaacaGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACt                                                                                <  1:5447445/93‑1 (MQ=255)

             caacaGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTg                                                                               <  2:3424739/93‑1 (MQ=255)

               acaGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGc                                                                             <  1:2121080/93‑1 (MQ=255)

                cagcagCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCg                                                                            >  2:2870591/1‑93 (MQ=255)

                cagcagCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCg                                                                            >  2:2967781/1‑93 (MQ=255)

                  gcagcGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGa                                                                          >  2:802218/1‑93 (MQ=255)

                  gcagcGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGa                                                                          >  2:6077433/1‑93 (MQ=255)

                   cagcGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAg                                                                         <  1:99798/93‑1 (MQ=255)

                     gcggcgTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAggt                                                                       >  1:5043839/1‑91 (MQ=255)

                     gcggcgTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGtt                                                                       <  1:1247168/93‑1 (MQ=255)

                     gcggcgTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGtt                                                                       >  2:6343877/1‑93 (MQ=255)

                      gggcgTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGttt                                                                      <  2:5669000/92‑1 (MQ=255)

                       ggcgTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGtttt                                                                     >  2:3101389/1‑93 (MQ=255)

                        gcgTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTa                                                                    >  1:1560333/1‑93 (MQ=255)

                          gTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTAcc                                                                  >  2:3606961/1‑93 (MQ=255)

                          gTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTAcc                                                                  <  2:6262647/93‑1 (MQ=255)

                          gTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTAcc                                                                  >  2:194704/1‑93 (MQ=255)

                          gTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGGTTTAcc                                                                  >  2:4616572/1‑93 (MQ=255)

                          gTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGCTTTAcc                                                                  >  1:2333660/1‑93 (MQ=255)

                           tCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCg                                                                 >  2:3513645/1‑93 (MQ=255)

                           tCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCg                                                                 >  1:4507650/1‑93 (MQ=255)

                            cACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACACCTGGCGGAGTTTTCCCGa                                                                >  1:3339119/1‑93 (MQ=255)

                               gATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCg                                                             >  2:4791843/1‑93 (MQ=255)

                                 tatGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCa                                                           <  2:2505730/93‑1 (MQ=255)

                                  atGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCAt                                                          >  1:3288273/1‑93 (MQ=255)

                                  atGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCAt                                                          >  1:4336579/1‑93 (MQ=255)

                                   tGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGCCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATa                                                         <  2:6511587/93‑1 (MQ=255)

                                    ggggTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGGTTTACCGATCGCAtca                                                        >  1:5137278/1‑91 (MQ=255)

                                     gggTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATaaa                                                       <  1:6066076/93‑1 (MQ=255)

                                      ggTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAAt                                                      >  2:5938462/1‑93 (MQ=255)

                                      ggTCAGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAAt                                                      <  1:2685718/93‑1 (MQ=255)

                                       gTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATa                                                     <  1:6409500/93‑1 (MQ=255)

                                       gTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATa                                                     >  2:1904901/1‑93 (MQ=255)

                                        tCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATAc                                                    <  1:136954/93‑1 (MQ=255)

                                           gTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACtgt                                                 >  2:2367198/1‑93 (MQ=255)

                                            tACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACtgtg                                                >  1:1011008/1‑93 (MQ=255)

                                            tACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGGGCAAAACTGGCGGCGGGTTACCGATTGCGTCAATACtgtg                                                >  1:716107/1‑93 (MQ=255)

                                             aCGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTgg                                               <  2:659402/93‑1 (MQ=255)

                                               gTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGtt                                             >  2:320219/1‑93 (MQ=255)

                                                 cAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTAc                                           >  2:3416087/1‑93 (MQ=255)

                                                  aGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACt                                          >  1:2854931/1‑93 (MQ=255)

                                                   ggTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTg                                         >  2:4537920/1‑93 (MQ=255)

                                                   ggTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTg                                         >  2:4430083/1‑93 (MQ=255)

                                                    gTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGc                                        >  1:6272609/1‑93 (MQ=255)

                                                    gTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGc                                        <  1:5601295/93‑1 (MQ=255)

                                                    gTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGc                                        <  1:2333732/93‑1 (MQ=255)

                                                      cATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCgg                                      <  2:5036988/93‑1 (MQ=255)

                                                       aTTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCggg                                     >  2:600872/1‑93 (MQ=255)

                                                        ttGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGc                                    >  2:5382171/1‑93 (MQ=255)

                                                        ttGATCTCGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGc                                    >  2:1463544/1‑93 (MQ=255)

                                                            tctcGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGa                                >  2:1025062/1‑93 (MQ=255)

                                                            tctcGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGa                                >  1:1461320/1‑93 (MQ=255)

                                                              tcGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGaaa                              >  1:2465648/1‑93 (MQ=255)

                                                              tcGATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGaaa                              <  2:6222575/93‑1 (MQ=255)

                                                              tcGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGaaa                              <  2:5751480/93‑1 (MQ=255)

                                                                gATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATg                            >  1:4292990/1‑93 (MQ=255)

                                                                gATGTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATg                            >  1:5719822/1‑93 (MQ=255)

                                                                 aTGTTGGGCTGTCTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGt                           <  1:1843994/93‑1 (MQ=255)

                                                                   gTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTAt                         <  1:5613213/93‑1 (MQ=255)

                                                                   gTTGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTAt                         <  2:4485622/93‑1 (MQ=255)

                                                                    ttGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATc                        <  1:5248718/93‑1 (MQ=255)

                                                                     tGGGCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATcc                       <  1:2068159/93‑1 (MQ=255)

                                                                       ggCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGa                     <  1:1988348/93‑1 (MQ=255)

                                                                        gCTGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAg                    >  2:5176497/1‑93 (MQ=255)

                                                                          tGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCa                  >  1:1343509/1‑93 (MQ=255)

                                                                          tGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCa                  <  1:3884503/93‑1 (MQ=255)

                                                                          tGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCa                  >  2:4099481/1‑93 (MQ=255)

                                                                          tGTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCa                  >  1:4021354/1‑93 (MQ=255)

                                                                          tGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCa                  >  2:3855067/1‑93 (MQ=255)

                                                                           gTTTCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCaa                 <  1:1901425/93‑1 (MQ=255)

                                                                            tttCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCaaa                <  2:4711960/93‑1 (MQ=255)

                                                                            tctCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCaaa                >  1:4248220/3‑93 (MQ=255)

                                                                             ttCAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCaaaa               >  2:5168897/1‑93 (MQ=255)

                                                                               cAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAAcc             <  1:6426978/93‑1 (MQ=255)

                                                                               cAACTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAAcc             >  1:5774215/1‑93 (MQ=255)

                                                                                 aCTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGa           >  2:21607/1‑93 (MQ=255)

                                                                                  cTGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGaa          >  2:5875843/1‑93 (MQ=255)

                                                                                   tGGCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAAt         <  1:4359862/93‑1 (MQ=255)

                                                                                    gaccGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATg        <  1:2393019/91‑1 (MQ=255)

                                                                                     gCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATgg       <  2:1569450/93‑1 (MQ=255)

                                                                                     gCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATgg       >  1:5870914/1‑93 (MQ=255)

                                                                                     gCCGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATgg       <  1:4432772/93‑1 (MQ=255)

                                                                                      ccGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATggg      >  2:365321/1‑93 (MQ=255)

                                                                                      ccGGGCGAGGCGTACAACTAGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATggg      >  1:6495578/1‑93 (MQ=255)

                                                                                       cGGGCGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGc     >  2:2408631/1‑93 (MQ=255)

                                                                                         gACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGa   >  1:2980885/1‑93 (MQ=255)

TGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGAT  >  NC_000913/4194639‑4194821

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