BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,253,742	Δ1 bp	coding (551/552 nt)	yhaNM ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,253,741	0	G	.	99.2%	586.0 / ‑3.4	131	coding (552/552 nt)	yhaNM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base . (70/60); minor base C (1/0); total (71/60)
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.97e-01

Only 100 of 137 total aligned reads displayed.

TAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGA  >  NC_000913/3253650‑3253829

  aCGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGttt                                                                                       >  1:3495351/1‑93 (MQ=255)

    gCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTg                                                                                    <  2:2173633/93‑1 (MQ=255)

    gCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTg                                                                                    <  2:1117043/93‑1 (MQ=255)

    gCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTg                                                                                    >  1:958702/1‑93 (MQ=255)

      acaATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGa                                                                                  <  2:3072403/93‑1 (MQ=255)

      acaATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGa                                                                                  >  2:3633583/1‑93 (MQ=255)

       caATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGaa                                                                                 >  2:499063/1‑93 (MQ=255)

        aaTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGaaa                                                                                <  1:2699773/93‑1 (MQ=255)

        aaTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGaaa                                                                                <  1:580734/93‑1 (MQ=255)

         aTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAg                                                                               >  1:1281020/1‑93 (MQ=255)

         aTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAg                                                                               >  1:4558273/1‑93 (MQ=255)

         aTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAg                                                                               >  1:2522754/1‑93 (MQ=255)

            tAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCAc                                                                            <  1:1868886/93‑1 (MQ=255)

             aGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCAcc                                                                           <  1:3845592/93‑1 (MQ=255)

              gCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCg                                                                          >  2:3782494/1‑93 (MQ=255)

                ttCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTc                                                                        <  2:1531125/93‑1 (MQ=255)

                ttCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTc                                                                        >  2:3386059/1‑93 (MQ=255)

                 tCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCa                                                                       <  1:877753/93‑1 (MQ=255)

                   gcgGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGc                                                                     <  1:3880876/93‑1 (MQ=255)

                   gcgGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGc                                                                     <  1:1737231/93‑1 (MQ=255)

                    cgGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACCTCAATGATGAACTACAGGAGCGCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCg                                                                    >  2:4226148/1‑93 (MQ=255)

                     ggAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCgg                                                                   >  1:577528/1‑93 (MQ=255)

                     ggAATCGAGAATAAAGCGGATCGCCGCAAACCGGACATCATTAACGCACTTCAGGATCTCAGGCAGCG‑CAGTTCTGTAAAGCACCGTCAGCgg                                                                   >  2:3228521/1‑93 (MQ=255)

                      gAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGga                                                                  >  1:772479/1‑93 (MQ=255)

                      gAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGga                                                                  <  2:1241853/93‑1 (MQ=255)

                       aaTCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGgag                                                                 >  1:1696652/1‑93 (MQ=255)

                        aTCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGgaga                                                                >  1:4416117/1‑93 (MQ=255)

                         tCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGAc                                                               <  2:2899814/93‑1 (MQ=255)

                          cGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGAct                                                              >  2:832439/1‑93 (MQ=255)

                           gagaATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGActc                                                             >  1:4481686/1‑93 (MQ=255)

                           gagaATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGActc                                                             >  2:569643/1‑93 (MQ=255)

                           gagaATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGActc                                                             >  2:758767/1‑93 (MQ=255)

                             gaATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCtt                                                           <  1:383979/93‑1 (MQ=255)

                              aaTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTg                                                          >  2:1794889/1‑93 (MQ=255)

                                 aaaGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGGAAGCACCGTCAGCGGAGACTCTTGCTc                                                       >  1:324348/1‑93 (MQ=255)

                                 aaaGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTc                                                       >  1:324347/1‑93 (MQ=255)

                                   aGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGc                                                     >  2:1640821/1‑93 (MQ=255)

                                    gCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGcc                                                    <  2:3181693/93‑1 (MQ=255)

                                     cGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGccc                                                   >  2:720951/1‑93 (MQ=255)

                                      ggATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGTAGACTCTTGCTCGcccc                                                  <  1:693143/93‑2 (MQ=255)

                                       gATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGccccc                                                 <  2:2725193/93‑3 (MQ=255)

                                       gATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGccccc                                                 <  1:3142987/93‑3 (MQ=255)

                                         tCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTg                                               >  2:1296031/1‑93 (MQ=255)

                                         tCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTg                                               <  1:3992164/93‑1 (MQ=255)

                                         tCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTg                                               <  2:4250559/93‑1 (MQ=255)

                                         tCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTg                                               <  1:2766702/93‑1 (MQ=255)

                                           gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                             >  1:140221/1‑93 (MQ=255)

                                           gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                             >  1:1548490/1‑93 (MQ=255)

                                           gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                             <  2:2054139/93‑1 (MQ=255)

                                            ccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCca                                            >  1:3967410/1‑93 (MQ=255)

                                             cgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcac                                           >  2:2683886/1‑93 (MQ=255)

                                              gcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcaca                                          >  1:2891367/1‑93 (MQ=255)

                                              gcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcaca                                          >  1:410734/1‑93 (MQ=255)

                                               cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacac                                         >  1:2611276/1‑93 (MQ=255)

                                               cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacac                                         <  2:2369874/93‑1 (MQ=255)

                                               cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTAAGCGGAGACTCTTGCTCGCCCCATGCcacac                                         >  2:1508165/1‑93 (MQ=255)

                                                aaaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacaca                                        >  1:3333272/1‑93 (MQ=255)

                                                 aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCacacac                                       >  1:2392508/1‑93 (MQ=255)

                                                   cGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGc                                     >  2:2786429/1‑93 (MQ=255)

                                                    gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGcc                                    <  1:4747675/93‑1 (MQ=255)

                                                     ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                   >  2:4392421/1‑93 (MQ=255)

                                                     ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                   <  1:1181333/93‑1 (MQ=255)

                                                     ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                   <  1:848565/93‑1 (MQ=255)

                                                     ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                   >  1:1310616/1‑93 (MQ=255)

                                                     ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                   >  1:4349286/1‑93 (MQ=255)

                                                      gACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctg                                  <  2:838210/93‑1 (MQ=255)

                                                       gcTTTATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgc                                 <  1:2986017/92‑1 (MQ=255)

                                                       aCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACATGCctgc                                 <  2:4608636/93‑1 (MQ=255)

                                                       aCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgc                                 <  2:1874687/93‑1 (MQ=255)

                                                        cTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgct                                >  2:4554551/1‑93 (MQ=255)

                                                        cTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgct                                >  2:3758484/1‑93 (MQ=255)

                                                         ttCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgctg                               >  1:1442137/1‑93 (MQ=255)

                                                         ttCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgctg                               <  2:2934718/93‑1 (MQ=255)

                                                         ttCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgctg                               <  2:2699665/93‑1 (MQ=255)

                                                          tCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTgg                              >  1:3416340/1‑93 (MQ=255)

                                                               gACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTaaa                         >  2:2609197/1‑93 (MQ=255)

                                                               gACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTaaa                         <  1:2481780/93‑1 (MQ=255)

                                                                 cGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAaca                       <  2:4472898/93‑1 (MQ=255)

                                                                  gAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacac                      <  2:3413814/93‑1 (MQ=255)

                                                                    aCTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacacca                    >  1:2773904/1‑93 (MQ=255)

                                                                     cTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacaccac                   <  1:291255/93‑1 (MQ=255)

                                                                         aGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAAcaccacacca               >  1:3024759/1‑93 (MQ=255)

                                                                           gATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATc             <  2:1065774/93‑1 (MQ=255)

                                                                           gATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATc             >  1:638689/1‑93 (MQ=255)

                                                                            aTCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg            >  1:2941810/1‑93 (MQ=255)

                                                                            aTCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg            >  1:10681/1‑93 (MQ=255)

                                                                               tcAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGc         >  2:4113297/1‑93 (MQ=255)

                                                                               tcAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGc         >  1:2081915/1‑93 (MQ=255)

                                                                               tcAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGc         >  2:225586/1‑93 (MQ=255)

                                                                               ccagccagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGc         <  2:1432107/92‑1 (MQ=255)

                                                                                  gccagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtc      >  1:3011430/1‑93 (MQ=255)

                                                                                   ccagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtct     <  1:1912397/93‑1 (MQ=255)

                                                                                   ccagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtct     <  1:2656199/93‑1 (MQ=255)

                                                                                    cagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtctc    <  1:365736/93‑1 (MQ=255)

                                                                                    cagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtctc    <  1:3907522/93‑1 (MQ=255)

                                                                                    cagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtctc    >  1:451717/1‑93 (MQ=255)

                                                                                      gcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCgg  >  1:226909/1‑92 (MQ=255)

                                                                                      gcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa  >  1:226910/1‑93 (MQ=255)

                                                                                      gcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa  <  1:3216469/93‑1 (MQ=255)

                                                                                      gcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa  <  1:3087371/93‑1 (MQ=255)

TAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGA  >  NC_000913/3253650‑3253829

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