BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,158,202	+C	coding (961/975 nt)	ygiRQ ←	Fe‑S oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,158,201	1	.	C	99.2%	844.0 / ‑1.5	235	coding (962/975 nt)	ygiRQ	Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (100/133); total (100/135)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 261 total aligned reads displayed.

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158109‑3158288

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                       <  1:783110/93‑2 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                       <  1:4648872/93‑2 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                       <  2:3501618/93‑2 (MQ=255)

   cTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cct                                                                                      >  2:747268/1‑91 (MQ=255)

   cTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cct                                                                                      <  1:4443497/93‑3 (MQ=255)

   cTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cct                                                                                      >  2:3366834/1‑91 (MQ=255)

     gTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cctgg                                                                                    >  1:967443/1‑89 (MQ=255)

     gTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cctgg                                                                                    >  2:3019002/1‑89 (MQ=255)

       ccATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggt                                                                                   <  2:6201587/93‑1 (MQ=255)

        cATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggtt                                                                                  <  2:4571178/93‑1 (MQ=255)

        cATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggtt                                                                                  >  2:2852473/1‑93 (MQ=255)

         aTCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttg                                                                                 <  2:1278891/93‑1 (MQ=255)

         aTCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttg                                                                                 >  2:4859511/1‑93 (MQ=255)

         aTCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttg                                                                                 <  1:3238827/93‑1 (MQ=255)

         aTCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttg                                                                                 <  1:7539/93‑1 (MQ=255)

          tctcTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttgg                                                                                >  1:1800038/1‑93 (MQ=255)

          tctcTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTCggttgg                                                                                >  1:2442019/1‑93 (MQ=255)

           ctctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttggt                                                                               <  1:5636646/93‑1 (MQ=255)

           ctctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGGGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttggt                                                                               >  2:6564581/1‑93 (MQ=255)

             ctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGttt                                                                             <  1:1505500/93‑1 (MQ=255)

              ttCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTc                                                                            <  1:3720581/93‑1 (MQ=255)

              ttCGGTAGAAAGCGGAATAGCAGGCGGGGTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGGTGGTTTc                                                                            >  1:2610016/1‑93 (MQ=255)

              ttCGGTAGAAAGCGGAATAGCAGGCGGGATGATCCACACATAGCGGTCGCCGAGCTTTTGTATCACTGCGCAGGCACA‑GACTGGGCGGGtttg                                                                            >  1:3325341/1‑92 (MQ=255)

               tCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCg                                                                           >  2:4978458/1‑93 (MQ=255)

               tCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCg                                                                           <  2:5607075/93‑1 (MQ=255)

                   tAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtggt                                                                       <  1:2538966/93‑1 (MQ=255)

                    aGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtggtg                                                                      <  2:6105424/93‑1 (MQ=255)

                     gAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGc                                                                     <  2:3286499/93‑1 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCa                                                                    <  1:6567560/93‑1 (MQ=255)

                       aaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAg                                                                   >  2:88121/1‑93 (MQ=255)

                         gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGaa                                                                 >  1:4270715/1‑93 (MQ=255)

                         gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGaa                                                                 <  2:6338277/93‑1 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATag                                                              >  2:594366/1‑92 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAc                                                              <  2:478154/93‑1 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAc                                                              <  1:824775/93‑1 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAc                                                              <  2:1832054/93‑1 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAc                                                              <  2:5076260/93‑1 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAc                                                              >  1:5364809/1‑93 (MQ=255)

                             aaTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAcg                                                             <  2:589824/93‑1 (MQ=255)

                              aTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAcgc                                                            >  2:5166549/1‑93 (MQ=255)

                              aTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAcgc                                                            <  1:4738033/93‑1 (MQ=255)

                               tAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                                           <  1:292768/93‑1 (MQ=255)

                                   aGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                                       >  2:3097190/1‑93 (MQ=255)

                                   aGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                                       <  1:3602297/93‑1 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     <  1:6634070/93‑1 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                   >  2:4432403/1‑93 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                   <  1:5032774/93‑1 (MQ=255)

                                       cgcttGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACACC‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                   <  1:899751/90‑1 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                  >  1:4355422/1‑93 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                  <  1:2991074/93‑1 (MQ=255)

                                          ttGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                                                >  2:2030183/1‑93 (MQ=255)

                                           tGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGt                                               <  1:3427147/93‑1 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTAc                                             >  1:2021328/1‑93 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTAc                                             <  2:694954/93‑1 (MQ=255)

                                              tCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCGCA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCTAAGCATGGGCGTACa                                            <  1:3522858/93‑1 (MQ=255)

                                                cacacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                          <  1:1437818/93‑1 (MQ=255)

                                                cacacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                          <  1:2748391/93‑1 (MQ=255)

                                                 acacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCa                                         >  1:2312472/1‑93 (MQ=255)

                                                  cacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCAc                                        >  1:6130771/1‑93 (MQ=255)

                                                   acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACt                                       >  2:4894682/1‑93 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      <  2:2794684/93‑1 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      <  1:3879843/93‑1 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      >  1:2840552/1‑93 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      <  2:4157339/93‑1 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      <  2:6309339/93‑1 (MQ=255)

                                                     aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                     >  2:2268571/1‑93 (MQ=255)

                                                       aGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTAt                                   <  2:360594/93‑1 (MQ=255)

                                                       aGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTAt                                   >  1:6355683/1‑93 (MQ=255)

                                                         cGGTCGCCTTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTGCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCg                                 <  1:2244974/93‑1 (MQ=255)

                                                         cGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCg                                 >  2:5348043/1‑93 (MQ=255)

                                                         cGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCg                                 >  1:5918566/1‑93 (MQ=255)

                                                          ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc                                >  1:6011505/1‑93 (MQ=255)

                                                             cGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCt                             <  1:1154526/93‑1 (MQ=255)

                                                             cGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCt                             <  1:5727165/93‑1 (MQ=255)

                                                              gCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCtt                            >  1:2115772/1‑93 (MQ=255)

                                                              gCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCtt                            <  1:4124006/93‑1 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           <  1:164114/93‑1 (MQ=255)

                                                                cGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCa                          <  1:4811436/93‑1 (MQ=255)

                                                                 gtgtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc                         <  2:3445173/93‑1 (MQ=255)

                                                                   gtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtt                       <  2:5583049/93‑1 (MQ=255)

                                                                   gtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtt                       <  2:1687700/93‑1 (MQ=255)

                                                                       ttGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtctc                   <  2:2369651/93‑1 (MQ=255)

                                                                       ttGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtctc                   <  1:1907865/93‑1 (MQ=255)

                                                                         gCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGa                 <  2:6487885/93‑1 (MQ=255)

                                                                           aTCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaaa               >  1:5483616/1‑93 (MQ=255)

                                                                            tCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaag              >  1:4522647/1‑93 (MQ=255)

                                                                             cAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaaga             >  2:1471620/1‑93 (MQ=255)

                                                                             cAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaaga             <  1:5803238/93‑1 (MQ=255)

                                                                              aaTGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGATGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            <  2:4308322/93‑1 (MQ=255)

                                                                              aaTGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            <  2:2681148/93‑1 (MQ=255)

                                                                              aaTGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            <  1:3928210/93‑1 (MQ=255)

                                                                                 gcgcGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgc         <  1:1032437/93‑1 (MQ=255)

                                                                                 gcgcGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgc         <  2:3758428/93‑1 (MQ=255)

                                                                                 gcgcGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgc         <  1:5152556/93‑1 (MQ=255)

                                                                                  cgcgGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgca        >  2:2498456/1‑93 (MQ=255)

                                                                                    cgGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcagc      <  2:6451839/93‑1 (MQ=255)

                                                                                    cgGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcagc      <  2:2628439/93‑1 (MQ=255)

                                                                                      ggCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCaa    >  2:1928379/1‑93 (MQ=255)

                                                                                       gCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAac   >  2:4511295/1‑93 (MQ=255)

                                                                                        caca‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAaca  <  1:137062/93‑1 (MQ=255)

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158109‑3158288

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