BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,486,147	(C)_7→8	coding (1485/1503 nt)	yjgR ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,486,140	1	.	C	96.0%	446.1 / 10.5	148	coding (1492/1503 nt)	yjgR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/4); new base C (95/47); total (97/51)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 178 total aligned reads displayed.

TAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACAC  >  NC_000913/4486049‑4486226

tAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTt‑ct‑‑c                                                                                        >  1:1380325/1‑93 (MQ=255)

  aaaTTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑ccc                                                                                      <  2:1668525/93‑1 (MQ=255)

   aaTTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑cccc                                                                                     >  1:1990303/1‑93 (MQ=255)

    aTTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑ccccc                                                                                    >  1:3284125/1‑93 (MQ=255)

      tCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑ccccc‑cc                                                                                 <  2:2722441/93‑1 (MQ=255)

        aaTATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑ccccc‑ccca                                                                               >  2:1137283/1‑91 (MQ=255)

         atatatTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑ccccc‑cccag                                                                              <  1:3897125/93‑4 (MQ=255)

          tatatTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑ccccc‑cccagc                                                                             >  1:4198919/1‑89 (MQ=255)

            tatTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCaa                                                                            >  2:2772123/1‑93 (MQ=255)

            tatTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCaa                                                                            >  1:3928438/1‑93 (MQ=255)

            tatTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCaa                                                                            >  1:4411677/1‑93 (MQ=255)

              ttGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAAc                                                                          >  2:2568287/1‑93 (MQ=255)

              ttGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAAc                                                                          >  2:523577/1‑93 (MQ=255)

                 cAGTGGGGTGGGCATGTTAAGACGCACCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTGCCTTCTT‑CTC‑CCCCC‑CCAGCAAACTcc                                                                       <  2:4650204/93‑1 (MQ=38)

                  aGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTccc                                                                      >  1:744037/1‑93 (MQ=255)

                  aGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTccc                                                                      <  1:147648/93‑1 (MQ=255)

                  aGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTccc                                                                      >  2:3084530/1‑93 (MQ=255)

                   gtgtggtgcgCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGCCAATTCTCGCCTGACTTACCTTCTTCCC‑‑CCCCC‑CCAGCAAACTcccc                                                                     <  2:758267/93‑1 (MQ=25)

                   gtggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTcccc                                                                     <  2:3660089/93‑1 (MQ=255)

                   gtggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTcccc                                                                     <  1:4576568/93‑1 (MQ=255)

                    tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCa                                                                    >  2:874895/1‑93 (MQ=255)

                     ggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCaa                                                                   <  2:3429142/93‑1 (MQ=255)

                     ggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCACCAAACTCCCCaa                                                                   >  1:2688777/1‑93 (MQ=255)

                       ggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACa                                                                 <  1:556326/93‑1 (MQ=255)

                       ggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACa                                                                 <  2:3048251/93‑1 (MQ=255)

                       ggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACa                                                                 <  1:2762597/93‑1 (MQ=255)

                        gtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACAt                                                                <  1:4102560/93‑1 (MQ=255)

                        gtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACAt                                                                >  1:2366328/1‑93 (MQ=255)

                          gggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATcc                                                              >  2:1250723/1‑93 (MQ=255)

                          gggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATcc                                                              >  1:241109/1‑93 (MQ=255)

                            gCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCcg                                                            >  2:3653141/1‑93 (MQ=255)

                            gCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCcg                                                            >  1:591402/1‑93 (MQ=255)

                            gCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCcg                                                            >  2:2625627/1‑93 (MQ=255)

                               tGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGt                                                         >  1:1210996/1‑93 (MQ=255)

                                gATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTa                                                        <  2:740483/93‑1 (MQ=255)

                                 aTAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTAc                                                       >  1:2651750/1‑93 (MQ=255)

                                 aTAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTAc                                                       >  1:1234196/1‑93 (MQ=255)

                                  tAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACa                                                      >  2:1532188/1‑93 (MQ=255)

                                  tAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACa                                                      >  1:1715624/1‑93 (MQ=255)

                                  tAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACa                                                      >  1:2278861/1‑93 (MQ=255)

                                  tAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACa                                                      >  1:4176257/1‑93 (MQ=255)

                                    aGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAAt                                                    <  2:8272/93‑1 (MQ=255)

                                       cgcgCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTCCCCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCt                                                  <  1:4540304/93‑1 (MQ=255)

                                        gcgcCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGa                                                >  1:4361922/1‑93 (MQ=255)

                                           ccAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTc                                             >  1:2193571/1‑93 (MQ=255)

                                           ccAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTc                                             >  1:932056/1‑93 (MQ=255)

                                              gCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTc                                          <  2:728144/93‑1 (MQ=255)

                                               cGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCa                                         <  2:2116057/93‑1 (MQ=255)

                                                gTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCAc                                        >  2:2038699/1‑93 (MQ=255)

                                                  ctcATCCGGCAAAGACAAATCTCCCCTGACTTACCTGCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACtt                                      <  1:2937056/91‑1 (MQ=255)

                                                   gCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTg                                     >  1:572131/1‑93 (MQ=255)

                                                   gCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTg                                     <  1:3995110/93‑1 (MQ=255)

                                                   gCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTgg                                    <  2:4082059/93‑1 (MQ=255)

                                                     accAGGAAAAGAAAAATCTCGCATGACTTACCTTGTC‑CTG‑CCCCC‑CCAGCAGCCTCCCCAACATCCCGCGTACAATCTGATTCGCCCCTTGGc                                   <  2:4362414/91‑1 (MQ=255)

                                                     aTCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGc                                   >  1:3795569/1‑93 (MQ=255)

                                                      tCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCg                                  >  2:394779/1‑93 (MQ=255)

                                                      tCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATACCGAGTACAATCTGATTCGTCACTTGGCg                                  >  2:1940117/1‑93 (MQ=255)

                                                       cAGCCAAAGACAAAGCTCGCCCGACTTACCTTCTTCCC‑‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGa                                 <  1:479012/93‑1 (MQ=255)

                                                        aGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAg                                >  1:3034/1‑93 (MQ=255)

                                                        aGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAg                                >  1:1865785/1‑93 (MQ=255)

                                                          gCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAgcg                              <  2:358479/93‑1 (MQ=255)

                                                             aaGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAgcggcg                           >  2:1371070/1‑93 (MQ=255)

                                                             aaGACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAgcggcg                           >  2:3149423/1‑93 (MQ=255)

                                                               gACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCt                         >  2:1831754/1‑93 (MQ=255)

                                                               gACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCt                         <  1:1396892/93‑1 (MQ=255)

                                                               gACAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCt                         >  1:2598500/1‑93 (MQ=255)

                                                                aCAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCtt                        <  1:4221178/93‑1 (MQ=255)

                                                                 cAAATCTCGCCTGCCTTACCTTCTT‑CT‑‑CCCCCTCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCttt                       <  1:2378470/93‑1 (MQ=255)

                                                                 cAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCttt                       <  1:2244162/93‑1 (MQ=255)

                                                                 cAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCttt                       <  1:1203575/93‑1 (MQ=255)

                                                                 cAAATCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCttt                       <  2:4267036/93‑1 (MQ=255)

                                                                  aaaTCTCGCCTTACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCtttt                      <  2:3009635/93‑1 (MQ=255)

                                                                  aaaTCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCtttt                      >  2:3396991/1‑93 (MQ=255)

                                                                  aaaTCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCtttt                      >  2:964469/1‑93 (MQ=255)

                                                                  aaaTCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCtttt                      <  1:2036817/93‑1 (MQ=255)

                                                                   aaTCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTg                     >  2:2334478/1‑93 (MQ=255)

                                                                   aaTCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTg                     >  1:2067834/1‑93 (MQ=255)

                                                                   aaTCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTg                     >  1:4400673/1‑93 (MQ=255)

                                                                    aTCTCGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTgg                    <  1:3100034/93‑1 (MQ=255)

                                                                       tcGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCa                 >  2:2175132/1‑93 (MQ=255)

                                                                       tcGCCTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCa                 >  2:3885708/1‑93 (MQ=255)

                                                                        cGCCTGACTTACCTTCTT‑CT‑‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCAtt               >  2:4133666/1‑93 (MQ=255)

                                                                        cGCCTGACTTACCTTCTC‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCAt                <  2:1579219/93‑1 (MQ=255)

                                                                          ccTGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTg              >  1:3873927/1‑93 (MQ=255)

                                                                            tGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGtt            >  1:3445303/1‑93 (MQ=255)

                                                                            tGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGtt            >  2:4355432/1‑93 (MQ=255)

                                                                            tGACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGAGGCGCTTTTGGCCATTGtt            >  2:1447419/1‑93 (MQ=255)

                                                                             gACTTCCCTTCTTCCC‑‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTAGGCCATTGttt           <  2:1245792/93‑1 (MQ=255)

                                                                             gACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGttt           >  1:3409362/1‑93 (MQ=255)

                                                                             gACTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGttt           >  1:954972/1‑93 (MQ=255)

                                                                              aCTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTg          >  1:3628563/1‑93 (MQ=255)

                                                                               cTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGc         >  2:3766675/1‑93 (MQ=255)

                                                                               cTTACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGc         >  2:4294001/1‑93 (MQ=255)

                                                                                ttACCTTCTT‑CTC‑CCCCC‑CCCGCAAAATCCCCAACAACCCCCGCACAAACTGATTCGTCACTTGGCGAGCGGGGCTTTTGGTCATTGTTTGca        >  2:2465569/1‑93 (MQ=255)

                                                                                ttACCTTCTT‑CTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGca        >  2:3284153/1‑93 (MQ=255)

                                                                                   ccTTCTT‑CTC‑CCCCC‑CCAGCCAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGTGCGGCGCTTTTGTTCATTCTTTGcacca     >  2:1015956/1‑93 (MQ=255)

                                                                                     ttctt‑cTC‑CCCCC‑CCAGCAAACTCCCCACCATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGcaccaca   <  1:1341406/93‑1 (MQ=255)

                                                                                     ttctt‑cTC‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGcaccaca   <  1:2825592/93‑1 (MQ=255)

                                                                                      tctt‑ctc‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTTcaccacac  >  2:1032524/1‑93 (MQ=255)

                                                                                      tctt‑ctc‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGcaccacac  >  2:3561724/1‑93 (MQ=255)

TAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT‑CT‑‑CCCCC‑CCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACAC  >  NC_000913/4486049‑4486226

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