BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,266,308	T→C	S85P (TCT→CCT)	yhaC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,266,308	0	T	C	99.4%	514.7 / ‑4.9	166	S85P (TCT→CCT)	yhaC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (96/69); total (97/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 183 total aligned reads displayed.

ACGTTACAACCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAA  >  NC_000913/3266218‑3266397

 cGTTACAACCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATgtctcgcctc                                                                                        <  2:915060/93‑1 (MQ=255)

  gTTACAACCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCctct                                                                                       >  1:1422229/1‑93 (MQ=255)

  gTTACAACCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCctct                                                                                       <  2:1182213/93‑1 (MQ=255)

     aCAACCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCg                                                                                    >  1:2837805/1‑93 (MQ=255)

      cAACCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGt                                                                                   >  1:1281490/1‑93 (MQ=255)

      cAACCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGt                                                                                   <  2:587258/93‑1 (MQ=255)

      cAACCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGt                                                                                   <  1:960626/93‑1 (MQ=255)

       aaCCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTc                                                                                  <  1:3463209/93‑1 (MQ=255)

       aaCCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTc                                                                                  >  1:1899912/1‑93 (MQ=255)

        aCCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCt                                                                                 >  1:3347116/1‑93 (MQ=255)

            ggtggAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGAt                                                                             >  2:3494565/1‑93 (MQ=255)

             gtggAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGAtt                                                                            <  2:3208060/93‑1 (MQ=255)

                  aaaTGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTa                                                                       <  1:2362940/93‑1 (MQ=255)

                   aaTGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAg                                                                      >  1:314949/1‑93 (MQ=255)

                   aaTGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAg                                                                      <  1:514953/93‑1 (MQ=255)

                   aaTGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAg                                                                      >  2:3062537/1‑93 (MQ=255)

                    aTGCATTATCACTGGTGATTTTCAAAAAGAACTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAgg                                                                     >  2:2877133/1‑93 (MQ=255)

                      gCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGtt                                                                   <  2:960485/93‑1 (MQ=255)

                         ttATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTtata                                                                >  2:1935581/1‑93 (MQ=255)

                         ttATCACTGGGGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGGTTTTTTAGGTtata                                                                >  2:3834083/1‑93 (MQ=255)

                          tATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATaa                                                               <  1:2707100/93‑1 (MQ=255)

                          tATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATaa                                                               <  1:2856131/93‑1 (MQ=255)

                          tATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATaa                                                               <  2:429079/93‑1 (MQ=255)

                            tCACTGGTGATTTTCAAAAAGAATTAGCGGACTTACAAAATGCAATTCTTTATAATGTCTCGCCTCTTTGTCTTGATTTTTTAGTCTATAACg                                                             >  2:2990956/1‑93 (MQ=255)

                             cACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGc                                                            <  2:2241600/93‑1 (MQ=255)

                             cACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGc                                                            <  2:4313158/93‑1 (MQ=255)

                               cTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGccc                                                          <  2:4341697/93‑1 (MQ=255)

                                tGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGTGTATAACGCCCa                                                         >  1:4118036/1‑93 (MQ=255)

                                 ggTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCaa                                                        >  1:4489566/1‑93 (MQ=255)

                                   tgtttttCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAat                                                      >  2:653124/4‑93 (MQ=255)

                                    gATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAatt                                                     >  2:55790/1‑93 (MQ=255)

                                     aTTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAatta                                                    >  1:2008785/1‑93 (MQ=255)

                                      ttttCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAattat                                                   >  1:873123/1‑93 (MQ=255)

                                      ttttCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAattat                                                   >  1:715209/1‑93 (MQ=255)

                                       tttCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGCCTGGATTTTTTAGGTTATAACGCCCAAattatt                                                  >  2:327051/1‑93 (MQ=255)

                                          cAAAAAGAATTAGTGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCaa                                               <  1:3829123/93‑1 (MQ=255)

                                          cAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCaa                                               <  1:4208479/93‑1 (MQ=255)

                                          cAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCaa                                               <  1:4198493/93‑1 (MQ=255)

                                           aaaaaGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAAc                                              >  2:1424028/1‑93 (MQ=255)

                                            aaaaGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACg                                             >  1:553029/1‑93 (MQ=255)

                                             aaaGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGa                                            >  1:4460479/1‑93 (MQ=255)

                                             aaaGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGa                                            >  2:2696496/1‑93 (MQ=255)

                                              aaGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGAt                                           >  1:2919790/1‑93 (MQ=255)

                                               aGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATc                                          >  1:2609901/1‑93 (MQ=255)

                                               aGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATc                                          <  1:679583/93‑1 (MQ=255)

                                                gAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCg                                         <  2:1242592/93‑1 (MQ=255)

                                                gAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCg                                         <  1:2424132/93‑1 (MQ=255)

                                                gAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCg                                         >  1:304155/1‑93 (MQ=255)

                                                 aaTTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCgg                                        >  2:1151247/1‑93 (MQ=255)

                                                 aaTTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCgg                                        >  1:1064949/1‑93 (MQ=255)

                                                  aTTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGa                                       >  2:397150/1‑93 (MQ=255)

                                                  aTTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGa                                       >  2:780816/1‑93 (MQ=255)

                                                   ttAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGAc                                      <  2:186703/93‑1 (MQ=255)

                                                    tAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAGATTATTCAACGATCGGACa                                     <  2:4576873/93‑1 (MQ=255)

                                                    tAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACa                                     >  1:3710455/1‑93 (MQ=255)

                                                       cGTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATa                                  <  1:1564437/93‑1 (MQ=255)

                                                        gTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATat                                 >  2:3202985/1‑92 (MQ=255)

                                                        gTACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATAc                                 >  1:3253094/1‑93 (MQ=255)

                                                         tACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACt                                <  1:2934903/93‑1 (MQ=255)

                                                         tACTTACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACt                                >  2:1541625/1‑93 (MQ=255)

                                                            ttACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTtgt                             <  2:134046/93‑1 (MQ=255)

                                                            ttACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTtgt                             >  1:1191245/1‑93 (MQ=255)

                                                            ttACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTtgt                             >  1:1770534/1‑93 (MQ=255)

                                                            ttACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTtgt                             <  2:4101149/93‑1 (MQ=255)

                                                            ttACAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTtgt                             <  2:22038/93‑1 (MQ=255)

                                                              aCAAAATGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGa                           >  2:1462214/1‑93 (MQ=255)

                                                                aaaaTGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAAc                         >  1:1584654/1‑93 (MQ=255)

                                                                aaaaTGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAAc                         <  2:3111758/93‑1 (MQ=255)

                                                                aaaaTGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAAc                         >  1:3322760/1‑93 (MQ=255)

                                                                aaaaTGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAAc                         >  2:1118362/1‑93 (MQ=255)

                                                                aaaaTGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGAGAAc                         >  2:4047002/1‑93 (MQ=255)

                                                                   aTGTAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTa                      >  1:2058575/1‑93 (MQ=255)

                                                                      tAATTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCa                   <  1:1362371/93‑1 (MQ=255)

                                                                       aaTTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCAt                  >  1:4032819/1‑93 (MQ=255)

                                                                       aaTTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCAt                  >  1:1434857/1‑93 (MQ=255)

                                                                        aTTCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCAtt                 >  1:252733/1‑93 (MQ=255)

                                                                         ttCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTa                >  2:2305416/1‑93 (MQ=255)

                                                                          tCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTaa               <  1:478201/93‑1 (MQ=255)

                                                                          tCTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTaa               <  1:721424/93‑1 (MQ=255)

                                                                           cTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAAt              >  1:1645492/1‑93 (MQ=255)

                                                                           cTTTATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAAt              >  1:3405965/1‑93 (MQ=255)

                                                                            tttATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATg             >  1:4171169/1‑93 (MQ=255)

                                                                            tttATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATg             >  2:2895095/1‑93 (MQ=255)

                                                                             ttATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGa            <  1:581091/93‑1 (MQ=255)

                                                                             ttATAATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGa            >  1:3608826/1‑93 (MQ=255)

                                                                              tataATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGaa           >  1:3554681/1‑93 (MQ=255)

                                                                                taATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAAcc         <  2:3971192/93‑1 (MQ=255)

                                                                                taATGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAAcc         >  1:2172674/1‑93 (MQ=255)

                                                                                 aaTGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCg        >  2:3162089/1‑93 (MQ=255)

                                                                                  aTGTCTCGCCTCTTCGTCTGAATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGt       <  1:3131995/93‑1 (MQ=255)

                                                                                   tGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGtt      <  2:2406115/93‑1 (MQ=255)

                                                                                   tGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGtt      >  2:498253/1‑93 (MQ=255)

                                                                                   tGTCTCGCCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGtt      >  1:2531796/1‑93 (MQ=255)

                                                                                    gtctcgcctcTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTa     <  2:2662225/93‑1 (MQ=255)

                                                                                     tctcgcctcTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTaa    <  2:4268206/93‑1 (MQ=255)

                                                                                     tctcgcctcTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTaa    >  2:3451573/1‑93 (MQ=255)

                                                                                     tctcgcctcTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTaa    <  1:1353421/93‑1 (MQ=255)

                                                                                      ctcgcctcTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTaaa   >  2:4671153/1‑93 (MQ=255)

                                                                                      ctcgcctcTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTaaa   <  1:3699306/93‑1 (MQ=255)

                                                                                       tcgcctcTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTaaaa  <  2:4059206/93‑1 (MQ=255)

ACGTTACAACCTGGTGGAAAATGCATTATCACTGGTGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAA  >  NC_000913/3266218‑3266397

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