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	100.0%	1235.2 / NA	314	L133I (CTA→ATA)	yjjI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base T (162/152); total (162/152)

Only 100 of 332 total aligned reads displayed.

AGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGA  >  NC_000913/4614594‑4614789

aggcGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                                 <  1:5072772/98‑1 (MQ=255)

aGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                                 <  2:4067715/101‑1 (MQ=255)

aGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                                 <  2:2918915/101‑1 (MQ=255)

aGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                                 <  2:204114/101‑1 (MQ=255)

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

    gggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCt                                                                                             <  2:4627643/101‑1 (MQ=255)

    gggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCt                                                                                             <  1:408286/101‑1 (MQ=255)

     ggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                                            <  2:1148895/101‑1 (MQ=255)

     ggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGGTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                                            >  2:5153871/1‑101 (MQ=255)

      gCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTaa                                                                                           >  1:5597581/1‑101 (MQ=255)

      gCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTaa                                                                                           <  1:1608227/101‑1 (MQ=255)

        cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACa                                                                                         <  1:1240199/101‑1 (MQ=255)

        cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGGAACGCCAGAAAAGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACa                                                                                         >  2:6425784/1‑101 (MQ=255)

         gATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGGACATCGATCTCGTCTTGTGTTATAATTCTAACAt                                                                                        >  1:1774833/1‑101 (MQ=255)

           tatTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                                                      <  1:6575042/101‑1 (MQ=255)

            atTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACg                                                                                     <  1:1346251/101‑1 (MQ=255)

                 cGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGc                                                                                <  2:2393244/101‑1 (MQ=255)

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

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

                   tGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGGTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaa                                                                              >  1:5685302/1‑101 (MQ=255)

                    gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaac                                                                             >  1:5719665/1‑101 (MQ=255)

                    gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaac                                                                             <  2:5247374/101‑1 (MQ=255)

                    gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaac                                                                             <  1:3121159/101‑1 (MQ=255)

                     cATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaaca                                                                            >  2:5015345/1‑101 (MQ=255)

                     cATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaaca                                                                            <  1:2098942/101‑1 (MQ=255)

                       tAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                                          <  1:1131467/101‑1 (MQ=255)

                             cGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                                                    <  2:3928199/101‑1 (MQ=255)

                             cGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                                                    <  2:3008200/101‑1 (MQ=255)

                             cGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGGTGCAACAACGCATcc                                                                    >  1:1175/1‑101 (MQ=255)

                              gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCa                                                                   <  1:5913388/101‑1 (MQ=255)

                              gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCa                                                                   >  1:3532218/1‑101 (MQ=255)

                               tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAg                                                                  <  2:5330441/101‑1 (MQ=255)

                                cTGGCAGGGGTCTCTCGGGGGGACGTGAGAAACGTTTTATTCGAACCACGATCTCGTCCTTCGTTATTATTTTTAAAAACGGTTGCAATAACGCATCCAGt                                                                 >  2:2329953/1‑101 (MQ=255)

                                 tGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                                                <  2:221586/101‑1 (MQ=255)

                                    cAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGGTGCAACAACGCATCCAGTTGcc                                                             >  2:7001156/1‑101 (MQ=255)

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

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

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

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

                                          tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                       >  1:3946144/1‑101 (MQ=255)

                                          tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAGGATAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                       <  2:3332568/101‑1 (MQ=255)

                                           cTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTa                                                      >  1:4264486/1‑101 (MQ=255)

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

                                              tCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAc                                                   <  1:7059010/101‑1 (MQ=255)

                                              tCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAc                                                   <  2:2830316/101‑1 (MQ=255)

                                                gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCg                                                 >  2:6321130/1‑101 (MQ=255)

                                                gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCg                                                 >  1:3073932/1‑101 (MQ=255)

                                                 aGGTAACGCCAGAAACGTTTTATTCGAGCGGCAATCTGGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCgg                                                >  2:6897150/1‑101 (MQ=255)

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

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

                                                     aaCGCCAGAAACGTTTTATTCGAACAGCGGTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                                            >  1:6870535/1‑101 (MQ=255)

                                                       cGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGa                                          >  2:4463870/1‑101 (MQ=255)

                                                        gCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                                         >  2:5408235/1‑101 (MQ=255)

                                                        gCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                                         <  1:4666728/101‑1 (MQ=255)

                                                         ccAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtg                                        >  1:4858415/1‑101 (MQ=255)

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

                                                           aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgtg                                      >  2:130482/1‑101 (MQ=255)

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

                                                            gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCAACCAGTTGCCCCAGGTAGACCGGCATCGATGTGa                                     <  2:6451638/101‑1 (MQ=255)

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

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

                                                                  ttttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACg                               <  2:6615448/101‑1 (MQ=255)

                                                                  ttttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACg                               <  1:3571158/101‑1 (MQ=255)

                                                                  ttttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACg                               <  2:4212193/101‑1 (MQ=255)

                                                                     tATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTa                            <  2:1039878/101‑1 (MQ=255)

                                                                     tATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTa                            <  1:2295409/101‑1 (MQ=255)

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

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

                                                                      aTTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAc                           >  1:4601809/1‑101 (MQ=255)

                                                                      aTTCGAACATCGATCTCGTCTTGTGATATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAc                           <  1:661569/101‑1 (MQ=255)

                                                                       ttCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACg                          >  2:6768534/1‑101 (MQ=255)

                                                                       ttCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACg                          <  1:764909/101‑1 (MQ=255)

                                                                       ttCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACg                          >  1:1900911/1‑101 (MQ=255)

                                                                         cGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtg                        <  2:1588429/101‑1 (MQ=255)

                                                                         cGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtg                        >  1:2381827/1‑101 (MQ=255)

                                                                          gAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGGATCGATGTGACCGACGGTACgtgg                       >  1:6411757/1‑101 (MQ=255)

                                                                           aaCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtggt                      >  2:4300453/1‑101 (MQ=255)

                                                                           aaCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtggt                      <  2:3104881/101‑1 (MQ=255)

                                                                            aCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtggtg                     <  1:2470859/101‑1 (MQ=255)

                                                                             cATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGtggtgg                    <  2:190836/101‑1 (MQ=255)

                                                                                cGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaa                 <  2:2612451/101‑1 (MQ=255)

                                                                                 gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa                >  1:4404369/1‑101 (MQ=255)

                                                                                 gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa                <  1:3811531/101‑1 (MQ=255)

                                                                                 gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa                <  1:5255555/101‑1 (MQ=255)

                                                                                 gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa                >  1:4764154/1‑101 (MQ=255)

                                                                                   tctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGa              <  1:5414932/101‑1 (MQ=255)

                                                                                   tctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGa              <  1:3299226/101‑1 (MQ=255)

                                                                                    ctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGaa             <  2:5719665/101‑1 (MQ=255)

                                                                                     tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt            <  1:81736/101‑1 (MQ=255)

                                                                                     tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt            <  2:6817432/101‑1 (MQ=255)

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

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

                                                                                       gTCTTGTGTTATAATTCTAACATACGGATGCAACAACACATCCAGTTGCCCCAGGTAGTCCGGCATCGATGTGACCGGCGGTACGTGGTGGTAAAGAAtcg          <  2:5779864/101‑3 (MQ=255)

                                                                                       gTCTTGTGCTATAATTCTAACATACGGTTGCAACAACGCATCCAGTGGCACCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGaaggt          >  2:3277461/1‑98 (MQ=255)

                                                                                        tCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGt         <  1:7001156/101‑1 (MQ=255)

                                                                                         cTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGGGGGGGGAAAGAATGGTc        >  2:255961/1‑101 (MQ=255)

                                                                                            gtgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagc     >  1:6766747/1‑101 (MQ=255)

                                                                                             ngtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagca    <  2:3374059/100‑1 (MQ=255)

                                                                                              gtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagcag   <  2:366702/101‑1 (MQ=255)

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

AGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGA  >  NC_000913/4614594‑4614789

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