BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,025,451	T→C	T202A (ACG→GCG)	b1956 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,025,451	0	T	C	83.8%	606.0 / 113.2	271	T202A (ACG→GCG)	b1956	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (21/23); new base C (106/121); total (127/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

Only 100 of 274 total aligned reads displayed.

GCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGA  >  NC_000913/2025361‑2025541

gCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTg                                                                                          <  2:5782903/93‑1 (MQ=255)

gCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTg                                                                                          <  2:1067787/93‑1 (MQ=255)

   aaGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGaa                                                                                       >  1:1374580/1‑93 (MQ=255)

   aaGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGaa                                                                                       <  1:4901449/93‑1 (MQ=255)

   aaGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGaa                                                                                       >  1:2797804/1‑93 (MQ=255)

       gAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTCAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACAtt                                                                                   <  1:6074749/93‑1 (MQ=255)

         aCCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTgc                                                                                 >  2:291140/1‑93 (MQ=255)

         aCCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTgc                                                                                 <  2:2059517/93‑1 (MQ=255)

         aCCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTgc                                                                                 <  2:631427/93‑1 (MQ=255)

          ccAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTgcg                                                                                >  2:1504257/1‑93 (MQ=255)

          ccAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTgcg                                                                                >  1:7363885/1‑93 (MQ=255)

           cAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTgcgc                                                                               >  2:5663432/1‑93 (MQ=255)

            aGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGcgcg                                                                              >  1:4839446/1‑93 (MQ=255)

             gCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGt                                                                             >  1:103843/1‑93 (MQ=255)

             gCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGt                                                                             >  2:7016471/1‑93 (MQ=255)

               gTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTaa                                                                           <  1:3322537/93‑1 (MQ=255)

                tACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTaaa                                                                          <  2:875599/93‑1 (MQ=255)

                tACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTaaa                                                                          >  1:4330884/1‑93 (MQ=255)

                 aCCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAAc                                                                         <  1:157959/93‑1 (MQ=255)

                    gcgcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAg                                                                      <  1:2490437/93‑1 (MQ=255)

                    gcgcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAg                                                                      <  1:5374525/93‑1 (MQ=255)

                     cgcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAgg                                                                     <  1:2385996/93‑1 (MQ=255)

                      gcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGt                                                                    <  2:4776893/93‑1 (MQ=255)

                      gcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGt                                                                    <  2:113214/93‑1 (MQ=255)

                      gcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGt                                                                    <  1:5972823/93‑1 (MQ=255)

                       cgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTc                                                                   <  2:6254142/93‑1 (MQ=255)

                       cgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTc                                                                   <  1:6357665/93‑1 (MQ=255)

                         gTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCgg                                                                 >  1:5824395/1‑93 (MQ=255)

                          tGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGc                                                                <  2:3660078/93‑1 (MQ=255)

                          tGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGc                                                                >  2:2936868/1‑93 (MQ=255)

                          ggACGATTTTCTCGTTGTGAATGGCCAATAAACCAAGGATGGGCAACAAAGCCATAATAAAACGCTGGAACATTGCGCGTAAACAAGGTAGGc                                                                >  1:5456240/2‑93 (MQ=255)

                            aCGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTg                                                              <  2:5513234/93‑1 (MQ=255)

                             cGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGc                                                             <  2:2151165/93‑1 (MQ=255)

                              gATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCg                                                            >  1:4125811/1‑93 (MQ=255)

                              gATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCg                                                            <  2:6235290/93‑1 (MQ=255)

                               aTTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGt                                                           >  1:5210539/1‑93 (MQ=255)

                               aTTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGt                                                           >  1:3987891/1‑93 (MQ=255)

                                ttttCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTc                                                          >  1:319653/1‑93 (MQ=255)

                                ttttCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTc                                                          >  2:6030887/1‑93 (MQ=255)

                                   tctcGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGaa                                                       >  1:5537791/1‑93 (MQ=255)

                                     tcGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAAc                                                     <  2:4357431/93‑1 (MQ=255)

                                     tcGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAAc                                                     <  2:7202329/93‑1 (MQ=255)

                                       gTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGt                                                   >  2:4921081/1‑93 (MQ=255)

                                         tgtgAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTaa                                                 <  1:2221542/93‑1 (MQ=255)

                                         tgtgAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTaa                                                 >  1:889300/1‑93 (MQ=255)

                                          gtgAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTaaa                                                <  2:1350756/93‑1 (MQ=255)

                                           tgAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTaaaa                                               <  2:7293153/93‑1 (MQ=255)

                                            gAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAAt                                              <  2:3946484/93‑1 (MQ=255)

                                            gAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAAt                                              <  2:7104026/93‑1 (MQ=255)

                                            gAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAAt                                              <  2:3546113/93‑1 (MQ=255)

                                            gAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAAt                                              <  2:3312322/93‑1 (MQ=255)

                                             aaTGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATc                                             >  1:1843342/1‑93 (MQ=255)

                                             aaTGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATc                                             <  1:5622078/93‑1 (MQ=255)

                                              aTGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATcc                                            <  1:3081312/93‑1 (MQ=255)

                                              aTGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATcc                                            >  1:424587/1‑93 (MQ=255)

                                              aTGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATcc                                            <  2:3154021/93‑1 (MQ=255)

                                               tGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCg                                           <  1:4651196/93‑1 (MQ=255)

                                               tGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCg                                           >  2:4411712/1‑93 (MQ=255)

                                               tGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCg                                           >  1:3835579/1‑93 (MQ=255)

                                                 gCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGc                                         <  1:658289/93‑1 (MQ=255)

                                                   cAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCAc                                       <  2:6826627/93‑1 (MQ=255)

                                                     aTAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCAcgc                                     >  1:1915991/1‑93 (MQ=255)

                                                       aaaCCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCga                                   <  1:838366/93‑1 (MQ=255)

                                                       aaaCCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCga                                   <  2:5343544/93‑1 (MQ=255)

                                                       aaaCCAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCga                                   <  1:6873227/93‑1 (MQ=255)

                                                          ccAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGac                                >  2:5401779/1‑93 (MQ=255)

                                                          ccAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACATAAAATCCGGCACGCGAGac                                >  2:5165268/1‑93 (MQ=255)

                                                           cAGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGaca                               >  1:7518401/1‑93 (MQ=255)

                                                            aGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacat                              >  2:4118971/1‑93 (MQ=255)

                                                            aGGGTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacat                              <  1:4707253/93‑1 (MQ=255)

                                                             gggTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacata                             <  2:2462847/93‑1 (MQ=255)

                                                             gggTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacata                             <  1:3597238/93‑1 (MQ=255)

                                                              ggTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacatac                            <  2:3425850/93‑1 (MQ=255)

                                                              ggTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacatac                            >  2:1822773/1‑93 (MQ=255)

                                                               gTTGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacataca                           <  2:3817696/93‑1 (MQ=255)

                                                                ttGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacatacat                          <  2:6519783/93‑1 (MQ=255)

                                                                ttGGGTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacatacat                          <  2:2622235/93‑1 (MQ=255)

                                                                   ggTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGc                       <  1:1100944/93‑1 (MQ=255)

                                                                   ggTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGc                       <  1:7264171/93‑1 (MQ=255)

                                                                    gTGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGcc                      <  2:5890621/93‑1 (MQ=255)

                                                                     tGACATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCt                     <  2:4389252/93‑1 (MQ=255)

                                                                      gACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTg                    >  1:7034969/1‑93 (MQ=255)

                                                                       aCATAGCCGTAATAACGCGCTGGAACATTGCGCGTTAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGt                   >  2:4524145/1‑93 (MQ=255)

                                                                       aCATAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGt                   <  2:2033110/93‑1 (MQ=255)

                                                                         aTAGCCGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTc                 >  2:2408841/1‑93 (MQ=255)

                                                                             ccGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAAcc             >  1:1499816/1‑93 (MQ=255)

                                                                              cGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCa            >  2:39815/1‑93 (MQ=255)

                                                                              cGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCa            <  1:7323799/93‑1 (MQ=255)

                                                                              cGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCa            <  1:2617999/93‑1 (MQ=255)

                                                                              cGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCa            >  1:3481948/1‑93 (MQ=255)

                                                                              cGTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCa            <  1:376543/93‑1 (MQ=255)

                                                                               gTAATAACGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCAt           >  2:91611/1‑93 (MQ=255)

                                                                                taataaCGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATc          >  1:6698104/1‑93 (MQ=255)

                                                                                taataaCGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATc          >  2:6891163/1‑93 (MQ=255)

                                                                                  ataaCGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcga        <  1:3623129/93‑1 (MQ=255)

                                                                                    aaCGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcgacg      >  1:2760843/1‑93 (MQ=255)

                                                                                     aCGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcgacga     <  2:4237786/93‑1 (MQ=255)

                                                                                     aCGCGCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcgacga     <  2:6152792/93‑1 (MQ=255)

                                                                                      cgcgCTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACgag    >  2:1577729/1‑93 (MQ=255)

                                                                                        cgTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAAGCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACgagga  <  2:4691721/93‑1 (MQ=255)

GCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGA  >  NC_000913/2025361‑2025541

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