BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,694,167	C→T	G27G (GGC→GGT)	yhjS →	putative protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,694,167	0	C	T	84.3%	727.6 / 95.5	286	G27G (GGC→GGT)	yhjS	putative protease
Reads supporting (aligned to +/- strand): ref base C (29/16); new base T (126/115); total (155/131)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-03

Only 100 of 291 total aligned reads displayed.

AAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGTCGC  >  NC_000913/3694077‑3694256

aaaTGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGgtgt                                                                                         <  1:1221327/93‑1 (MQ=255)

  aTGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCt                                                                                       <  2:3235621/93‑1 (MQ=255)

   tGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCtg                                                                                      <  2:6266606/93‑1 (MQ=255)

   tGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCtg                                                                                      <  2:1126686/93‑1 (MQ=255)

     gTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCtggt                                                                                    <  2:4369922/93‑1 (MQ=255)

     gTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCtggt                                                                                    >  1:5091935/1‑93 (MQ=255)

       gACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCtggtgg                                                                                  <  2:37878/93‑1 (MQ=255)

       gACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCtggtgg                                                                                  <  2:3364240/93‑1 (MQ=255)

        aCATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTggtggt                                                                                 <  1:1897860/93‑1 (MQ=255)

          aTGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGttt                                                                               <  1:7316490/93‑1 (MQ=255)

           tGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTa                                                                              >  2:1303043/1‑93 (MQ=255)

            gAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTaa                                                                             <  2:6440314/93‑1 (MQ=255)

               ggACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGt                                                                          >  1:625812/1‑93 (MQ=255)

               ggACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGt                                                                          <  2:602420/93‑1 (MQ=255)

                gACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTc                                                                         >  2:4703283/1‑93 (MQ=255)

                 aCATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCg                                                                        <  1:4368732/93‑1 (MQ=255)

                   aTTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGAt                                                                      <  1:5972610/93‑1 (MQ=255)

                   aTTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGCTTAACGTCGAt                                                                      >  2:4263242/1‑93 (MQ=255)

                    ttGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATc                                                                     <  1:68067/93‑1 (MQ=255)

                    ttGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATc                                                                     <  1:642014/93‑1 (MQ=255)

                    ttGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATc                                                                     >  1:7072298/1‑93 (MQ=255)

                         gACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCAt                                                                <  1:3994906/93‑1 (MQ=255)

                         gACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCAt                                                                <  2:3957202/93‑1 (MQ=255)

                          aCCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATg                                                               >  2:1548453/1‑93 (MQ=255)

                           cccTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGa                                                              >  1:214733/1‑93 (MQ=255)

                            ccTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGaa                                                             <  2:2342903/93‑1 (MQ=255)

                            ccTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGaa                                                             >  2:276989/1‑93 (MQ=255)

                             cTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAg                                                            >  1:4093808/1‑93 (MQ=255)

                              tGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGa                                                           >  1:7487069/1‑93 (MQ=255)

                              tGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGa                                                           <  1:2905137/93‑1 (MQ=255)

                              tGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGa                                                           <  1:1294154/93‑1 (MQ=255)

                              tGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGAGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGa                                                           <  2:4525713/93‑1 (MQ=255)

                               gTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGAt                                                          >  1:5471949/1‑93 (MQ=255)

                               gTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGAt                                                          >  1:6126710/1‑93 (MQ=255)

                                tATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATg                                                         <  1:1647882/93‑1 (MQ=255)

                                 aTTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGc                                                        <  2:5656879/93‑1 (MQ=255)

                                  ttCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCt                                                       >  1:1988359/1‑93 (MQ=255)

                                   tctcTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTa                                                      >  1:5589269/1‑93 (MQ=255)

                                    ctctATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTAt                                                     <  1:172824/93‑1 (MQ=255)

                                     tctATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATc                                                    <  1:2568927/93‑1 (MQ=255)

                                      ctATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCa                                                   >  2:3314736/1‑93 (MQ=255)

                                       tATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAg                                                  >  2:139111/1‑93 (MQ=255)

                                        aTCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGt                                                 >  1:7169727/1‑93 (MQ=255)

                                        aTCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGt                                                 >  2:1796391/1‑93 (MQ=255)

                                        aTCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGt                                                 >  2:1194864/1‑93 (MQ=255)

                                         tCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTc                                                >  1:3909324/1‑93 (MQ=255)

                                         tCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTc                                                <  2:1990745/93‑1 (MQ=255)

                                          cGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCt                                               <  2:4918028/93‑1 (MQ=255)

                                           ggTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTg                                              >  1:6900305/1‑93 (MQ=255)

                                            gTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTgg                                             >  1:5669082/1‑93 (MQ=255)

                                             tATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGc                                            <  2:2537551/93‑1 (MQ=255)

                                             tATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGc                                            >  1:2250458/1‑93 (MQ=255)

                                              aTCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCg                                           >  1:1112339/1‑93 (MQ=255)

                                              aTCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCg                                           >  2:4096953/1‑93 (MQ=255)

                                               tctcATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGa                                          >  2:2151299/1‑93 (MQ=255)

                                               tctcATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGa                                          >  1:2812859/1‑93 (MQ=255)

                                               tctcATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGa                                          >  1:5055208/1‑93 (MQ=255)

                                                ctcATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGaa                                         <  2:5057094/93‑1 (MQ=255)

                                                 tcatcaTTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAAt                                        >  2:3485892/1‑93 (MQ=255)

                                                 tcatcaTTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAAt                                        >  2:7298492/1‑93 (MQ=255)

                                                  catcatTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATc                                       >  2:6922123/1‑93 (MQ=255)

                                                    tcatTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCaa                                     >  1:5186351/1‑93 (MQ=255)

                                                    tcatTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCaa                                     >  2:7388695/1‑93 (MQ=255)

                                                     catTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCaaa                                    >  2:5202025/1‑93 (MQ=255)

                                                     catTATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCaaa                                    <  1:5921904/93‑1 (MQ=255)

                                                       ttATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACa                                  >  2:6482706/1‑93 (MQ=255)

                                                        tATGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACaa                                 <  2:1029322/93‑1 (MQ=255)

                                                         aTGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAAt                                >  1:1258964/1‑93 (MQ=255)

                                                         aTGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAAt                                >  2:421930/1‑93 (MQ=255)

                                                          tGGGATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAAtt                               >  2:67410/1‑93 (MQ=255)

                                                           gggATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTg                              >  2:4037526/1‑93 (MQ=255)

                                                            ggATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGc                             <  1:3025121/93‑1 (MQ=255)

                                                            ggATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGc                             <  1:4715000/93‑1 (MQ=255)

                                                             gATGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCa                            <  2:7183469/93‑1 (MQ=255)

                                                             gATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCa                            <  2:793227/93‑1 (MQ=255)

                                                              aTGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCAt                           <  2:5122724/93‑1 (MQ=255)

                                                               tGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATc                          >  2:2687129/1‑93 (MQ=255)

                                                               tGAGCTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATc                          >  2:312144/1‑93 (MQ=255)

                                                               tGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATc                          >  1:3788837/1‑93 (MQ=255)

                                                               tGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATc                          <  2:7492265/93‑1 (MQ=255)

                                                                  gCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCa                       >  1:5798290/1‑93 (MQ=255)

                                                                  gCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCa                       <  2:6637690/93‑1 (MQ=255)

                                                                   cTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAg                      <  1:887996/93‑1 (MQ=255)

                                                                   cTGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAg                      <  1:6080140/93‑1 (MQ=255)

                                                                    tGCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAgg                     >  2:1840034/1‑93 (MQ=255)

                                                                     gCGACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGc                    <  1:5282531/93‑1 (MQ=255)

                                                                       gACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTg                  >  1:6602630/1‑93 (MQ=255)

                                                                       gACATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTg                  >  1:5680241/1‑93 (MQ=255)

                                                                        aCATATGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGa                 <  1:3682814/93‑1 (MQ=255)

                                                                           tatGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAAc              <  2:5468159/93‑1 (MQ=255)

                                                                           tatGCCAGCAGGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAAc              <  2:894120/93‑1 (MQ=255)

                                                                           tatGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAAc              >  1:2543679/1‑93 (MQ=255)

                                                                                 agcagGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGcaca        <  2:3606184/93‑1 (MQ=255)

                                                                                  gcagGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCacac       <  1:3368783/93‑1 (MQ=255)

                                                                                  gcagGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCacac       <  1:3160470/93‑1 (MQ=255)

                                                                                  gcagGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCacac       >  1:1475320/1‑93 (MQ=255)

                                                                                   cagGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACg      >  2:833727/1‑93 (MQ=255)

                                                                                    agGCGGTGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGt     >  2:1668470/1‑93 (MQ=255)

                                                                                       cggtgTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGTcgc  <  1:3318024/93‑1 (MQ=255)

                                                                                       cggcgTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGTcgc  <  1:4930361/93‑1 (MQ=255)

AAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGTCGC  >  NC_000913/3694077‑3694256

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: `‑`