BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,613,561	C→T	R27R (CGC→CGT)	yjjV →	putative DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,613,561	0	C	T	99.4%	567.5 / ‑4.8	166	R27R (CGC→CGT)	yjjV	putative DNase
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (87/77); total (87/79)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 171 total aligned reads displayed.

AACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTG  >  NC_000913/4613464‑4613657

aaCAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcgtgc                                                                                                <  2:3751062/100‑1 (MQ=255)

 aCAGTGAGGTTCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcgtgcg                                                                                               <  2:1681922/100‑1 (MQ=255)

  cAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCgg                                                                                              <  2:1147227/100‑1 (MQ=255)

  cAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCgg                                                                                              >  2:2123277/1‑100 (MQ=255)

   aaTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGc                                                                                             <  2:1043154/100‑1 (MQ=255)

   aaTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGc                                                                                             <  1:1202701/100‑1 (MQ=255)

    aTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGca                                                                                            <  1:875081/100‑1 (MQ=255)

     tGAGGTTCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGcac                                                                                           <  1:1525384/100‑1 (MQ=255)

     tGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGcac                                                                                           >  1:651642/1‑100 (MQ=255)

         gATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGc                                                                                       >  1:287467/1‑100 (MQ=255)

          aTCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCg                                                                                      >  2:1837621/1‑100 (MQ=255)

            cTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCggg                                                                                    <  1:1715859/100‑1 (MQ=255)

              ggCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCg                                                                                  >  1:387132/1‑100 (MQ=255)

                cTTGATTTACCGTTTTTTCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTa                                                                                >  2:3869555/1‑100 (MQ=255)

                  tGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAgg                                                                              <  2:1092813/100‑1 (MQ=255)

                   gATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGc                                                                             <  1:3820156/100‑1 (MQ=255)

                    aTTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCa                                                                            >  2:2323681/1‑100 (MQ=255)

                      ttGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAg                                                                          >  2:3431554/1‑100 (MQ=255)

                       tGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGa                                                                         >  1:1189424/1‑100 (MQ=255)

                        gCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGAt                                                                        >  2:2961995/1‑100 (MQ=255)

                        gCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGAt                                                                        <  2:548783/100‑1 (MQ=255)

                         ccGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATc                                                                       <  1:2902514/100‑1 (MQ=255)

                         ccGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATc                                                                       <  1:2390346/100‑1 (MQ=255)

                          cGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCa                                                                      >  2:2541952/1‑100 (MQ=255)

                           gTTTTCTCGACCCCCACTGCCATTTTGATTTACCGACGTTTAGTGGCGATGAAGCGGCCAGCCTGCAACATGCGGCACAAGCGGGCGCAGGCAAGATCAt                                                                     >  2:383024/1‑100 (MQ=255)

                            ttttATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCAtt                                                                    <  1:1013829/100‑1 (MQ=255)

                             tttATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTg                                                                   <  2:1693926/100‑1 (MQ=255)

                             tttATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTg                                                                   <  2:2241802/100‑1 (MQ=255)

                             tttATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTg                                                                   <  1:1188511/100‑1 (MQ=255)

                             tttATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTg                                                                   <  2:450187/100‑1 (MQ=255)

                             tttATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTg                                                                   <  1:1348307/100‑1 (MQ=255)

                             tttATCGACACCCACAGCAATTTTGATTTCCCGCCGTGTAGTGTCGATGAAGGGGCTAGCCTGGAACGTGCGGCACAAGCGGGCGTAGGCAAGATCAGTg                                                                   <  2:898313/100‑1 (MQ=255)

                                 tCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTcc                                                               >  2:1985179/1‑100 (MQ=255)

                                  cGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATCGTTCCg                                                              <  1:2027742/100‑1 (MQ=255)

                                   gACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCgg                                                             >  2:3658786/1‑100 (MQ=255)

                                   gACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCgg                                                             >  1:2337511/1‑100 (MQ=255)

                                    acacCCGCTGTCATTTTCCTAGCCCACTGTGTAGTGGCGATGAAGATGCCGGCTTGGAACGTGCCGCATGAGCGGGCGTAGGCAAGATCATTGTTCCGGc                                                            <  2:1810095/100‑1 (MQ=255)

                                    acacCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGc                                                            <  2:1446876/100‑1 (MQ=255)

                                    acacCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGc                                                            <  1:3533907/100‑1 (MQ=255)

                                     ctcccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCa                                                           <  1:2647472/98‑1 (MQ=255)

                                     cacCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCa                                                           >  2:2814091/1‑100 (MQ=255)

                                       cccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCaag                                                         >  2:2607110/1‑99 (MQ=255)

                                       cccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAc                                                         >  2:1295029/1‑100 (MQ=255)

                                       cccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAc                                                         >  2:3119489/1‑100 (MQ=255)

                                       cccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAc                                                         >  2:2974852/1‑100 (MQ=255)

                                       cccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAc                                                         <  1:2806379/100‑1 (MQ=255)

                                        ccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACa                                                        >  1:88635/1‑100 (MQ=255)

                                        ccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACa                                                        >  1:1443358/1‑100 (MQ=255)

                                        ccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACa                                                        <  2:3096908/100‑1 (MQ=255)

                                           cTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACagag                                                     >  2:2537005/1‑100 (MQ=255)

                                            tGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAgg                                                    <  1:1308946/100‑1 (MQ=255)

                                            tGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAgg                                                    <  1:3228814/100‑1 (MQ=255)

                                             gCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGc                                                   >  2:139626/1‑100 (MQ=255)

                                              ccATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCg                                                  <  1:1106394/100‑1 (MQ=255)

                                              ccATTTGGATTTCCCGCCGTGTAGTGGCGATGATGTGGCCAGCCTGCTACGTGCGGCACCCGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCg                                                  <  2:75675/100‑1 (MQ=255)

                                                aTTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGa                                                >  1:2224172/1‑100 (MQ=255)

                                                     gATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAtt                                           <  1:2866071/100‑1 (MQ=255)

                                                     gATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAtt                                           <  1:3570795/100‑1 (MQ=255)

                                                      aTTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAttt                                          <  1:506189/100‑1 (MQ=255)

                                                      aTTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAttt                                          >  1:2228163/1‑100 (MQ=255)

                                                      aTTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAttt                                          <  2:3520387/100‑1 (MQ=255)

                                                      aTTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAttt                                          >  1:370693/1‑100 (MQ=255)

                                                       tttCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAtttt                                         >  1:300467/1‑100 (MQ=255)

                                                       tttCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAtttt                                         <  1:1674044/100‑1 (MQ=255)

                                                        ttCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCCACAGAGGCGGAAAAttttc                                        >  1:3388831/1‑99 (MQ=255)

                                                        ttCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTg                                        <  1:2586607/100‑1 (MQ=255)

                                                        ttCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTg                                        >  1:4006306/1‑100 (MQ=255)

                                                        ttCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTg                                        >  1:555435/1‑100 (MQ=255)

                                                            cgccgTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCg                                    >  1:752976/1‑100 (MQ=255)

                                                            cgccgTTGAGTGGAGATGAAGAGGCCAGCATGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCg                                    >  2:3300546/1‑100 (MQ=255)

                                                              ccgTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCgtg                                  <  2:2392884/100‑1 (MQ=255)

                                                              ccgTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCgtg                                  <  1:1573007/100‑1 (MQ=255)

                                                                gTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCgtgtg                                <  1:230168/100‑1 (MQ=255)

                                                                  ttAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGtt                              <  1:3117369/100‑1 (MQ=255)

                                                                     gTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGc                           <  2:3972502/100‑1 (MQ=255)

                                                                     gTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGc                           >  2:2466845/1‑100 (MQ=255)

                                                                       ggCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCAt                         >  1:321197/1‑100 (MQ=255)

                                                                       ggCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCAt                         >  1:3182338/1‑100 (MQ=255)

                                                                       ggCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCAt                         >  1:3057640/1‑100 (MQ=255)

                                                                       ggCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCAt                         <  1:2943049/100‑1 (MQ=255)

                                                                       ggCGATGAAGAGGCCAGCCGGCAACGTGCGGCACAACCGGACGTAGGCAAGATCAGTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCAt                         >  2:2835796/1‑100 (MQ=255)

                                                                        gctaTGAAGAGGCCAGCCTGCAACGTGCGGCCCAAGCGGGCGGCGGGAAGATCGTTTTTAGGGCAACATAGGCGGAAACTATTGGCCGTGTGTTGGCGtt                        >  2:3783596/4‑100 (MQ=255)

                                                                          gATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAg                      >  1:130081/1‑100 (MQ=255)

                                                                          gATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGAGAGTAGGCAAGATCATTGTTCCGGCAACCGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAg                      >  2:1992295/1‑100 (MQ=255)

                                                                            tGAGGAGGCGGGCCTGCAAGGGGGGGCACAAGCGGGGCGGGGAAAGATCATGGTTCCGGAAACAGAGGCGGAAAATTTTGCCCGTGTGTGGGCATTAGCg                    <  2:3634833/100‑1 (MQ=255)

                                                                            tGAAGAGGCCGGCCTGCAACGTGCGTCACAGGCGGGGCTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCg                    <  2:3634689/100‑1 (MQ=255)

                                                                            tGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCg                    >  1:4062551/1‑100 (MQ=255)

                                                                               agagGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGaa                 >  2:3755607/1‑100 (MQ=255)

                                                                                  ggCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAAt              <  2:1223501/100‑1 (MQ=255)

                                                                                  ggCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCCGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAAt              <  2:2572365/100‑1 (MQ=255)

                                                                                     cAGCCTGCAACGTGCGGCACAAGCGGGTGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTAt           >  1:417027/1‑100 (MQ=255)

                                                                                     cAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTAt           >  1:3775289/1‑100 (MQ=255)

                                                                                      aGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATc          >  2:560316/1‑100 (MQ=255)

                                                                                      aGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATc          <  1:1902926/100‑1 (MQ=255)

                                                                                       gCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCa         <  2:2860452/100‑1 (MQ=255)

                                                                                        ccTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCaa        <  2:1004334/100‑1 (MQ=255)

                                                                                         cTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAAc       <  1:1139053/100‑1 (MQ=255)

                                                                                           gCAACGTGCGGCACAAGCGGGCGTAGGCAAGAACATTGATCCGGCAGCAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATGAGCGGAAAATTATCAACCg     >  1:3620221/1‑100 (MQ=255)

                                                                                             aaCGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCt   >  1:2768621/1‑100 (MQ=255)

                                                                                              aCGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTg  >  1:1744218/1‑100 (MQ=255)

AACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTG  >  NC_000913/4613464‑4613657

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