BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,444,741	A→G	E1018G (GAG→GGG)	ytfN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,444,741	0	A	G	86.9%	705.4 / 96.0	289	E1018G (GAG→GGG)	ytfN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (18/20); new base G (119/132); total (137/152)
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.83e-01

Only 100 of 294 total aligned reads displayed.

TTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAGCGCGG  >  NC_000913/4444649‑4444832

  gggCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGc                                                                                           <  1:3217237/93‑1 (MQ=255)

   ggCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGcc                                                                                          <  2:6591749/93‑1 (MQ=255)

   ggCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGcc                                                                                          <  1:254035/93‑1 (MQ=255)

      gcgcATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGaa                                                                                       >  1:1453196/1‑93 (MQ=255)

       cgcATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGaaa                                                                                      <  2:3216006/93‑1 (MQ=255)

        gcATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGaaaa                                                                                     >  2:2369130/1‑93 (MQ=255)

        gcATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGaaaa                                                                                     <  1:5586050/93‑1 (MQ=255)

        gcATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGaaaa                                                                                     <  2:6560397/93‑1 (MQ=255)

        gcATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGaaaa                                                                                     <  1:950492/93‑1 (MQ=255)

         cATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAAc                                                                                    <  1:5137823/93‑1 (MQ=255)

           tCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACgg                                                                                  <  2:4009716/93‑1 (MQ=255)

           tCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACgg                                                                                  <  1:2670818/93‑1 (MQ=255)

           tCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACgg                                                                                  >  1:6430420/1‑93 (MQ=255)

              tggtgCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGt                                                                               <  2:3942939/93‑1 (MQ=255)

              tggtgCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGt                                                                               <  1:508746/93‑1 (MQ=255)

               ggtgCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTc                                                                              <  2:268019/93‑1 (MQ=255)

                gtgCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCg                                                                             <  2:7251456/93‑1 (MQ=255)

                gtgCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCg                                                                             >  2:5295660/1‑93 (MQ=255)

                 tgCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGa                                                                            <  1:2731869/93‑1 (MQ=255)

                    aCGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTc                                                                         >  1:901415/1‑93 (MQ=255)

                    aCGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTc                                                                         <  1:755160/93‑1 (MQ=255)

                     cGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTcc                                                                        <  2:3833878/93‑1 (MQ=255)

                     cGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTcc                                                                        >  1:2893914/1‑93 (MQ=255)

                      gATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCg                                                                       >  2:3575464/1‑93 (MQ=255)

                      gATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCg                                                                       >  1:4661026/1‑93 (MQ=255)

                       aTCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGa                                                                      <  1:3011704/93‑1 (MQ=255)

                       aTCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGa                                                                      >  1:3782535/1‑93 (MQ=255)

                       aTCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGa                                                                      >  1:585498/1‑93 (MQ=255)

                        tCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGAt                                                                     >  1:207918/1‑93 (MQ=255)

                        tCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGAt                                                                     <  1:3487390/93‑1 (MQ=255)

                         cTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGAtt                                                                    >  1:1743205/1‑93 (MQ=255)

                          tGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTa                                                                   <  2:5496226/93‑1 (MQ=255)

                           gCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTaa                                                                  >  1:4925916/1‑93 (MQ=255)

                              ggAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAg                                                               <  1:1084933/93‑1 (MQ=255)

                              ggAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAg                                                               <  2:16336/93‑1 (MQ=255)

                              ggAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAg                                                               <  2:6320348/93‑1 (MQ=255)

                              ggAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAg                                                               >  2:4073945/1‑93 (MQ=255)

                                  aGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAAc                                                           >  1:3925854/1‑93 (MQ=255)

                                   gcgcAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAAcc                                                          >  2:4661330/1‑93 (MQ=255)

                                   gcgcAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAAcc                                                          <  2:7047261/93‑1 (MQ=255)

                                   gcgcAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAAcc                                                          <  1:3394537/93‑1 (MQ=255)

                                    cgcAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCt                                                         <  1:5273904/93‑1 (MQ=255)

                                       aGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGAt                                                      <  1:2271131/93‑1 (MQ=255)

                                        gTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGAtt                                                     <  1:6278517/93‑1 (MQ=255)

                                         tAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTg                                                    >  1:1972142/1‑93 (MQ=255)

                                          aGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGt                                                   >  1:5208148/1‑93 (MQ=255)

                                           ggCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTc                                                  <  2:1965254/93‑1 (MQ=255)

                                              gTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCAc                                               <  2:951549/93‑1 (MQ=255)

                                              gTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCAc                                               >  2:228619/1‑93 (MQ=255)

                                              gTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCAc                                               >  2:6723849/1‑93 (MQ=255)

                                               tctcCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACg                                              >  2:7166367/1‑93 (MQ=255)

                                                ctcCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGt                                             <  1:3725928/93‑1 (MQ=255)

                                                 tcCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGtt                                            >  2:2212584/1‑93 (MQ=255)

                                                 tcCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGtt                                            <  2:621223/93‑1 (MQ=255)

                                                 tcCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGtt                                            <  2:4507359/93‑1 (MQ=255)

                                                   cAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTgg                                          >  2:2458618/1‑93 (MQ=255)

                                                   cAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTgg                                          <  1:5052623/93‑1 (MQ=255)

                                                   cAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTgg                                          >  1:6920741/1‑93 (MQ=255)

                                                    aGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGc                                         >  2:3366585/1‑93 (MQ=255)

                                                    aGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGc                                         >  1:4023554/1‑93 (MQ=255)

                                                     gCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGca                                        >  1:4384580/1‑93 (MQ=255)

                                                       gATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaac                                      >  1:4852641/1‑93 (MQ=255)

                                                         tgtgGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaacaa                                    >  1:2555872/1‑93 (MQ=255)

                                                         tgtgGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaacaa                                    <  2:6526593/93‑1 (MQ=255)

                                                          gtggtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAt                                   >  2:100363/1‑93 (MQ=255)

                                                          gtggtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAt                                   >  2:5113468/1‑93 (MQ=255)

                                                          gtggtgATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAt                                   <  2:5619093/93‑1 (MQ=255)

                                                           tggtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtg                                  >  2:4850906/1‑93 (MQ=255)

                                                           tggtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtg                                  >  2:131180/1‑93 (MQ=255)

                                                            ggtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtgt                                 <  2:605318/93‑1 (MQ=255)

                                                             gtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtgtg                                >  1:1389189/1‑93 (MQ=255)

                                                             gtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtgtg                                >  2:1223877/1‑93 (MQ=255)

                                                             gtgATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtgtg                                <  2:2725509/93‑1 (MQ=255)

                                                              ngATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTgc                               <  2:3435811/92‑1 (MQ=255)

                                                                 tGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCa                            >  2:363673/1‑93 (MQ=255)

                                                                 tGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCa                            <  2:808323/93‑1 (MQ=255)

                                                                 tGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCa                            <  2:1908922/93‑1 (MQ=255)

                                                                 tGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCa                            >  1:5999499/1‑93 (MQ=255)

                                                                  gCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCAt                           <  2:4272707/93‑1 (MQ=255)

                                                                  gCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCAt                           <  2:6689596/93‑1 (MQ=255)

                                                                   cTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCAtt                          >  1:5486038/1‑93 (MQ=255)

                                                                    ttAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTg                         >  2:5384390/1‑93 (MQ=255)

                                                                      aaCGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGAc                       <  1:6810628/93‑1 (MQ=255)

                                                                      aaCGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGAc                       <  2:6501777/93‑1 (MQ=255)

                                                                        cGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGc                     >  1:5644818/1‑93 (MQ=255)

                                                                         gATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGcc                    <  2:582557/93‑1 (MQ=255)

                                                                         gATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGcc                    <  2:6954708/93‑1 (MQ=255)

                                                                            aaCCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCttt                 <  2:3728619/93‑1 (MQ=255)

                                                                            aaCCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCttt                 <  2:2772562/93‑1 (MQ=255)

                                                                              ccTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTgg               >  1:6154007/1‑93 (MQ=255)

                                                                              ccTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTgg               >  1:7142559/1‑93 (MQ=255)

                                                                               cTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGc              <  1:6953497/93‑1 (MQ=255)

                                                                                tGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGcc             >  1:1615552/1‑93 (MQ=255)

                                                                                tGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGcc             <  2:2101660/93‑1 (MQ=255)

                                                                                 gCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCt            >  1:4124467/1‑93 (MQ=255)

                                                                                 gCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCt            >  1:125504/1‑93 (MQ=255)

                                                                                  cAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTg           >  1:1666818/1‑93 (MQ=255)

                                                                                   aaCCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGa          <  2:7487777/93‑1 (MQ=255)

                                                                                    aCCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGaa         >  2:5075899/1‑93 (MQ=255)

                                                                                    aCCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGaa         <  2:887680/93‑1 (MQ=255)

TTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAGCGCGG  >  NC_000913/4444649‑4444832

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