BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,614,692	G→T	L133I (CTA→ATA)	yjjI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,614,692	0	G	T	99.8%	2156.4 / NA	543	L133I (CTA→ATA)	yjjI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base T (273/268); total (275/268)
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.99e-01

Only 100 of 572 total aligned reads displayed.

AGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGA  >  NC_000913/4614594‑4614789

aGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGGTaaaa                                                                                                 >  2:1357873/1‑101 (MQ=255)

  aCGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAtt                                                                                               <  1:7560703/101‑1 (MQ=255)

    gggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGGTATAATTCt                                                                                             >  2:151353/1‑101 (MQ=255)

     ggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGACCTCGGCTTGTGTTATAATTCTa                                                                                            >  1:12485677/1‑101 (MQ=255)

        nnntatTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACa                                                                                         <  2:14556819/98‑1 (MQ=255)

        cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTannn                                                                                         >  2:11957082/1‑98 (MQ=255)

        cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGGTCTGTCGAGGTAACGCCCGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTTACa                                                                                         >  1:8714742/1‑101 (MQ=255)

        cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGGTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGGGTTATAATTGTAACa                                                                                         >  2:6141296/1‑101 (MQ=255)

         gATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACAt                                                                                        <  1:11037357/101‑1 (MQ=255)

          atatTGGCGTGCATAAAGGCGTCTGGCAGGGATCTGCCGAGGAAACGCCAGAAACGATTTATTCGAAAATCGGTCTCGTCTTGTGGTGTAATTCTAACATa                                                                                       >  2:2913070/1‑101 (MQ=255)

           tatTGGCGTGCATAAAGGCGTCTGGCAGGGGTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                                                      >  1:6913477/1‑101 (MQ=255)

            atTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGGAACGCCAGAAACGTTTTATTCGAACATCGATCGCGGCTTGTGTTATAATTCTAAAATACg                                                                                     >  2:4911936/1‑101 (MQ=255)

             ttGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACgg                                                                                    >  1:3378025/1‑101 (MQ=255)

             ttGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGGAACGCCAGAAAAGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATannn                                                                                    >  2:13382248/1‑98 (MQ=255)

               ggCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGtt                                                                                  <  1:5726904/101‑1 (MQ=255)

               ggCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGtt                                                                                  >  1:10002609/1‑101 (MQ=255)

                gCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGGTg                                                                                 >  2:5229485/1‑101 (MQ=255)

                gCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGGAACGCCAGAACCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACggnnn                                                                                 >  2:7519316/1‑98 (MQ=255)

                  gTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGca                                                                               >  1:14571045/1‑101 (MQ=255)

                  gTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGca                                                                               <  2:5432680/101‑1 (MQ=255)

                   tGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaa                                                                              <  2:5699541/101‑1 (MQ=255)

                   tGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTAGTGTTATAATTCTAACATACGGTTGcaa                                                                              <  1:11765979/101‑1 (MQ=255)

                   tGCATAAAGGCGTCTGGCAGCGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCAATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcac                                                                              <  2:1943591/101‑2 (MQ=255)

                    gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACTGTTGcaac                                                                             >  1:7392857/1‑101 (MQ=255)

                     nnnaaaGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaaca                                                                            <  2:9711889/98‑1 (MQ=255)

                       tAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacnnn                                                                          >  2:9195492/1‑98 (MQ=255)

                       nnnaGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                                          <  2:10775484/98‑1 (MQ=255)

                        aaaGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                                                         <  1:12078162/101‑1 (MQ=255)

                        aaaGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                                                         <  1:7199251/101‑1 (MQ=255)

                          aGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACCTACGGTTGCAACAACGCa                                                                       <  2:1748814/101‑1 (MQ=255)

                          aGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCa                                                                       <  1:151166/101‑1 (MQ=255)

                           ggCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCAt                                                                      <  1:2500587/101‑1 (MQ=255)

                            gCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATc                                                                     >  1:542862/1‑101 (MQ=255)

                             cGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                                                    >  2:4667726/1‑101 (MQ=255)

                              nnntGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcaa                                                                   <  2:9473842/98‑3 (MQ=255)

                                cTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGt                                                                 <  1:5143769/101‑1 (MQ=255)

                                  nnnaGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                                               <  2:8065697/98‑1 (MQ=255)

                                  gtcAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                                               <  2:365706/99‑1 (MQ=255)

                                   gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGc                                                              <  1:4872489/101‑1 (MQ=255)

                                    cAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACAGCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCACTTGcc                                                             >  1:2336767/1‑101 (MQ=255)

                                     aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGccc                                                            >  1:13914696/1‑101 (MQ=255)

                                     aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGccc                                                            >  1:8385488/1‑101 (MQ=255)

                                     aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCACCCAGTTGccc                                                            >  2:3496084/1‑101 (MQ=255)

                                      gggTTCTGTCGAGGTAACGCCAGAACCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGcccc                                                           <  2:433515/101‑1 (MQ=255)

                                      gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGcccc                                                           <  2:6328269/101‑1 (MQ=255)

                                       ggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCa                                                          >  1:2224910/1‑101 (MQ=255)

                                        gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGcccnnn                                                         >  2:11563924/1‑98 (MQ=255)

                                        gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                                         >  1:12908402/1‑101 (MQ=255)

                                          tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                       <  1:5976931/101‑1 (MQ=255)

                                           cTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTa                                                      <  2:262800/101‑1 (MQ=255)

                                            tGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAg                                                     <  1:303755/101‑1 (MQ=255)

                                             gTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGtnnn                                                    >  2:14217605/1‑98 (MQ=255)

                                             gTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGa                                                    <  2:593695/101‑1 (MQ=255)

                                               nnnggTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGACCCAGGTAGAcc                                                  <  2:7661918/98‑1 (MQ=255)

                                               cGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAgnnn                                                  >  2:12840069/1‑98 (MQ=255)

                                                nnngTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCg                                                 <  2:12605687/98‑1 (MQ=255)

                                                nnngTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCACCCAGTTGCCCCAGGTAGACCg                                                 <  2:13227644/98‑1 (MQ=255)

                                                gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGannn                                                 >  2:10219739/1‑98 (MQ=255)

                                                gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGannn                                                 >  2:8566467/1‑98 (MQ=255)

                                                gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCg                                                 <  1:2752438/101‑1 (MQ=255)

                                                gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCg                                                 <  1:8070921/101‑1 (MQ=255)

                                                  ggTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                                               <  1:992779/101‑1 (MQ=255)

                                                    tAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCggnnn                                             >  2:8995281/1‑98 (MQ=255)

                                                     aaCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                                            <  1:7053561/101‑1 (MQ=255)

                                                       cGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTGTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGGAGAACGGCATCGa                                          >  1:13144628/1‑101 (MQ=255)

                                                          nnnaaaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                                       <  2:12289336/98‑1 (MQ=255)

                                                          cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                                       >  1:1012322/1‑101 (MQ=255)

                                                          cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGATGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                                       >  1:3802220/1‑101 (MQ=255)

                                                           nnnaaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAACGGCATCTAtgtg                                      <  2:11300119/98‑1 (MQ=255)

                                                           aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgtg                                      <  1:4359275/101‑1 (MQ=255)

                                                           aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgtg                                      <  1:9218518/101‑1 (MQ=255)

                                                            gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGa                                     >  2:3791288/1‑101 (MQ=255)

                                                              aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAcc                                   <  1:10800126/101‑1 (MQ=255)

                                                               aCGTTTTATTCGAACCTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCg                                  <  1:6586395/101‑1 (MQ=255)

                                                               aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGCCCGGCATCGATGTGACCg                                  <  2:6614961/101‑1 (MQ=255)

                                                               aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCg                                  >  1:12755505/1‑101 (MQ=255)

                                                                cGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAcnnn                                 >  2:10699098/1‑98 (MQ=255)

                                                                 gTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGAc                                >  2:2442977/1‑101 (MQ=255)

                                                                  nnntATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACg                               <  2:7939611/98‑1 (MQ=255)

                                                                  ggttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCgnnn                               >  2:12113230/3‑98 (MQ=255)

                                                                     tATTCGAACAGCCATCTCGTCTTGTGTTATAATTATAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTa                            <  2:7069286/101‑1 (MQ=255)

                                                                      aTTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAc                           <  1:1155860/101‑1 (MQ=255)

                                                                          nnncATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtgg                       <  2:7269983/98‑1 (MQ=255)

                                                                            aCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAACGGCATCGATGTGACCGACGGTACgtggcg                     >  1:8936045/1‑99 (MQ=255)

                                                                            aCATCCATCTCGTCTTGTGTTATAATTCTAACACACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCACCGATGTGTCCGCCGGTACgtggtg                     <  2:218411/101‑1 (MQ=255)

                                                                                  aTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAg               <  1:3646534/101‑1 (MQ=255)

                                                                                   tctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGAACGCGGTGGTaannn              >  2:14052791/1‑98 (MQ=255)

                                                                                       gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATgg          >  1:4632082/1‑101 (MQ=255)

                                                                                       gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATgg          <  1:1846309/101‑1 (MQ=255)

                                                                                       gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATgg          >  1:9595076/1‑101 (MQ=255)

                                                                                       gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGATGCCCCAGGTAGACCGGCCTCGATGTGACCGACGGTACGTGGTGGCAAAGAATgg          >  1:3261744/1‑101 (MQ=255)

                                                                                         cTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTc        >  2:1088930/1‑101 (MQ=255)

                                                                                         cTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTc        >  1:11144362/1‑101 (MQ=255)

                                                                                          ttGTGTTATAATTCTGACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTca       >  1:13404828/1‑101 (MQ=255)

                                                                                          ttGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTca       <  1:6393174/101‑1 (MQ=255)

                                                                                             tcttATAATTCTAACATACGGTTGCAACAACTCATCCAATTGCCCCAGGTAGACCGGCAACGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagca    <  2:3124012/99‑1 (MQ=255)

                                                                                              gtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagcag   >  2:3361608/1‑101 (MQ=255)

                                                                                              gtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGCGGTAAAGAATGGTcagcag   >  1:4206580/1‑101 (MQ=255)

                                                                                               ttATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCaga  <  1:8639139/101‑1 (MQ=255)

                                                                                               ttATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCaga  >  1:13687784/1‑101 (MQ=255)

AGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGA  >  NC_000913/4614594‑4614789

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