BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	618,600	A→G	H375R (CAC→CGC)	fepE →	FepE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	618,600	0	A	G	81.4%	530.4 / 118.7	258	H375R (CAC→CGC)	fepE	FepE
Reads supporting (aligned to +/- strand): ref base A (22/26); new base G (104/106); total (126/132)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 259 total aligned reads displayed.

CGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGC  >  NC_000913/618511‑618692

cgcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCtt                                                                                           >  1:5534793/1‑93 (MQ=255)

 gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTa                                                                                          <  1:3892958/93‑1 (MQ=255)

 gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTa                                                                                          <  2:7425905/93‑1 (MQ=255)

 gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTa                                                                                          <  1:838902/93‑1 (MQ=255)

 gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTa                                                                                          <  1:305646/93‑1 (MQ=255)

 gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTa                                                                                          <  1:2201431/93‑1 (MQ=255)

    ttGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGtt                                                                                       >  1:476793/1‑93 (MQ=255)

    ttGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGtt                                                                                       <  2:7146156/93‑1 (MQ=255)

     tGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGttt                                                                                      >  2:6109104/1‑93 (MQ=255)

       aTCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTaa                                                                                    >  2:4603730/1‑93 (MQ=255)

       aTCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTaa                                                                                    >  1:1506588/1‑93 (MQ=255)

       aTCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTaa                                                                                    <  1:5315393/93‑1 (MQ=255)

         cggcggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACt                                                                                  <  1:6049783/93‑1 (MQ=255)

           gcggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTAt                                                                                >  2:2578492/1‑93 (MQ=255)

           gcggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTAt                                                                                >  1:6903065/1‑93 (MQ=255)

            cggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTAtt                                                                               <  1:7430612/93‑1 (MQ=255)

              ggATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTg                                                                             >  2:3426577/1‑93 (MQ=255)

               gATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGc                                                                            <  2:3321799/93‑1 (MQ=255)

               gATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGc                                                                            >  1:5815806/1‑93 (MQ=255)

                 tggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCgg                                                                          <  2:4780429/93‑1 (MQ=255)

                  ggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGt                                                                         <  1:7090428/93‑1 (MQ=255)

                  ggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGt                                                                         <  2:4705744/93‑1 (MQ=255)

                   gtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGtt                                                                        >  2:4750011/1‑93 (MQ=255)

                    tggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTg                                                                       >  2:5536421/1‑93 (MQ=255)

                     ggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGa                                                                      >  1:402156/1‑93 (MQ=255)

                     ggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGa                                                                      <  2:6757548/93‑1 (MQ=255)

                      gCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGaa                                                                     >  2:5640580/1‑93 (MQ=255)

                      gCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGaa                                                                     >  2:2181556/1‑93 (MQ=255)

                      gCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGaa                                                                     <  2:1535086/93‑1 (MQ=255)

                       cTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAg                                                                    <  2:6615148/93‑1 (MQ=255)

                       cTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAg                                                                    <  2:4219110/93‑1 (MQ=255)

                        ttGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAgg                                                                   <  2:470230/93‑1 (MQ=255)

                        ttGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAgg                                                                   >  2:3933667/1‑93 (MQ=255)

                        ttGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAgg                                                                   >  2:2349947/1‑93 (MQ=255)

                        ttGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAgg                                                                   >  2:6325127/1‑93 (MQ=255)

                         tgtgGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGt                                                                  >  1:1791754/1‑93 (MQ=255)

                           tggtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGc                                                                <  2:412804/93‑1 (MQ=255)

                            ggtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGcc                                                               <  1:463869/93‑1 (MQ=255)

                            ggtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGcc                                                               <  1:2715689/93‑1 (MQ=255)

                              tggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGt                                                             <  1:4552167/93‑1 (MQ=255)

                              tggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGt                                                             >  1:4169990/1‑93 (MQ=255)

                              tggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGt                                                             <  1:1084485/93‑1 (MQ=255)

                               ggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGtt                                                            <  2:1866707/93‑1 (MQ=255)

                                 cGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCg                                                          >  1:5460609/1‑93 (MQ=255)

                                  gTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGt                                                         >  1:4447478/1‑93 (MQ=255)

                                  gTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGt                                                         >  2:6863523/1‑93 (MQ=255)

                                  gTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGt                                                         >  1:1659984/1‑93 (MQ=255)

                                   tGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTc                                                        <  2:161930/93‑1 (MQ=255)

                                   tGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTc                                                        <  1:6586096/93‑1 (MQ=255)

                                    gCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTcc                                                       >  1:2409322/1‑93 (MQ=255)

                                     cTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCg                                                      >  2:2980822/1‑93 (MQ=255)

                                      tGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGa                                                     <  2:3600629/93‑1 (MQ=255)

                                       gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAg                                                    >  1:4529960/1‑93 (MQ=255)

                                       gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAg                                                    >  1:6533183/1‑93 (MQ=255)

                                       gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAg                                                    >  2:319140/1‑93 (MQ=255)

                                        ttGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGc                                                   <  1:5614896/93‑1 (MQ=255)

                                        ttGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGc                                                   >  2:3819169/1‑93 (MQ=255)

                                        ttGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGc                                                   <  2:3658349/93‑1 (MQ=255)

                                        ttGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGc                                                   <  2:3132878/93‑1 (MQ=255)

                                           cgcTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGc                                                <  1:25376/93‑1 (MQ=255)

                                             cTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcac                                              >  2:6552813/1‑93 (MQ=255)

                                              tATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcacc                                             >  1:3680821/1‑93 (MQ=255)

                                              tATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcacc                                             >  1:1024816/1‑93 (MQ=255)

                                               aTGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcacca                                            <  1:7502773/93‑1 (MQ=255)

                                               aTGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcacca                                            >  2:5391829/1‑93 (MQ=255)

                                                   gATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAg                                        <  2:2279858/93‑1 (MQ=255)

                                                    aTGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgc                                       <  2:7101807/93‑1 (MQ=255)

                                                     tGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcg                                      <  1:2866842/93‑1 (MQ=255)

                                                      ggCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcgg                                     <  2:1449453/93‑1 (MQ=255)

                                                       gCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcggc                                    <  2:4989435/93‑1 (MQ=255)

                                                          tCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGtt                                 >  1:6709015/1‑93 (MQ=255)

                                                          tCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGtt                                 >  1:4199841/1‑93 (MQ=255)

                                                          tCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGtt                                 >  1:7418610/1‑93 (MQ=255)

                                                           ccAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTc                                <  2:1241976/93‑1 (MQ=255)

                                                              gAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCgg                             >  2:649843/1‑93 (MQ=255)

                                                               aaaaCAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGc                            <  2:4397000/93‑1 (MQ=255)

                                                                 aaCAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCaa                          <  2:4572556/93‑1 (MQ=255)

                                                                  aCAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAAc                         <  1:3445753/93‑1 (MQ=255)

                                                                  aCAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAAc                         <  2:2758574/93‑1 (MQ=255)

                                                                   cAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACt                        <  2:2094492/93‑1 (MQ=255)

                                                                    aGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTg                       <  2:776683/93‑1 (MQ=255)

                                                                     ggATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTgg                      >  2:4620664/1‑93 (MQ=255)

                                                                      gATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGa                     <  2:4082073/93‑1 (MQ=255)

                                                                      gATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGa                     >  1:1430638/1‑93 (MQ=255)

                                                                       aTGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGAt                    <  1:1703429/93‑1 (MQ=255)

                                                                        tGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATc                   <  1:23104/93‑1 (MQ=255)

                                                                          ccATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGt                 >  2:1062322/1‑93 (MQ=255)

                                                                           cATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTc                >  2:1951172/1‑93 (MQ=255)

                                                                            atgatgGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCa               >  1:3711291/1‑93 (MQ=255)

                                                                             tgatgGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCaa              <  1:3283097/93‑1 (MQ=255)

                                                                             tgatgGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCaa              <  1:6164350/93‑1 (MQ=255)

                                                                                    aGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCAt       <  1:5133016/93‑1 (MQ=255)

                                                                                    aGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCAt       >  2:582919/1‑93 (MQ=255)

                                                                                     gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCTGGCAACTGGATCGTCAATGATCATg      >  2:5202153/1‑93 (MQ=255)

                                                                                     gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg      <  2:4356316/93‑1 (MQ=255)

                                                                                     gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg      <  1:3390052/93‑1 (MQ=255)

                                                                                     gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg      <  2:1314190/93‑1 (MQ=255)

                                                                                      aCCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGc     <  2:3985968/93‑1 (MQ=255)

                                                                                      aCCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGc     <  2:3250549/93‑1 (MQ=255)

                                                                                         aCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAgc  <  2:7375075/93‑1 (MQ=255)

CGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGC  >  NC_000913/618511‑618692

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