BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,158,827	G→T	G61V (GGC→GTC)	yijC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,158,827	0	G	T	99.7%	1298.8 / NA	372	G61V (GGC→GTC)	yijC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base T (157/214); minor base A (1/0); total (158/214)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 381 total aligned reads displayed.

GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCA  >  NC_000913/4158737‑4158917

gccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCa                                                                                          >  2:3343746/1‑93 (MQ=255)

  cgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAtt                                                                                        >  1:3962825/1‑93 (MQ=255)

  cgTTCGCTGGTGGAACCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAtt                                                                                        <  1:1220374/93‑1 (MQ=255)

   gtgcGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTg                                                                                       <  1:2852282/90‑1 (MQ=255)

   gTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTg                                                                                       <  2:3265875/93‑1 (MQ=255)

   gTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTg                                                                                       >  2:7399281/1‑93 (MQ=255)

    ttCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGc                                                                                      <  1:117433/93‑1 (MQ=255)

    ttCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGc                                                                                      >  1:4067009/1‑93 (MQ=255)

       gCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTcc                                                                                   >  1:3622416/1‑93 (MQ=255)

       ccTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTcc                                                                                   <  2:6387916/92‑1 (MQ=255)

        cTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTccc                                                                                  <  1:2281290/93‑1 (MQ=255)

         tggtggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCa                                                                                 <  1:5460106/93‑1 (MQ=255)

          ggtggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAc                                                                                >  1:2376701/1‑93 (MQ=255)

              gAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACctct                                                                            <  1:6979592/93‑1 (MQ=255)

               aaGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtt                                                                           >  2:358774/1‑93 (MQ=255)

               aaGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtt                                                                           >  2:6277794/1‑93 (MQ=255)

                 gCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtttt                                                                         >  2:4246849/1‑93 (MQ=255)

                 gCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtttt                                                                         <  1:3479624/93‑1 (MQ=255)

                   cGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTa                                                                       >  1:3245370/1‑93 (MQ=255)

                   cGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTa                                                                       <  1:2100355/93‑1 (MQ=255)

                   cGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTa                                                                       >  2:752664/1‑93 (MQ=255)

                    gCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTAt                                                                      <  2:6947950/93‑1 (MQ=255)

                         tAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGcg                                                                 <  1:6596871/93‑2 (MQ=255)

                         tAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCa                                                                 <  1:4168710/93‑1 (MQ=255)

                         tAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCa                                                                 <  1:721361/93‑1 (MQ=255)

                          aGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                                <  1:5664675/93‑1 (MQ=255)

                          aGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                                <  1:2958675/93‑1 (MQ=255)

                          aGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGGGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                                <  2:2488249/93‑1 (MQ=255)

                            ccAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTTAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAttt                                                              <  1:6405271/93‑1 (MQ=255)

                             cAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                             >  2:5785036/1‑93 (MQ=255)

                             cAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                             >  1:7325356/1‑93 (MQ=255)

                             cAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                             >  2:1183519/1‑93 (MQ=255)

                              aaTTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTcc                                                            <  2:5032555/93‑1 (MQ=255)

                              aaTTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTcc                                                            <  1:2531868/93‑1 (MQ=255)

                              aaTTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTcc                                                            <  2:4024320/93‑1 (MQ=255)

                               aTTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                                                           >  2:3459978/1‑93 (MQ=255)

                               aTTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                                                           >  1:3129999/1‑93 (MQ=255)

                               aTTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                                                           >  2:3110762/1‑93 (MQ=255)

                                  aaGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGa                                                        <  1:4870616/93‑1 (MQ=255)

                                  aaGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGa                                                        <  2:1026343/93‑1 (MQ=255)

                                   aGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc                                                       <  1:4478250/93‑1 (MQ=255)

                                    gTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACg                                                      <  1:4614983/93‑1 (MQ=255)

                                    gTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGCGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACg                                                      >  2:6461432/1‑93 (MQ=255)

                                     tGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt                                                     >  2:5525718/1‑93 (MQ=255)

                                      gCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCCCCTCTTTTTATCGGCATTTCCGCGACGTa                                                    <  1:1388153/93‑1 (MQ=255)

                                      gCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTa                                                    <  2:1937398/93‑1 (MQ=255)

                                        tGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGa                                                  <  1:6703891/93‑1 (MQ=255)

                                        tGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGa                                                  <  2:2090679/93‑1 (MQ=255)

                                            cGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGaa                                              <  2:1691346/93‑1 (MQ=255)

                                             gCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAAc                                             >  1:6059486/1‑93 (MQ=255)

                                             gCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAAc                                             <  1:1750212/93‑1 (MQ=255)

                                              cAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                            >  2:3367232/1‑93 (MQ=255)

                                              cAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                            <  2:4416990/93‑1 (MQ=255)

                                              cAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                            >  1:2149360/1‑93 (MQ=255)

                                               aGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg                                           >  1:5008856/1‑93 (MQ=255)

                                                gCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTgg                                          <  2:7754224/93‑1 (MQ=255)

                                                gCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTgg                                          >  1:919381/1‑93 (MQ=255)

                                                 cTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTggg                                         >  1:2669589/1‑93 (MQ=255)

                                                  ttCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGt                                        >  1:1622182/1‑93 (MQ=255)

                                                   tCGCCAGCCTGAGTTTGCGTGAAGTGGCGAGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTc                                       <  2:5368765/93‑1 (MQ=255)

                                                     gCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTg                                     <  1:4504368/93‑1 (MQ=255)

                                                        aGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGAcc                                  >  2:3987962/1‑93 (MQ=255)

                                                         gCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTTTCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCa                                 >  1:4299433/1‑93 (MQ=255)

                                                         gCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCa                                 >  2:1441207/1‑93 (MQ=255)

                                                         gCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCa                                 >  1:5770882/1‑93 (MQ=255)

                                                         gCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCa                                 >  1:5906553/1‑93 (MQ=255)

                                                         gCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCa                                 >  2:7008382/1‑93 (MQ=255)

                                                         gCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCa                                 <  2:2002618/93‑1 (MQ=255)

                                                           cTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATg                               <  2:2560778/93‑1 (MQ=255)

                                                           cTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATg                               >  1:4783028/1‑93 (MQ=255)

                                                            tGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATgg                              >  1:7580094/1‑93 (MQ=255)

                                                            tGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATgg                              >  1:5289115/1‑93 (MQ=255)

                                                             gAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGt                             <  2:2449425/93‑1 (MQ=255)

                                                             gAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGt                             >  1:3935338/1‑93 (MQ=255)

                                                              aGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGtt                            <  1:2583530/93‑1 (MQ=255)

                                                              aGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGtt                            <  2:263288/93‑1 (MQ=255)

                                                              aGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGtt                            <  2:6214910/93‑1 (MQ=255)

                                                               gTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtg                           >  2:6510303/1‑93 (MQ=255)

                                                               gTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtg                           >  2:1103550/1‑93 (MQ=255)

                                                                  tGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatg                        >  2:1831986/1‑93 (MQ=255)

                                                                   gCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatga                       <  1:1030355/93‑1 (MQ=255)

                                                                   gCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatga                       <  2:5183858/93‑1 (MQ=255)

                                                                    cGTGAAGTGGCTCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag                      <  1:3721547/93‑1 (MQ=255)

                                                                    cGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag                      <  2:4763769/93‑1 (MQ=255)

                                                                    cGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag                      <  1:5235731/93‑1 (MQ=255)

                                                                     gTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgaga                     <  2:4068179/93‑1 (MQ=255)

                                                                        aaGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCg                  >  2:6094609/1‑93 (MQ=255)

                                                                        aaGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCg                  >  2:3807647/1‑93 (MQ=255)

                                                                         aGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCgg                 <  2:727407/93‑1 (MQ=255)

                                                                           tGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGtt               <  1:7888375/93‑1 (MQ=255)

                                                                            ggCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGttt              >  2:665019/1‑93 (MQ=255)

                                                                             gcgcGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTa             <  1:4925557/93‑1 (MQ=255)

                                                                             gcgcGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTa             <  2:6859702/93‑1 (MQ=255)

                                                                              cgcgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTaa            <  1:5585481/93‑1 (MQ=255)

                                                                              cgcgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTaa            <  2:2849631/93‑1 (MQ=255)

                                                                                 gTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGc         <  1:6468733/93‑1 (MQ=255)

                                                                                 gTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGc         <  2:5941039/93‑1 (MQ=255)

                                                                                 gTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGc         <  2:4702738/93‑1 (MQ=255)

                                                                                      gCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGc    <  1:7774822/93‑1 (MQ=255)

                                                                                        gggCATTGCGCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCa  <  2:2359396/93‑1 (MQ=255)

GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCA  >  NC_000913/4158737‑4158917

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