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	97.9%	668.4 / ‑1.0	192	G439S (GGC→AGC)	yjeF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (3/1); new base A (93/95); total (96/96)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 210 total aligned reads displayed.

GCGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCG  >  NC_000913/4392868‑4393047

gCGGTGCTGAAAGGTGTCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCggcagc                                                                                         >  2:181736/1‑93 (MQ=255)

gCGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCggcagc                                                                                         >  2:793081/1‑93 (MQ=255)

gCGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCggcagc                                                                                         <  2:176473/93‑1 (MQ=255)

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

   gTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGAATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATg                                                                                      >  2:3920463/1‑93 (MQ=255)

    tGCTGAAAGGTCCCGGAACCGGGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATgg                                                                                     <  2:1340714/93‑1 (MQ=255)

     gCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATggg                                                                                    <  2:33259/93‑1 (MQ=255)

          aaGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGAtg                                                                               >  2:345351/1‑93 (MQ=255)

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

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

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

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

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

                ccGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctct                                                                         <  1:4639312/93‑1 (MQ=255)

                ccGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctct                                                                         >  2:4399532/1‑93 (MQ=255)

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

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

                  ggAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctctct                                                                       >  2:4181049/1‑93 (MQ=255)

                  ggAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGctctct                                                                       >  2:4557128/1‑93 (MQ=255)

                   gacccGTTGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTg                                                                      <  1:3328495/90‑1 (MQ=255)

                   gAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTg                                                                      >  2:226968/1‑93 (MQ=255)

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

                    aaCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTgg                                                                     >  1:2290968/1‑93 (MQ=255)

                     aCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGGTGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCTATGTGCTCTCTGGt                                                                    >  2:4059498/1‑93 (MQ=255)

                     aCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGt                                                                    >  2:3828935/1‑93 (MQ=255)

                      ccGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGta                                                                   <  2:4007861/93‑1 (MQ=255)

                        gTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGtatt                                                                 <  2:2006219/93‑1 (MQ=255)

                         tGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCCTGGGCGATGTGCTCTCTGGtatta                                                                >  1:670633/1‑93 (MQ=255)

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

                           gTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTattatt                                                              <  1:3822362/93‑1 (MQ=255)

                              tccgccCCTCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCCGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGgc                                                           <  2:4041487/92‑1 (MQ=255)

                              gccgccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGgc                                                           >  1:155383/1‑93 (MQ=255)

                              gccgccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGgc                                                           <  1:1506109/93‑1 (MQ=255)

                                  cccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAt                                                       <  1:223489/93‑1 (MQ=255)

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

                                  cccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGAATGGCGAGCGGCAGCATGGGCGATGTGATCTCTGGTATTATTTGCGCAt                                                       >  2:608832/1‑93 (MQ=255)

                                   ccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGTCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCAtt                                                      <  1:265720/93‑1 (MQ=255)

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

                                    cATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTg                                                     >  1:4709867/1‑93 (MQ=255)

                                    cATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTg                                                     >  1:1823661/1‑93 (MQ=255)

                                    cATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTg                                                     >  1:4677572/1‑93 (MQ=255)

                                     ctCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGc                                                    <  1:1300013/92‑1 (MQ=255)

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

                                       ccTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCtt                                                  <  2:2663645/93‑1 (MQ=255)

                                        cTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTg                                                 <  2:2968774/93‑1 (MQ=255)

                                        cTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTg                                                 >  2:4550324/1‑93 (MQ=255)

                                          gACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTggg                                               >  2:625988/1‑93 (MQ=255)

                                          gACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTggg                                               <  2:2129110/93‑1 (MQ=255)

                                          gACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTggg                                               <  2:751480/93‑1 (MQ=255)

                                           aCGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGc                                              <  1:3164682/93‑1 (MQ=255)

                                           aCGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGc                                              <  2:4661453/93‑1 (MQ=255)

                                            cGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCa                                             >  2:2346609/1‑93 (MQ=255)

                                              cTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaa                                           >  1:1396372/1‑93 (MQ=255)

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

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

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

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

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

                                                ttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaaaa                                         >  2:1389891/1‑93 (MQ=255)

                                                ttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaaaa                                         >  2:3058128/1‑93 (MQ=255)

                                                ttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCaaaaa                                         <  1:377318/93‑1 (MQ=255)

                                                  aGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACt                                       >  1:922955/1‑93 (MQ=255)

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

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

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

                                                      attattGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCg                                   >  1:4431731/1‑93 (MQ=255)

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

                                                        tattGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGcc                                 >  1:1464734/1‑93 (MQ=255)

                                                        tattGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGcc                                 >  1:395272/1‑93 (MQ=255)

                                                         tttGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCg                                <  1:4023676/92‑1 (MQ=255)

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

                                                           tGATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTa                              <  1:86151/93‑1 (MQ=255)

                                                            gATGCCGGAAATGCAGGCATGGCGAGCTGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTat                             <  1:3287816/93‑1 (MQ=255)

                                                            gATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTat                             >  1:3227304/1‑93 (MQ=255)

                                                            gATGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTat                             >  1:4338688/1‑93 (MQ=255)

                                                             aTGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTatg                            <  2:528229/93‑1 (MQ=255)

                                                             aTGCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTatg                            >  2:3609214/1‑93 (MQ=255)

                                                               gCCGGAAATGCAGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTatgat                          >  1:1917155/1‑93 (MQ=255)

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

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

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

                                                                          aGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGc               >  1:2987091/1‑93 (MQ=255)

                                                                          aGGCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGc               <  2:2387191/93‑1 (MQ=255)

                                                                           ggCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCa              <  2:570950/93‑1 (MQ=255)

                                                                           ggCATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCa              <  2:2214399/93‑1 (MQ=255)

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

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

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

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

                                                                             cATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAgg            <  1:4262567/93‑1 (MQ=255)

                                                                             cATGGCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAgg            <  2:3559178/93‑1 (MQ=255)

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

                                                                                 gCGAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt        <  1:4587659/93‑1 (MQ=255)

                                                                                   gAGCGGCAGTATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAGCTGTCGCCGTATGATGCAGCCTGTGCTGGCTgtgt      <  2:277379/93‑1 (MQ=255)

                                                                                   gAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTgtgt      <  2:1322369/93‑1 (MQ=255)

                                                                                   gAGCGGCAGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTgtgt      >  1:188246/1‑93 (MQ=255)

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

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

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

                                                                                    aGCGGCAGCATGGGCGATGTGCTATCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTc     <  2:2596872/93‑1 (MQ=255)

GCGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCG  >  NC_000913/4392868‑4393047

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