BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	526,302	(G)_6→7	coding (3818/4281 nt)	rhsD →	RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	526,296	1	.	G	95.9%	454.9 / 11.4	147	coding (3812/4281 nt)	rhsD	RhsD
Reads supporting (aligned to +/- strand): ref base . (2/4); new base G (93/47); total (96/51)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 172 total aligned reads displayed.

GAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGA  >  NC_000913/526204‑526383

gAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑c                                                                                             <  1:4614703/93‑1 (MQ=255)

 aaTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑g                                                                                          <  2:4164648/93‑1 (MQ=255)

  aTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg                                                                                         <  2:2850815/93‑1 (MQ=255)

  aTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg                                                                                         <  2:1475124/93‑1 (MQ=255)

   tttATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg‑g                                                                                       >  1:4227894/1‑93 (MQ=255)

    ttATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg‑gg                                                                                      <  2:3260097/93‑1 (MQ=255)

     tataTCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg‑ggg                                                                                     <  2:4058239/93‑1 (MQ=255)

      atatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg‑gggg                                                                                    >  2:3601427/1‑93 (MQ=255)

       tatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg‑ggggg                                                                                   >  1:3626973/1‑92 (MQ=255)

       tatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg‑ggggg                                                                                   >  2:1738463/1‑92 (MQ=255)

        atCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg‑ggggga                                                                                  >  1:3748265/1‑91 (MQ=255)

        atCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGAC‑C‑‑gg‑gggggg                                                                                  >  1:3149092/1‑91 (MQ=255)

         tCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑gg‑gggggag                                                                                 >  2:4575746/1‑90 (MQ=255)

            gTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGttt                                                                              >  1:160881/1‑93 (MQ=255)

              aTCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTtgt                                                                            >  2:2873512/1‑93 (MQ=255)

              aTCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTtgt                                                                            >  2:541360/1‑93 (MQ=255)

              aTCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTtgt                                                                            >  2:4288189/1‑93 (MQ=255)

               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑GGTTtgtg                                                                           >  2:2426259/1‑93 (MQ=255)

               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTtgtg                                                                           >  1:4544080/1‑93 (MQ=255)

                ccTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTgtgg                                                                          >  2:2174429/1‑93 (MQ=255)

                 cTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTgtggt                                                                         >  1:2202791/1‑93 (MQ=255)

                  tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTgtggtg                                                                        >  2:2053431/1‑93 (MQ=255)

                  tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTgtggtg                                                                        >  1:3410967/1‑93 (MQ=255)

                   ttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGgtgt                                                                       <  2:3043600/93‑1 (MQ=255)

                    tAAATCCACTGCAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGtt                                                                      <  1:910072/93‑1 (MQ=255)

                    tAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGtt                                                                      >  2:751558/1‑93 (MQ=255)

                    tAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGtt                                                                      >  1:4502112/1‑93 (MQ=255)

                      aaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTATTCt                                                                    <  2:4094335/93‑1 (MQ=255)

                        tCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCttt                                                                  >  2:19629/1‑93 (MQ=255)

                          cACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCa                                                                >  1:1245044/1‑93 (MQ=255)

                           aCTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCAc                                                               >  2:2227593/1‑93 (MQ=255)

                            cTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACg                                                              >  2:3574759/1‑93 (MQ=255)

                            cTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACg                                                              <  2:1034913/93‑1 (MQ=255)

                             tACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGt                                                             >  1:2362563/1‑93 (MQ=255)

                               caacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGtat                                                           >  1:4377734/1‑93 (MQ=255)

                               caacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGtata                                                          >  1:4018072/1‑93 (MQ=255)

                                aacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGtata                                                          >  2:3171050/1‑93 (MQ=255)

                                aacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGtata                                                          <  1:2879051/93‑1 (MQ=255)

                                 acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTataa                                                         >  2:2980401/1‑93 (MQ=255)

                                  caaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTataat                                                        <  2:2542801/93‑1 (MQ=255)

                                  caaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTataat                                                        <  1:4456047/93‑1 (MQ=255)

                                   aaaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTataata                                                       <  2:3145356/93‑1 (MQ=255)

                                   aaaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTataata                                                       <  1:2152703/93‑1 (MQ=255)

                                   aaaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAg                                                      <  1:4030831/93‑1 (MQ=255)

                                      ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAggc                                                    >  2:4057274/1‑92 (MQ=255)

                                      ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGa                                                    <  1:3990859/93‑1 (MQ=255)

                                      ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGa                                                    >  1:1827583/1‑93 (MQ=255)

                                          cccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCag                                                >  1:2871794/1‑92 (MQ=255)

                                          cccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCaa                                                >  2:1419215/1‑93 (MQ=255)

                                           ccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAg                                               <  2:779120/93‑1 (MQ=255)

                                           ccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAg                                               >  2:1034976/1‑93 (MQ=255)

                                            cTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAga                                              >  1:4330682/1‑92 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑GGGTTGGGGTGTTCTTTCACGTATAATAGGACCAAgca                                             >  1:3356307/1‑91 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                             >  2:3566814/1‑93 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                             >  1:2094642/1‑93 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CA‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                             >  2:4535903/1‑93 (MQ=255)

                                               tGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                           >  1:2402198/1‑93 (MQ=255)

                                               tGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                           >  1:1294446/1‑93 (MQ=255)

                                                 ggATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAtt                                         >  2:2204924/1‑93 (MQ=255)

                                                 ggATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAtt                                         <  1:1098332/93‑1 (MQ=255)

                                                  gATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAttt                                        >  2:369196/1‑93 (MQ=255)

                                                   attattGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTg                                       >  2:853369/1‑93 (MQ=255)

                                                      attGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATa                                    >  2:2818008/1‑93 (MQ=255)

                                                       ttGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                   >  1:2162616/1‑93 (MQ=255)

                                                        tGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGt                                  <  2:1809434/93‑1 (MQ=255)

                                                           aGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                               <  2:2774554/93‑1 (MQ=255)

                                                           aGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                               <  1:3522407/93‑1 (MQ=255)

                                                            gACTTGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                              >  1:1835508/1‑93 (MQ=255)

                                                              cTTGGGATGCTGGCAGATCTGGAGCATGTA‑C‑‑GG‑GGGGCAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCa                            <  1:1390991/93‑1 (MQ=255)

                                                                tGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGa                          <  1:708252/93‑1 (MQ=255)

                                                                tGGGATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGa                          <  1:3525712/93‑1 (MQ=255)

                                                                 gggATGATGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAt                         >  1:1073704/1‑93 (MQ=255)

                                                                   gatgatGCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                       >  2:1070316/1‑93 (MQ=255)

                                                                   gatgatGCCAGATCTGGAGCATATA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                       <  1:1530836/93‑1 (MQ=255)

                                                                     tgatgCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctg                     >  1:1326033/1‑93 (MQ=255)

                                                                     tgatgCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGGGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctg                     >  1:242856/1‑93 (MQ=255)

                                                                      gatgCCAGATCTGGAGCATGTA‑CGGGG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctg                     <  2:733083/93‑1 (MQ=255)

                                                                      gatgCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctgc                    <  2:2201159/93‑1 (MQ=255)

                                                                       atgCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCg                   <  2:3840951/93‑1 (MQ=255)

                                                                        tgCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGt                  >  2:1629217/1‑93 (MQ=255)

                                                                        tgCCAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGt                  >  1:3767760/1‑93 (MQ=255)

                                                                         gCCAGATCTGGAGCATGTA‑C‑‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                <  2:1594852/93‑1 (MQ=255)

                                                                          gcAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                <  1:2405052/92‑1 (MQ=255)

                                                                          ccAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGGGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                >  2:2333632/1‑93 (MQ=255)

                                                                           cAGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg               >  2:4687813/1‑93 (MQ=255)

                                                                            cgATCAGGAGCATGTA‑C‑‑GGAGGGG‑AGATTGTGGAGTTCAATCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAAGCTGAGTTTGa              <  2:4219710/92‑1 (MQ=255)

                                                                            aGATCTGGAGCCTGGACC‑‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGa              <  1:879856/93‑1 (MQ=255)

                                                                            aGATCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGa              >  2:1306335/1‑93 (MQ=255)

                                                                              aTCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATg            >  1:2253251/1‑93 (MQ=255)

                                                                              aTCTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATg            >  1:3871422/1‑93 (MQ=255)

                                                                                cTGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCt          <  2:3676783/93‑1 (MQ=255)

                                                                                 tGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCtt         >  1:2233632/1‑93 (MQ=255)

                                                                                 tGGAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTCAGATGCtt         >  1:4661492/1‑93 (MQ=255)

                                                                                  ggAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt        >  1:1799920/1‑93 (MQ=255)

                                                                                  ggAGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt        <  1:3593929/93‑1 (MQ=255)

                                                                                    aGCATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTga      >  1:1721332/1‑93 (MQ=255)

                                                                                      cATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaa    <  2:4515819/93‑1 (MQ=255)

                                                                                      cATGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaa    <  1:718637/93‑1 (MQ=255)

                                                                                       aTGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaag   <  1:2226389/93‑1 (MQ=255)

                                                                                        tGTA‑CG‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga  >  2:3401694/1‑93 (MQ=255)

GAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑‑GG‑GGGG‑AGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGA  >  NC_000913/526204‑526383

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