BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,392,958	G→A	G439S (GGC→AGC)	yjeF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,958	0	G	A	99.6%	843.4 / ‑5.7	231	G439S (GGC→AGC)	yjeF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (113/117); total (113/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 239 total aligned reads displayed.

CGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCG  >  NC_000913/4392869‑4393047

cGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGgcagca                                                                                        <  2:828832/93‑1 (MQ=255)

 ggTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCAt                                                                                       <  2:4398052/93‑1 (MQ=255)

  gTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATg                                                                                      <  1:722722/93‑1 (MQ=255)

   tGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATgg                                                                                     <  1:3991208/93‑1 (MQ=255)

      tGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCg                                                                                  >  2:4072337/1‑93 (MQ=255)

       gAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGa                                                                                 >  1:4524811/1‑93 (MQ=255)

       gAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGa                                                                                 <  2:717627/93‑1 (MQ=255)

       gAAAGGAGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGa                                                                                 >  2:3615513/1‑93 (MQ=255)

        aaaGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGAt                                                                                <  2:2053039/93‑1 (MQ=255)

        aaaGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGAt                                                                                >  1:2035395/1‑93 (MQ=255)

          aGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGAtgt                                                                              <  1:49852/93‑1 (MQ=255)

          aGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGAtgt                                                                              <  2:1180626/93‑1 (MQ=255)

          aGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGAtgt                                                                              >  2:4000127/1‑93 (MQ=255)

            gTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGc                                                                            <  1:2716523/93‑1 (MQ=255)

            gTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGc                                                                            <  1:1144002/93‑1 (MQ=255)

            gTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGc                                                                            <  2:1078300/93‑1 (MQ=255)

              gCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctc                                                                          <  2:395650/93‑1 (MQ=255)

                cGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctctc                                                                        <  1:606960/93‑1 (MQ=255)

                cGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctctc                                                                        >  2:1800750/1‑93 (MQ=255)

                  gagccGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTg                                                                      <  1:3333643/90‑1 (MQ=255)

                  gAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTg                                                                      <  2:4600246/93‑1 (MQ=255)

                   aaCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTgg                                                                     <  1:302978/93‑1 (MQ=255)

                   aaCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTgg                                                                     <  1:3060107/93‑1 (MQ=255)

                     ccGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGta                                                                   >  1:2939851/1‑93 (MQ=255)

                      cGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGtat                                                                  >  2:2565886/1‑93 (MQ=255)

                       gTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGtatt                                                                 >  2:46678/1‑93 (MQ=255)

                         ggTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGtattat                                                               <  2:3032187/93‑1 (MQ=255)

                           tCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTg                                                             <  2:1099187/93‑1 (MQ=255)

                              ccgccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGgcg                                                          <  1:1923147/93‑1 (MQ=255)

                                gccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCa                                                        >  1:3467230/1‑93 (MQ=255)

                                gccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCa                                                        >  2:480348/1‑93 (MQ=255)

                                gccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCa                                                        >  1:3733983/1‑93 (MQ=255)

                                 cccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAt                                                       >  1:1557449/1‑93 (MQ=255)

                                 cccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAt                                                       >  1:701441/1‑93 (MQ=255)

                                 cccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGGATTATTGGCGCAt                                                       >  2:3354319/1‑93 (MQ=255)

                                  ctaTCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAtt                                                      <  1:941770/91‑1 (MQ=255)

                                  ccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAtt                                                      <  1:3796258/93‑1 (MQ=255)

                                  ccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAtt                                                      <  1:178226/93‑1 (MQ=255)

                                  ccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAtt                                                      >  2:98139/1‑93 (MQ=255)

                                  ccATCCTGACGCTGTAGGCATTATTGATGCTGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTTCTCTCTGGTATTATTGGCGCAtt                                                      <  2:364151/93‑1 (MQ=255)

                                    aTCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGc                                                    >  2:3394086/1‑93 (MQ=255)

                                    aTCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGc                                                    <  1:1012933/93‑1 (MQ=255)

                                    aTCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGc                                                    <  2:866859/93‑1 (MQ=255)

                                    aTCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGc                                                    >  1:2318095/1‑93 (MQ=255)

                                     tCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCt                                                   >  2:1440680/1‑93 (MQ=255)

                                      ccTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCtt                                                  >  2:2889847/1‑93 (MQ=255)

                                        tGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTgg                                                <  1:4619976/93‑1 (MQ=255)

                                        tGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTgg                                                <  1:3188184/93‑1 (MQ=255)

                                        tGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTgg                                                <  1:3415947/93‑1 (MQ=255)

                                            gCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaa                                            <  1:1339755/93‑1 (MQ=255)

                                            gCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaa                                            >  2:344979/1‑93 (MQ=255)

                                              tttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaaa                                          >  2:510488/1‑93 (MQ=255)

                                               ttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATGATTGGCGCATTGCTTGGGCaaaaa                                         <  2:4518528/93‑1 (MQ=255)

                                                tAGGCATTATTGATGCCTGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAAc                                        <  1:1436488/93‑1 (MQ=255)

                                                tAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAAc                                        <  1:4101991/93‑1 (MQ=255)

                                                 aGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACt                                       <  1:3759706/93‑1 (MQ=255)

                                                  ggCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTg                                      >  1:3269100/1‑93 (MQ=255)

                                                    cATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTc                                    >  1:2334337/1‑93 (MQ=255)

                                                    cATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTc                                    <  2:2458303/93‑1 (MQ=255)

                                                    cATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTc                                    >  2:138582/1‑93 (MQ=255)

                                                      ttattGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGc                                  >  2:4177054/1‑93 (MQ=255)

                                                      ttattGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGc                                  <  2:3759528/93‑1 (MQ=255)

                                                        attGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCg                                >  2:1977182/1‑93 (MQ=255)

                                                         ttGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGt                               >  2:2136559/1‑93 (MQ=255)

                                                         ttGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGt                               <  2:3275590/93‑1 (MQ=255)

                                                           gATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTat                             <  2:2297770/93‑1 (MQ=255)

                                                              gCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAGAACTGTCGCCGTatgat                          <  2:1625592/93‑1 (MQ=255)

                                                              gCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTatgat                          <  1:4558481/93‑1 (MQ=255)

                                                              gCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTatgat                          <  2:1809536/93‑1 (MQ=255)

                                                                cGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGc                        <  1:1027664/93‑1 (MQ=255)

                                                                cGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGc                        >  2:2545078/1‑93 (MQ=255)

                                                                cGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGc                        <  1:3515060/93‑1 (MQ=255)

                                                                 ggAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCa                       >  1:2858607/1‑93 (MQ=255)

                                                                 ggAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCa                       <  1:4183416/93‑1 (MQ=255)

                                                                  gAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAg                      <  2:2277445/93‑1 (MQ=255)

                                                                  gAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAg                      >  1:3011885/1‑93 (MQ=255)

                                                                   aaaTGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGc                     <  2:4293314/93‑1 (MQ=255)

                                                                    aaTGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGcc                    >  2:1255089/1‑93 (MQ=255)

                                                                      tGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCtg                  >  2:1569261/1‑93 (MQ=255)

                                                                       gCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCtgt                 <  2:3900246/93‑1 (MQ=255)

                                                                        cAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCtgtg                <  1:2557756/93‑1 (MQ=255)

                                                                         aGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGc               >  2:2597381/1‑93 (MQ=255)

                                                                          ggCATGGCTAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCa              >  2:4548521/1‑93 (MQ=255)

                                                                           gCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAg             >  1:793529/1‑93 (MQ=255)

                                                                           gCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAg             >  1:2407693/1‑93 (MQ=255)

                                                                           gCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAg             <  2:3813216/93‑1 (MQ=255)

                                                                             aTGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGc           <  1:983350/93‑1 (MQ=255)

                                                                             aTGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGc           >  1:1306536/1‑93 (MQ=255)

                                                                              tGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCACAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCt          >  1:1827631/1‑93 (MQ=255)

                                                                              tGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCt          >  2:938942/1‑93 (MQ=255)

                                                                               ggCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtg         >  2:595612/1‑93 (MQ=255)

                                                                               ggCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtg         <  1:1063254/93‑1 (MQ=255)

                                                                               ggCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtg         <  2:58377/93‑1 (MQ=255)

                                                                                gCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt        >  2:3142677/1‑93 (MQ=255)

                                                                                   aGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTc     >  2:4134291/1‑93 (MQ=255)

                                                                                   aGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTc     <  2:2942514/93‑1 (MQ=255)

                                                                                    gcggcaGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcg    >  1:2681760/1‑93 (MQ=255)

                                                                                    gcggcaGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcg    >  2:2971332/1‑93 (MQ=255)

                                                                                     cggcagCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcgc   <  1:2113074/93‑1 (MQ=255)

                                                                                      ggcagcATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcgcg  <  1:3639848/93‑1 (MQ=255)

CGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCG  >  NC_000913/4392869‑4393047

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