BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,056,250	(T)_7→6	coding (1619/1644 nt)	yeeO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,056,244	0	T	.	99.5%	917.6 / ‑3.5	200	coding (1625/1644 nt)	yeeO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base . (112/87); total (112/88)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

TGCCCACCCACAATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGCTTTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCACGGCCCGCAC  >  NC_000913/2056153‑2056332

tGCCCACCCACAATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGCtt                                                                                         <  1:5327991/93‑1 (MQ=255)

  cccacccaCAATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGCtttt                                                                                       <  1:3837018/93‑1 (MQ=255)

  cccacccaCAATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGCtttt                                                                                       <  2:1262260/93‑1 (MQ=255)

  cccacccaCAATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGCtttt                                                                                       <  2:4552513/93‑1 (MQ=255)

    cacccacAATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGCtttttt                                                                                     >  1:5988174/1‑93 (MQ=255)

       ccacAATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGCttttttcac                                                                                  >  1:5129781/1‑90 (MQ=255)

         acaATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCacac                                                                               >  2:319615/1‑93 (MQ=255)

         acaATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCacac                                                                               <  2:3942874/93‑1 (MQ=255)

            aTGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACttt                                                                            <  2:2489992/93‑1 (MQ=255)

            aTGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACttt                                                                            <  2:4949319/93‑1 (MQ=255)

             tGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACtttt                                                                           >  2:598483/1‑93 (MQ=255)

              gTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTg                                                                          <  1:117706/93‑1 (MQ=255)

              gTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTg                                                                          <  1:1608572/93‑1 (MQ=255)

               tCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGc                                                                         <  1:757843/93‑1 (MQ=255)

               tCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGc                                                                         <  1:4564474/93‑1 (MQ=255)

               tCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGc                                                                         <  2:5973081/93‑1 (MQ=255)

                 aaaGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCgg                                                                       >  2:3273099/1‑93 (MQ=255)

                    gCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTc                                                                    <  1:2583012/93‑1 (MQ=255)

                     ccTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCg                                                                   <  2:5260783/93‑1 (MQ=255)

                     ccTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCg                                                                   <  1:3490292/93‑1 (MQ=255)

                      cTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCgg                                                                  <  2:4641882/93‑1 (MQ=255)

                       tGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCggg                                                                 <  2:2390858/93‑1 (MQ=255)

                        gTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGt                                                                <  1:3937829/93‑1 (MQ=255)

                         ttttttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTc                                                               >  2:1918685/1‑93 (MQ=255)

                          tttttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCg                                                              >  1:5899324/1‑93 (MQ=255)

                           ttttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGa                                                             <  2:2332561/93‑1 (MQ=255)

                           ttttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGa                                                             <  2:5048314/93‑1 (MQ=255)

                           ttttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGa                                                             >  2:887690/1‑93 (MQ=255)

                           ttttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGa                                                             <  1:1520398/93‑1 (MQ=255)

                           ttttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGGCGa                                                             >  1:2078689/1‑93 (MQ=255)

                            tttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAg                                                            <  1:4698768/93‑1 (MQ=255)

                            tttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAg                                                            >  1:3099324/1‑93 (MQ=255)

                             ttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAgg                                                           <  1:4041617/93‑1 (MQ=255)

                             ttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAgg                                                           <  1:484118/93‑1 (MQ=255)

                             ttAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAgg                                                           <  2:5539390/93‑1 (MQ=255)

                              tAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAggg                                                          <  1:3371809/93‑1 (MQ=255)

                              tAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAggg                                                          <  1:546130/93‑1 (MQ=255)

                               aaaTTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGt                                                         <  2:130687/93‑1 (MQ=255)

                               aaaTTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGt                                                         >  1:2260052/1‑93 (MQ=255)

                               aaaTTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGt                                                         <  1:952695/93‑1 (MQ=255)

                                aaTTTGAAATCGTTTGCTGTAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAgggga                                                        >  2:6125447/1‑91 (MQ=255)

                                aaTTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGTC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTa                                                        <  1:836709/93‑1 (MQ=255)

                                     gAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTcc                                                   <  1:2422300/93‑1 (MQ=255)

                                      aaaTCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTTGAGGGTATTTCCa                                                  >  1:5282742/1‑93 (MQ=255)

                                      aaaTCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTTCa                                                  >  2:3842825/1‑93 (MQ=255)

                                        aTCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATa                                                >  1:5022591/1‑93 (MQ=255)

                                          cGTTTGCTGGAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGc                                              >  2:2150145/1‑93 (MQ=255)

                                           gTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCCGGGTGAGGGTACTTCCATAGcc                                             >  2:241884/1‑93 (MQ=255)

                                             ttGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCaa                                           <  2:2819846/93‑1 (MQ=255)

                                             ttGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCaa                                           >  2:2951650/1‑93 (MQ=255)

                                              tGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTAtt                                                      >  2:5598812/1‑81 (MQ=255)

                                              tGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTAtt                                                      <  1:5598812/81‑1 (MQ=255)

                                              tGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAAc                                          <  1:6069795/93‑1 (MQ=255)

                                              tGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAAc                                          <  1:6117606/93‑1 (MQ=255)

                                              tGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGG‑TTTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAAc                                          <  2:3555364/93‑1 (MQ=255)

                                               gCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACg                                         <  2:3213601/93‑1 (MQ=255)

                                                 tGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTc                                       >  1:4907199/1‑93 (MQ=255)

                                                 tGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTc                                       >  1:6182317/1‑93 (MQ=255)

                                                 tGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTc                                       <  1:2785267/93‑1 (MQ=255)

                                                 tGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCAACGTc                                       >  1:509329/1‑93 (MQ=255)

                                                  gAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTcc                                      >  2:1546236/1‑93 (MQ=255)

                                                  gAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTcc                                      <  1:328405/93‑1 (MQ=255)

                                                   aaaTAATCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCAACGTCCa                                     >  2:2624936/1‑93 (MQ=255)

                                                     ataatCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAga                                   >  1:3244180/1‑92 (MQ=255)

                                                     ataatCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCCACGTGCAGt                                   >  1:3449893/1‑93 (MQ=255)

                                                     ataatCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGCCTTTCCCTAGCCAACGTCCAga                                   >  1:2951027/1‑92 (MQ=255)

                                                       aatCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCCACGTCCAGTaa                                 >  2:5772047/1‑93 (MQ=255)

                                                       aatCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTaa                                 <  2:2254240/93‑1 (MQ=255)

                                                         tCTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAAcc                               <  2:692179/93‑1 (MQ=255)

                                                          cTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCa                              <  1:1268213/93‑1 (MQ=255)

                                                          cTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCa                              <  1:5177646/93‑1 (MQ=255)

                                                          cTGCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCa                              <  2:3799359/93‑1 (MQ=255)

                                                            gCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACCCTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCAtt                            <  1:2526420/93‑1 (MQ=255)

                                                            gCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCAtt                            >  1:1759653/1‑93 (MQ=255)

                                                            gCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCAtt                            <  1:3749078/93‑1 (MQ=255)

                                                            gCATTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCAACGTCCAATAACCAtt                            >  2:1713959/1‑93 (MQ=255)

                                                              aTTTTGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCg                          <  2:3662415/93‑1 (MQ=255)

                                                               ttttGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCAACGTCCAGTAACCATTCGc                         >  1:2708115/1‑93 (MQ=255)

                                                               ttttGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAACCAACGTCCAGTAACCATTCGc                         >  2:3419751/1‑93 (MQ=255)

                                                               ttttGTCGTTTATCCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGc                         <  1:1042135/93‑1 (MQ=255)

                                                                tttGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGcc                        <  2:4822607/93‑1 (MQ=255)

                                                                 ttGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGccc                       <  1:5941606/93‑2 (MQ=255)

                                                                  tGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAg                      <  2:654541/93‑1 (MQ=255)

                                                                  tGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAg                      >  1:5157250/1‑93 (MQ=255)

                                                                  tGTCGTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAg                      >  1:658126/1‑93 (MQ=255)

                                                                      gTTTATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCCACGTCCAGTAACCATTCCCCAGTaaa                  >  2:4299546/1‑93 (MQ=255)

                                                                       tttATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTaaaa                 >  2:2703319/1‑93 (MQ=255)

                                                                        ttATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAAc                <  2:1402789/93‑1 (MQ=255)

                                                                         tATTCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACa               <  1:3898047/93‑1 (MQ=255)

                                                                            tCCGACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCa            >  2:4110035/1‑93 (MQ=255)

                                                                               gACACAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAACACgg         >  1:4317558/1‑93 (MQ=255)

                                                                                acacAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCACGGc        >  2:2783872/1‑93 (MQ=255)

                                                                                acacAACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCACCGc        >  1:4017862/1‑93 (MQ=255)

                                                                                 cacaACTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCACGGcc       <  1:1950998/93‑1 (MQ=255)

                                                                                  acaTCTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCACGGccc      <  2:4220497/93‑1 (MQ=255)

                                                                                   caACTGGC‑TTTTTTCACACTTTTGCGGCTCCGGTCGAGGGGATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGAACGGCCCg     >  1:2496698/1‑93 (MQ=255)

                                                                                    aaCTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCACGGCCCGc    >  2:1187541/1‑93 (MQ=255)

                                                                                     aCTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCACGGCCCGca   <  1:1172257/93‑1 (MQ=255)

                                                                                     aCTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTAGAGGGAATTTCCATAGCCAACGTCCAGCAACGCTTCTACAGTAAAACAGCACGGCCCGca   >  2:1883663/1‑93 (MQ=255)

                                                                                      cTGGC‑TTTTTTCACACTTTTGCGGCTCGGGTCGAGGGGATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCACGGCCCGcac  >  2:4382312/1‑93 (MQ=255)

TGCCCACCCACAATGTCAAAGCCTGTTTTTTAAATTTGAAATCGTTTGCTGAAATAATCTGCATTTTGTCGTTTATTCCGACACAACTGGCTTTTTTTCACACTTTTGCGGCTCGGGTCGAGGGTATTTCCATAGCCAACGTCCAGTAACCATTCGCCAGTAAAACAGCACGGCCCGCAC  >  NC_000913/2056153‑2056332

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