BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,491,023	A→G	Y318H (TAT→CAT)	idnD ←	L‑idonate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,491,023	0	A	G	93.6%	1173.2 / 52.2	420	Y318H (TAT→CAT)	idnD	L‑idonate dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (16/11); new base G (154/239); total (170/250)
Fisher's exact test for biased strand distribution p-value = 4.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

Only 100 of 447 total aligned reads displayed.

TTATTTATCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATATTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAAGAGC  >  NC_000913/4490932‑4491112

 tatttatCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGtt                                                                                         >  2:1806344/1‑93 (MQ=255)

 tatttatCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGtt                                                                                         <  1:9227351/93‑1 (MQ=255)

  atttatCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTc                                                                                        <  2:4310694/93‑1 (MQ=255)

   tttatCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCa                                                                                       >  1:6281147/1‑93 (MQ=255)

     tatCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGc                                                                                     >  2:4713567/1‑93 (MQ=255)

      atCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCa                                                                                    <  1:4606584/93‑1 (MQ=255)

       tctcTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCAc                                                                                   >  2:1174811/1‑93 (MQ=255)

       tctcTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCAc                                                                                   <  1:2938501/93‑1 (MQ=255)

          ctTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCTCTCa                                                                                <  2:159914/93‑1 (MQ=255)

          ctTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCa                                                                                <  2:6720227/93‑1 (MQ=255)

           ttAGAAAGCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAg                                                                               >  2:9013585/1‑93 (MQ=255)

           ttAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAg                                                                               >  1:9832701/1‑93 (MQ=255)

             aGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTa                                                                             >  1:6238462/1‑93 (MQ=255)

             aGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTa                                                                             <  1:7339429/93‑1 (MQ=255)

              gAAAACAAGCTGGACTTTTGCTGCCCGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTaa                                                                            >  1:7554428/1‑93 (MQ=255)

               aaaaCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTaaa                                                                           <  2:9579543/93‑1 (MQ=255)

               aaaaCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTaaa                                                                           <  2:8753729/93‑1 (MQ=255)

               aaaaCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTaaa                                                                           <  1:7047748/93‑1 (MQ=255)

               aaaaCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCCGGCAAACCTCGCGCCTTTTCCACGTCAGCGAGGGGATGCTCCGCACTCAGTaac                                                                           >  1:3449469/1‑92 (MQ=255)

                aaaCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAg                                                                          >  2:511968/1‑93 (MQ=255)

                 aaCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAgg                                                                         >  2:4592359/1‑93 (MQ=255)

                 aaCAAGCTGGACTTTTGCTGCCTGGGGTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGAGGTTCAGCACTCAGTaaacg                                                                         >  1:1355197/1‑91 (MQ=255)

                  aCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGc                                                                        <  1:8900440/93‑1 (MQ=255)

                  aCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATATTCAGCACTCAGTAAAGGc                                                                        >  2:439852/1‑93 (MQ=255)

                   cAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCa                                                                       <  1:10002172/93‑1 (MQ=255)

                    aaGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAg                                                                      >  1:5135302/1‑93 (MQ=255)

                     aGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATATTCAGCACTCAGTAAAGGCAgg                                                                     >  2:8113121/1‑92 (MQ=255)

                     aGCTGGACTTTTGCTGCCTGGGATTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCCGTGAAGGGTTGTTCAGCACTCAGTAAAGGCAgg                                                                     >  2:659389/1‑92 (MQ=255)

                      acTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAgag                                                                    <  1:8531948/92‑3 (MQ=255)

                       cgggACTTTTGCTGCCTGGGTTTTACACCCGGCGAACCGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTgg                                                                   <  1:865686/91‑1 (MQ=255)

                        tGGACTTTTGCTGCCTGGTTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGa                                                                  <  2:72075/93‑1 (MQ=255)

                        tGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGa                                                                  >  1:1222331/1‑93 (MQ=255)

                        tGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGa                                                                  <  2:2549096/93‑1 (MQ=255)

                         ggACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGAt                                                                 <  2:7753396/93‑1 (MQ=255)

                          gACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGAtt                                                                <  2:2206975/93‑1 (MQ=255)

                            cgtgtgctgcCTGGGTTTTAGCACCGGCGAAAGGTAGCGCCTCGTCCAGGTCAGTGAAGGGATATTCAGCACACAGTAAAGGAAGTGGATTGa                                                              <  2:7058320/89‑1 (MQ=255)

                             ttttGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGAt                                                             <  1:947034/93‑1 (MQ=255)

                              tttGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATa                                                            <  1:3905907/93‑1 (MQ=255)

                               ttGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATaa                                                           <  2:5956560/93‑1 (MQ=255)

                               ttGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATATTCAGCACTCAGTAAAGGCAGTGGATTGATaa                                                           <  2:2119559/93‑1 (MQ=255)

                                tgctgcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAAc                                                          <  2:8371140/93‑1 (MQ=255)

                                  ctgcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGc                                                        <  1:2266851/93‑1 (MQ=255)

                                    gcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCa                                                      <  2:3365356/93‑1 (MQ=255)

                                     ccTGGGTTTTATCACCGGCGAAACGTAGCTCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCAt                                                     <  1:4822449/93‑1 (MQ=255)

                                     ccTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCAt                                                     <  1:4186767/93‑1 (MQ=255)

                                        gggTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCg                                                  >  2:2944199/1‑93 (MQ=255)

                                         ggTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgc                                                 >  1:824067/1‑93 (MQ=255)

                                         ggTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgc                                                 <  2:4910764/93‑1 (MQ=255)

                                           ttttATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgcca                                               >  1:3749598/1‑93 (MQ=255)

                                           ttttATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgcca                                               <  1:8685217/93‑1 (MQ=255)

                                           ttttATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgcca                                               <  1:8774948/93‑1 (MQ=255)

                                           ttttATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgcca                                               <  1:1794370/93‑1 (MQ=255)

                                            tttATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgccag                                              <  1:3140248/93‑1 (MQ=255)

                                            tttATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgccag                                              >  1:4260683/1‑93 (MQ=255)

                                              tATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgccagcc                                            >  2:2739471/1‑93 (MQ=255)

                                              tATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgccagcc                                            <  1:3643042/93‑1 (MQ=255)

                                                tCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTAATAACGCCATTCGCCAGCCAt                                          <  1:4458001/93‑1 (MQ=255)

                                                 cACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGGGAAGGGATGTTAAGCACTCAGTAAAGGCAGTGGATTGATACCGCCATTCGCCAGCCATg                                         <  2:1714791/93‑1 (MQ=255)

                                                  cccGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCCCGCAGTAAAGGCGGTGGATTGATAACGCCATTCGCCAGCCATGa                                        <  2:9091205/92‑1 (MQ=255)

                                                  cccGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGa                                        <  1:8053892/92‑1 (MQ=255)

                                                  aCCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGa                                        >  2:3910611/1‑93 (MQ=255)

                                                   ccGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGac                                       <  2:9302572/93‑1 (MQ=255)

                                                    cGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGaca                                      <  2:1425656/93‑1 (MQ=255)

                                                    cGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATATTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGaca                                      <  1:7877193/93‑1 (MQ=255)

                                                     ggCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGacac                                     <  1:7971340/93‑1 (MQ=255)

                                                     ggCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGacac                                     <  1:6425089/93‑1 (MQ=255)

                                                     ggCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGacac                                     >  2:8202335/1‑93 (MQ=255)

                                                      gCGAAACCTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACt                                    <  2:376677/93‑1 (MQ=255)

                                                         aaaCGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCg                                 >  1:5598612/1‑93 (MQ=255)

                                                            cGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTa                              <  2:6034220/93‑1 (MQ=255)

                                                            cGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTa                              <  2:3478776/93‑1 (MQ=255)

                                                            cGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTa                              <  1:6662476/93‑1 (MQ=255)

                                                             gTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTAt                             <  1:9756166/93‑1 (MQ=255)

                                                             gTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTAt                             >  1:1315409/1‑93 (MQ=255)

                                                              gaGCTCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGAATTGATAACGCCATTCGCCAGCCATGACACTGCGGTAtt                            <  1:6786342/92‑1 (MQ=255)

                                                                 ctccTCTTCCAGGTCATTGAAGGGAGGTTCAGCGCTCAGTAAAGGCATTAGATTGATACCGCCATTCGCCAGCCATGACACTGCGGTATTaaa                         <  2:6327017/91‑1 (MQ=255)

                                                                   ccTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAAtt                       <  2:5578677/93‑1 (MQ=255)

                                                                   ccTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAAtt                       >  1:1902765/1‑93 (MQ=255)

                                                                   ccTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTAAGTAAAGGCAGTGGATTGATACCGCCATTCGCCAGCCATGACACTGCGGTATTAAAtt                       <  1:2725469/93‑1 (MQ=255)

                                                                   ccTCTTCCAGGTCAGGGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAAtt                       >  1:5375943/1‑93 (MQ=255)

                                                                    ctctTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTc                      >  1:8383412/1‑93 (MQ=255)

                                                                          cAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGt                <  1:1088825/93‑1 (MQ=255)

                                                                          cAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGt                <  1:5248348/93‑1 (MQ=255)

                                                                           aGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTa               <  2:6769125/93‑1 (MQ=255)

                                                                           aGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTa               <  1:4016308/93‑1 (MQ=255)

                                                                            ggTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTAGCTGGTaa              >  1:5475734/1‑93 (MQ=255)

                                                                              tCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTTCGGTATTAAATTCGCTGGTaaaa            >  2:6910521/1‑93 (MQ=255)

                                                                              tCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTaaaa            <  2:4550053/93‑1 (MQ=255)

                                                                              tCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTaaaa            >  1:5964091/1‑93 (MQ=255)

                                                                               cAGGGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAAc           >  2:2348794/1‑93 (MQ=255)

                                                                                 gTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACga         >  2:8693195/1‑92 (MQ=255)

                                                                                  tGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGa        >  2:7383721/1‑93 (MQ=255)

                                                                                   gAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGaa       <  1:5680822/93‑1 (MQ=255)

                                                                                    aaGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGaaa      <  2:5761101/93‑1 (MQ=255)

                                                                                      gggATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAaga    >  1:4988962/1‑93 (MQ=255)

                                                                                       ggATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAagag   >  2:5451022/1‑93 (MQ=255)

                                                                                       ggATGTTAAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAagag   <  2:8574694/93‑1 (MQ=255)

                                                                                        gATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAAGAGc  <  1:9293004/93‑1 (MQ=255)

                                                                                        gATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAAGAGc  <  1:2475655/93‑1 (MQ=255)

                                                                                        gATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAAGAGc  <  1:213952/93‑1 (MQ=255)

TTATTTATCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATATTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAAGAGC  >  NC_000913/4490932‑4491112

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: `‑`