BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,229,837	C→T	intergenic (‑214/‑153)	hlyE ← / → umuD	hemolysin E/UmuD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,229,837	0	C	T	82.0%	446.0 / 71.0	183	intergenic (‑214/‑153)	hlyE/umuD	hemolysin E/UmuD
Reads supporting (aligned to +/- strand): ref base C (32/1); new base T (111/39); total (143/40)
Fisher's exact test for biased strand distribution p-value = 2.15e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-02

Only 100 of 191 total aligned reads displayed.

GATAAATATCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATT  >  NC_000913/1229746‑1229929

gATAAATATCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGcg                                                                                             >  1:4613078/1‑93 (MQ=255)

gATAAATATCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGcg                                                                                             >  2:6193854/1‑93 (MQ=255)

 aTAAATATCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAgtgt                                                                                            <  2:3372315/93‑1 (MQ=255)

  tAAATATCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGcgtc                                                                                           >  2:5296485/1‑93 (MQ=255)

    aaTATCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGtgtcgt                                                                                         >  2:3570014/1‑93 (MQ=255)

    aaTATCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGtgtcgt                                                                                         >  1:2822129/1‑93 (MQ=255)

    aaTATCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGcgtcgt                                                                                         >  2:2778002/1‑93 (MQ=255)

       atCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGc                                                                                      >  1:7292862/1‑93 (MQ=255)

        tCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGcc                                                                                     >  1:7236626/1‑93 (MQ=255)

        tCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGcc                                                                                     >  1:3205398/1‑93 (MQ=255)

          aaaCAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAg                                                                                   >  2:4163705/1‑93 (MQ=255)

          aaaCAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAg                                                                                   >  1:3515702/1‑93 (MQ=255)

           aaCAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGa                                                                                  >  2:1822623/1‑93 (MQ=255)

           aaCAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGa                                                                                  >  2:6733840/1‑93 (MQ=255)

           aaCAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGa                                                                                  >  1:5444806/1‑93 (MQ=255)

           aaCAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGa                                                                                  >  2:5989444/1‑93 (MQ=255)

           aaCAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGa                                                                                  >  1:2962410/1‑93 (MQ=255)

            aCAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGaa                                                                                 <  1:7119192/93‑1 (MQ=255)

            aCAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGaa                                                                                 <  2:3110571/93‑1 (MQ=255)

             cAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAg                                                                                >  1:6784904/1‑93 (MQ=255)

             cAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAg                                                                                >  2:6720340/1‑93 (MQ=255)

             cAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAg                                                                                >  2:7326399/1‑93 (MQ=255)

              aaTGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAgg                                                                               >  2:5791347/1‑93 (MQ=255)

                  tttAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCAc                                                                           >  1:6432172/1‑93 (MQ=255)

                  tttAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCAc                                                                           <  2:5491938/93‑1 (MQ=255)

                  tttAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCAc                                                                           >  2:5030389/1‑93 (MQ=255)

                     aaTGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTg                                                                        >  1:1414786/1‑93 (MQ=255)

                     aaTGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTg                                                                        >  1:6135901/1‑93 (MQ=255)

                      aTGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGa                                                                       >  1:6248677/1‑93 (MQ=255)

                       tGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAg                                                                      >  1:4244870/1‑93 (MQ=255)

                       tGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAg                                                                      >  2:5185772/1‑93 (MQ=255)

                       tGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAg                                                                      >  2:4647405/1‑93 (MQ=255)

                        gTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGc                                                                     >  2:4333504/1‑93 (MQ=255)

                          cATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGcac                                                                   >  1:1860633/1‑93 (MQ=255)

                           aTTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGcaca                                                                  >  1:3394444/1‑93 (MQ=255)

                           aTTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGcaca                                                                  >  1:6962681/1‑93 (MQ=255)

                            ttATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACaa                                                                 >  1:2456458/1‑93 (MQ=255)

                              aTGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGa                                                               >  2:1731647/1‑93 (MQ=255)

                               tGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGAt                                                              >  2:3933413/1‑93 (MQ=255)

                                ggCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATa                                                             <  2:1648131/93‑1 (MQ=255)

                                 gCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATaa                                                            >  2:562072/1‑93 (MQ=255)

                                   gAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAaga                                                          >  1:1346531/1‑93 (MQ=255)

                                   gAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAaga                                                          >  1:5021841/1‑93 (MQ=255)

                                      tGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGaa                                                       >  2:7297758/1‑93 (MQ=255)

                                       gcctctattctatTTTTTAGCCGGGGGATATTTTTCTTTTCTGCTGGATGAGCGTCGTGGCCGGAAGGCAACGGGAGCACAGGTTAAGAGaag                                                      >  2:6960429/4‑92 (MQ=255)

                                       gCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAAc                                                      <  1:1198648/93‑1 (MQ=255)

                                       gCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAAc                                                      >  2:439176/1‑93 (MQ=255)

                                       gCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAAc                                                      >  1:3856259/1‑93 (MQ=255)

                                          tctattctatTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGaa                                                   >  1:6610160/1‑93 (MQ=255)

                                           ctattctattTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGaaa                                                  >  1:6906410/1‑93 (MQ=255)

                                           ctattctattTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGaaa                                                  >  2:223722/1‑93 (MQ=255)

                                            tattctattTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGaaaa                                                 >  1:5344623/1‑93 (MQ=255)

                                              ttctattTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAAt                                               <  1:5248624/93‑1 (MQ=255)

                                               tctattTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATc                                              >  1:7459207/1‑93 (MQ=255)

                                               tctattTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATc                                              >  2:271/1‑93 (MQ=255)

                                                ctattTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATca                                             >  2:935494/1‑93 (MQ=255)

                                                ctattTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATca                                             >  2:4864825/1‑93 (MQ=255)

                                                   ttttttAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATcagca                                          <  1:1931926/93‑1 (MQ=255)

                                                   ttttttAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATcagca                                          >  1:2477427/1‑93 (MQ=255)

                                                    tttttAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATcagcag                                         <  2:2338546/93‑1 (MQ=255)

                                                      tttAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGcc                                       <  2:3233831/93‑1 (MQ=255)

                                                      tttAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGcc                                       >  1:3828961/1‑93 (MQ=255)

                                                       ttAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCt                                      >  2:2547661/1‑93 (MQ=255)

                                                       ttAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCt                                      <  2:495879/93‑1 (MQ=255)

                                                            cGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCa                                 <  2:1346032/93‑1 (MQ=255)

                                                            cGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCa                                 <  1:6107149/93‑1 (MQ=255)

                                                            cGGGTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCa                                 >  2:2821723/1‑93 (MQ=255)

                                                              ggTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGc                               >  1:8014/1‑93 (MQ=255)

                                                              ggTGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGc                               >  1:2178408/1‑93 (MQ=255)

                                                                tGATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGa                             >  2:1543738/1‑93 (MQ=255)

                                                                 gATATTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGAc                            >  2:5096991/1‑93 (MQ=255)

                                                                  atatTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACa                           >  1:1159711/1‑93 (MQ=255)

                                                                   tatTTTTCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACaa                          >  1:2132654/1‑93 (MQ=255)

                                                                     tttttCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTTTGCAGCGACAAat                        >  1:4900312/1‑93 (MQ=255)

                                                                     tttttCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAat                        >  2:6138636/1‑93 (MQ=255)

                                                                      ttttCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAata                       <  1:3665841/93‑1 (MQ=255)

                                                                         tCATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTg                    >  1:5081487/1‑93 (MQ=255)

                                                                          cATTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGa                   <  1:4976571/93‑1 (MQ=255)

                                                                           aTTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGAt                  <  2:20481/93‑1 (MQ=255)

                                                                           aTTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGAt                  >  1:4428300/1‑93 (MQ=255)

                                                                           aTTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGAt                  >  1:7078805/1‑93 (MQ=255)

                                                                           aTTTCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGAt                  >  2:4626504/1‑93 (MQ=255)

                                                                            tttCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATa                 >  1:1066521/1‑93 (MQ=255)

                                                                            tttCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATa                 <  2:6872135/93‑1 (MQ=255)

                                                                            tttCCGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATa                 >  2:5720321/1‑93 (MQ=255)

                                                                              tCTGCTGGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGc               <  2:7046782/93‑1 (MQ=255)

                                                                                tgctgGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCt             >  2:3317555/1‑93 (MQ=255)

                                                                                 gctgGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTg            <  1:1563594/93‑1 (MQ=255)

                                                                                 gctgGATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTg            >  2:6609031/1‑93 (MQ=255)

                                                                                 gctgGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTg            >  2:1922970/1‑93 (MQ=255)

                                                                                   tgGATGAGTGTCGTTGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGaa          >  2:616196/1‑93 (MQ=255)

                                                                                     gATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATc        >  2:3703513/1‑93 (MQ=255)

                                                                                     gATGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATc        <  1:2062777/93‑1 (MQ=255)

                                                                                      aTGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCa       <  2:1190526/93‑1 (MQ=255)

                                                                                      aTGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCa       <  1:2195983/93‑1 (MQ=255)

                                                                                      aTGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCa       >  1:19422/1‑93 (MQ=255)

                                                                                       tGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAg      >  1:1034039/1‑93 (MQ=255)

                                                                                       tGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAg      <  1:767434/93‑1 (MQ=255)

                                                                                       tGAGTGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAg      >  1:5630115/1‑93 (MQ=255)

                                                                                       tGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAg      >  1:5458963/1‑93 (MQ=255)

GATAAATATCAAACAATGTTTAATGTCATTATGGCGAATGCTTCTATTCTATTTTTTAGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATT  >  NC_000913/1229746‑1229929

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