BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,175,853	G→A	V178I (GTT→ATT)	nusG →	NusG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,175,853	0	G	A	96.4%	573.0 / 3.4	167	V178I (GTT→ATT)	nusG	NusG
Reads supporting (aligned to +/- strand): ref base G (2/4); new base A (77/84); total (79/88)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 174 total aligned reads displayed.

GTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTT  >  NC_000913/4175763‑4175942

gttgAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGAtt                                                                                         <  2:1295030/93‑1 (MQ=255)

  tgAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGa                                                                                       <  1:3289670/93‑1 (MQ=255)

   gaagaaGTGGATTACGTGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCGGGACTTCAGCCAGAATGaa                                                                                      >  2:4217630/1‑93 (MQ=255)

   gaagaaGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCCGGTTGaa                                                                                      >  2:3106316/1‑93 (MQ=255)

   gaagaaGTGGATTACGAGAAATCTCGTATGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGaa                                                                                      >  2:705242/1‑93 (MQ=255)

     agaagTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGaaaa                                                                                    >  1:3344842/1‑93 (MQ=255)

     agaagTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGaaaa                                                                                    <  1:3716363/93‑1 (MQ=255)

      gaagTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGaaaaa                                                                                   >  1:4473832/1‑93 (MQ=255)

      gaagTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGaaaaa                                                                                   >  2:4462176/1‑93 (MQ=255)

      aaagTGGAGTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGCCCCCGGTAGAGCTGGACTTCAGCCAGATTGaaaaa                                                                                   <  2:3348474/92‑1 (MQ=255)

        agTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGc                                                                                 <  1:3389952/93‑1 (MQ=255)

           ggATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTa                                                                              <  2:2689509/93‑1 (MQ=255)

               tACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAAccc                                                                          <  1:4747831/93‑1 (MQ=255)

               tACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAAccc                                                                          <  2:1937469/93‑1 (MQ=255)

               tACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAAccc                                                                          <  1:1234601/93‑1 (MQ=255)

                 cGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAg                                                                        <  1:3894961/93‑1 (MQ=255)

                  gagaAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGc                                                                       <  1:4367517/93‑1 (MQ=255)

                  gagaAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGc                                                                       <  2:1802941/93‑1 (MQ=255)

                  gagaAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGc                                                                       >  1:2547111/1‑93 (MQ=255)

                    gaAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGa                                                                     >  1:1052689/1‑93 (MQ=255)

                      aaTCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATc                                                                   <  1:2251158/93‑1 (MQ=255)

                      aaTCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATc                                                                   <  2:1332347/93‑1 (MQ=255)

                      aaTCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATc                                                                   >  1:1525996/1‑93 (MQ=255)

                       aTCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCa                                                                  <  1:2364395/93‑1 (MQ=255)

                       aTCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCa                                                                  >  1:4258111/1‑93 (MQ=255)

                       aTCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCa                                                                  <  2:3154558/93‑1 (MQ=255)

                        tctcGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCaa                                                                 >  2:2769331/1‑93 (MQ=255)

                         ctcGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCaaa                                                                >  1:52830/1‑93 (MQ=255)

                          tcGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCGGCCAGATTGAAAAAGCCTACCCCAGAGATAaaaa                                                               >  1:2225953/1‑93 (MQ=255)

                           cGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCaaaaa                                                              <  2:3585899/93‑1 (MQ=255)

                           cGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCaaaaa                                                              <  1:1174966/93‑1 (MQ=255)

                            gTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCaaaaaa                                                             >  1:4589607/1‑93 (MQ=255)

                            gTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCaaaaaa                                                             <  1:4615781/93‑1 (MQ=255)

                            gTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCaaaaaa                                                             >  1:4347878/1‑93 (MQ=255)

                               tGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGCTCAAAAAAgcg                                                          >  1:174153/1‑93 (MQ=255)

                               tGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAgcg                                                          >  1:3540998/1‑93 (MQ=255)

                               tGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAgcg                                                          <  2:1182446/93‑1 (MQ=255)

                               tGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAgcg                                                          >  2:308081/1‑93 (MQ=255)

                               tGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAgcg                                                          >  1:2813078/1‑93 (MQ=255)

                                gAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAgcgg                                                         >  2:635922/1‑93 (MQ=255)

                                gAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAgcgg                                                         >  1:2803197/1‑93 (MQ=255)

                                   aGTGTCTGTTTCTATCTTCGGTCGTGCGCCCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCAAAAAAgcggcgt                                                      >  1:4045538/1‑92 (MQ=255)

                                    gtgtCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGAt                                                     <  2:3290958/93‑1 (MQ=255)

                                    gtgtCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGAt                                                     >  2:45520/1‑93 (MQ=255)

                                     tgtCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTACCCCGGCGATCAAAAAAGCGGCGAtt                                                    <  2:2365388/93‑1 (MQ=255)

                                     tgtCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGAtt                                                    >  1:901893/1‑93 (MQ=255)

                                     tgtCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGAtt                                                    <  2:3241682/93‑1 (MQ=255)

                                     tgtCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGAtt                                                    <  2:2730153/93‑1 (MQ=255)

                                     tgtCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGAtt                                                    >  2:4257740/1‑93 (MQ=255)

                                         tGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAAt                                                >  1:4606481/1‑93 (MQ=255)

                                         tGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAAt                                                <  1:1581937/93‑1 (MQ=255)

                                         tGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAAt                                                <  2:596540/93‑1 (MQ=255)

                                         tGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAAt                                                <  2:1406288/93‑1 (MQ=255)

                                           tttCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCg                                              >  2:340787/1‑93 (MQ=255)

                                            ttCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGt                                             <  1:2084599/93‑1 (MQ=255)

                                            ttCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGt                                             <  1:3477151/93‑1 (MQ=255)

                                            ttCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACCTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGt                                             <  1:3763172/93‑1 (MQ=255)

                                             tCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGtt                                            >  2:3463881/1‑93 (MQ=255)

                                             tCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGtt                                            <  1:4633388/93‑1 (MQ=255)

                                             tCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGtt                                            >  1:956664/1‑93 (MQ=255)

                                             tCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGtt                                            >  1:1055061/1‑93 (MQ=255)

                                               tATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGc                                          >  1:587067/1‑93 (MQ=255)

                                                aTCTTCGGTCGTGTGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGca                                         <  2:2068321/93‑1 (MQ=255)

                                                aTCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGca                                         >  2:2634589/1‑93 (MQ=255)

                                                 tCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGcac                                        >  2:1558194/1‑93 (MQ=255)

                                                  cTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGcaca                                       <  2:351884/93‑1 (MQ=255)

                                                   ttCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGcacag                                      >  2:1736888/1‑92 (MQ=255)

                                                   ttCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACaa                                      <  1:987005/93‑1 (MQ=255)

                                                    tCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAg                                     <  2:1492866/93‑1 (MQ=255)

                                                     cGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAgg                                    >  1:2302351/1‑93 (MQ=255)

                                                      ggTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGc                                   <  2:4415866/93‑1 (MQ=255)

                                                       gTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCg                                  >  1:104958/1‑93 (MQ=255)

                                                         cGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTg                                <  1:1618223/93‑1 (MQ=255)

                                                          gTGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTga                               <  1:1780810/93‑1 (MQ=255)

                                                           tGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTgag                              <  2:2963855/93‑1 (MQ=255)

                                                           tGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTgag                              >  2:4153381/1‑93 (MQ=255)

                                                           tGCGACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTgag                              >  1:4021015/1‑93 (MQ=255)

                                                            gCGACCCCGGTAGAGCTGGACTTCAGCCATATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTgaga                             <  2:1009381/93‑1 (MQ=255)

                                                              gACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGAtt                           <  2:2322328/93‑1 (MQ=255)

                                                              gACCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGAtt                           >  2:1057377/1‑93 (MQ=255)

                                                               aCCCCGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTg                          <  1:4212866/93‑1 (MQ=255)

                                                                ccccGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGACTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTgg                         <  1:2142538/93‑1 (MQ=255)

                                                                   cGGTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAAt                      <  1:1972592/93‑1 (MQ=255)

                                                                    ggTAGAGCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATa                     <  2:2501745/93‑1 (MQ=255)

                                                                        gagCTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAAt                 >  2:1585841/1‑93 (MQ=255)

                                                                           cTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTc              <  2:3072076/93‑1 (MQ=255)

                                                                           cTGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTc              >  1:3221138/1‑93 (MQ=255)

                                                                            tGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAAtttag             >  2:3015411/1‑91 (MQ=255)

                                                                            tGGACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTcg             >  2:2122940/1‑93 (MQ=255)

                                                                             ggACTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTcgc            <  1:4431368/93‑1 (MQ=255)

                                                                               aCTTCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCgcgc          <  1:3170661/93‑1 (MQ=255)

                                                                                  tCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCtt       <  2:1207048/93‑1 (MQ=255)

                                                                                  tCAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCtt       >  2:859756/1‑93 (MQ=255)

                                                                                   cAGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCttt      <  2:4578276/93‑1 (MQ=255)

                                                                                    aGCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCCGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCtttt     <  2:4050181/93‑1 (MQ=255)

                                                                                     gCCAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTg    <  2:3025112/93‑1 (MQ=255)

                                                                                      ccAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGt   >  2:3247225/1‑93 (MQ=255)

                                                                                       cAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGtt  >  1:3430599/1‑93 (MQ=255)

                                                                                       cAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGtt  >  1:3430600/1‑93 (MQ=255)

                                                                                       cAGATTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGtt  <  1:1400270/93‑1 (MQ=255)

GTTGAAGAAGTGGATTACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTT  >  NC_000913/4175763‑4175942

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