BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	619,925	T→C	T163A (ACC→GCC)	fepG ←	FepG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	619,925	0	T	C	99.2%	818.2 / NA	258	T163A (ACC→GCC)	fepG	FepG
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (120/135); total (121/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 261 total aligned reads displayed.

TTTTTGCCCACGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGA  >  NC_000913/619835‑620014

ttttttCCCACGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCAt                                                                                         <  1:92238/93‑1 (MQ=255)

tttttGCCCACGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCAt                                                                                         >  2:5565501/1‑93 (MQ=255)

tttttGCCCACGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCAt                                                                                         >  1:5556031/1‑93 (MQ=255)

 ttttGCCCACGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCAtt                                                                                        <  1:3481033/93‑1 (MQ=255)

   ttGCCCACGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTaa                                                                                      >  2:1143056/1‑93 (MQ=255)

   ttGCCCACGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTaa                                                                                      <  1:2303429/93‑1 (MQ=255)

    tGCCCACGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTaaa                                                                                     <  1:100211/93‑1 (MQ=255)

         ccGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGcca                                                                                <  2:6531492/92‑1 (MQ=255)

           gTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGccacc                                                                              >  2:2528886/1‑93 (MQ=255)

           gTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGccacc                                                                              <  1:5258789/93‑1 (MQ=255)

            tCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGccacca                                                                             <  2:6908661/93‑1 (MQ=255)

            tCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGccacca                                                                             <  2:1825575/93‑1 (MQ=255)

             cAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAg                                                                            <  1:1189207/93‑1 (MQ=255)

              aGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGc                                                                           >  2:4080313/1‑93 (MQ=255)

               ggCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCa                                                                          <  1:6527491/93‑1 (MQ=255)

               ggCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCa                                                                          <  2:4754236/93‑1 (MQ=255)

               ggCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCa                                                                          <  2:2249030/93‑1 (MQ=255)

               ggCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCa                                                                          >  1:7883785/1‑93 (MQ=255)

               ggCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCa                                                                          <  2:6398006/93‑1 (MQ=255)

                  cGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATgg                                                                       <  1:57805/93‑1 (MQ=255)

                   gTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGgc                                                                      <  2:5120743/93‑1 (MQ=255)

                    ttGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGgcg                                                                     >  1:3217306/1‑93 (MQ=255)

                     tGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGgcgc                                                                    >  2:3447157/1‑93 (MQ=255)

                       aGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGa                                                                  >  2:5948652/1‑93 (MQ=255)

                        gCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGaa                                                                 >  1:102387/1‑93 (MQ=255)

                        gCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGaa                                                                 <  2:3170620/93‑1 (MQ=255)

                         cGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAAc                                                                <  1:4674268/93‑1 (MQ=255)

                           aTCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGc                                                              >  1:850160/1‑93 (MQ=255)

                             ccGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccg                                                            <  1:2622137/93‑1 (MQ=255)

                              cGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGccga                                                           <  2:6587673/93‑1 (MQ=255)

                              cGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccga                                                           <  1:2747473/93‑1 (MQ=255)

                              cGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccga                                                           >  1:3633337/1‑93 (MQ=255)

                               ggCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGTCACCAGCATGGCGCGAACGccgaa                                                          >  1:7041434/1‑92 (MQ=255)

                               ggCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccgat                                                          >  2:2170205/1‑93 (MQ=255)

                               ggCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccgat                                                          <  2:160798/93‑1 (MQ=255)

                                 cATTCCACAAACTTGCTGTTATCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccgatac                                                        >  2:7025034/1‑93 (MQ=255)

                                  aTTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccgatacc                                                       <  1:7444049/93‑1 (MQ=255)

                                   ttCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccgataccg                                                      >  2:5716744/1‑93 (MQ=255)

                                   ttCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccgataccg                                                      >  1:4020115/1‑93 (MQ=255)

                                    tCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccgataccga                                                     <  2:2063673/93‑1 (MQ=255)

                                    tCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccgataccga                                                     <  2:1197792/93‑1 (MQ=255)

                                      cacaAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGccgataccgata                                                   <  1:4565239/93‑1 (MQ=255)

                                       acaAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGataa                                                  >  1:4336159/1‑93 (MQ=255)

                                       acaAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGataa                                                  >  1:356666/1‑93 (MQ=255)

                                       acaAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGACGCGAACGCCGATACCGataa                                                  >  2:6498531/1‑93 (MQ=255)

                                         aaaCCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGataata                                                <  2:2912975/93‑1 (MQ=255)

                                          aaCCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGAtaataa                                               <  2:824285/93‑1 (MQ=255)

                                            ccTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATc                                             >  1:827664/1‑93 (MQ=255)

                                            ccTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATc                                             >  1:7669405/1‑93 (MQ=255)

                                            ccTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATc                                             <  1:656204/93‑1 (MQ=255)

                                            ccTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATc                                             >  1:5725170/1‑93 (MQ=255)

                                             ctgctgTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCa                                            <  1:613163/93‑1 (MQ=255)

                                              tgctgTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCaa                                           >  1:4679711/1‑93 (MQ=255)

                                               gctgTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAAc                                          >  1:3114683/1‑93 (MQ=255)

                                               gctgTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAAc                                          <  1:5969276/93‑1 (MQ=255)

                                               gctgTTAACGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAAc                                          >  2:882074/1‑93 (MQ=255)

                                                  gTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGa                                       >  1:6054561/1‑93 (MQ=255)

                                                  gTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGa                                       <  1:7830484/93‑1 (MQ=255)

                                                  gTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGa                                       >  1:1827304/1‑93 (MQ=255)

                                                  gTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGa                                       >  2:3791075/1‑93 (MQ=255)

                                                    tAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGaaa                                     <  2:4409207/93‑1 (MQ=255)

                                                     aGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAg                                    <  1:1945956/93‑1 (MQ=255)

                                                          cGTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTc                               <  1:3855744/93‑1 (MQ=255)

                                                           gTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCg                              <  1:5679924/93‑1 (MQ=255)

                                                           gTTTCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCg                              >  1:2202632/1‑93 (MQ=255)

                                                            tttCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGa                             >  2:2189512/1‑93 (MQ=255)

                                                            tttCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGa                             <  1:5444255/93‑1 (MQ=255)

                                                            tttCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGa                             <  2:2831546/93‑1 (MQ=255)

                                                            tttCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGa                             <  1:2913190/93‑1 (MQ=255)

                                                            tttCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGa                             <  2:3201273/93‑1 (MQ=255)

                                                            tttCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGa                             <  1:6338215/93‑1 (MQ=255)

                                                             ttCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGAt                            <  2:4325588/93‑1 (MQ=255)

                                                             ttCTAAAGACGCTGTCAGCAACAGCCAGGCATTAGAGGCCACCAGCATGGCGCGAACGCCGATACCGATAAGAATCAACCGAAAGGTGTCGAt                            <  2:2727252/93‑1 (MQ=255)

                                                              tCTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATg                           >  2:4790384/1‑93 (MQ=255)

                                                               cTAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGc                          <  1:5290215/93‑1 (MQ=255)

                                                                tAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGcc                         >  2:6947135/1‑93 (MQ=255)

                                                                tAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGcc                         >  1:1309721/1‑93 (MQ=255)

                                                                tAAAGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGcc                         >  2:7885001/1‑93 (MQ=255)

                                                                 aaaGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCg                        <  2:5536929/93‑1 (MQ=255)

                                                                 aaaGACGCTTTCAGCAACAGCCAGGCATTAAAGCCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCg                        <  2:1746347/93‑1 (MQ=255)

                                                                  aaGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGt                       >  1:4830002/1‑93 (MQ=255)

                                                                  aaGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGt                       <  1:372504/93‑1 (MQ=255)

                                                                  aaGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGt                       <  1:6949671/93‑1 (MQ=255)

                                                                  aaGACGCTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGt                       >  2:7428481/1‑93 (MQ=255)

                                                                        cTTTCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCTATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGcc                 <  1:1617780/93‑1 (MQ=255)

                                                                           tCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAgg              >  1:4652908/1‑93 (MQ=255)

                                                                           tCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAgg              <  2:3621454/93‑1 (MQ=255)

                                                                           tCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAgg              >  2:7606754/1‑93 (MQ=255)

                                                                           tCAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAgg              <  2:1552044/93‑1 (MQ=255)

                                                                            cAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGc             <  2:4880763/93‑1 (MQ=255)

                                                                            cAGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGc             <  2:2959352/93‑1 (MQ=255)

                                                                             aGCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCg            >  1:2514043/1‑93 (MQ=255)

                                                                              gCAACAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGa           <  2:5392092/93‑1 (MQ=255)

                                                                                aaCAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGagc         <  1:295051/93‑1 (MQ=255)

                                                                                 aCAGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGagca        >  2:3883974/1‑93 (MQ=255)

                                                                                   aGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGagcagc      >  1:5354420/1‑93 (MQ=255)

                                                                                   aGCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGagcagc      <  1:5321334/93‑1 (MQ=255)

                                                                                    gCCAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGcc     >  2:3535854/1‑93 (MQ=255)

                                                                                     ccAGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCa    <  2:1375552/93‑1 (MQ=255)

                                                                                       aGGCATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGa  >  1:6741490/1‑93 (MQ=255)

TTTTTGCCCACGTCAGGCCGTTGAGCGATCCGGCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGA  >  NC_000913/619835‑620014

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