BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	526,303	A→G	G1273G (GGA→GGG)	rhsD →	RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	526,303	0	A	G	99.4%	517.6 / NA	162	G1273G (GGA→GGG)	rhsD	RhsD
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (83/78); minor base T (0/1); total (83/79)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

CAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑C‑‑GGG‑GGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAA  >  NC_000913/526214‑526392

     tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTtgtg                                                                                    >  2:5813869/1‑93 (MQ=255)

     tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑C‑‑GGG‑GGGGGTTTgtgg                                                                                   >  1:307468/1‑93 (MQ=255)

      ccTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTgtgg                                                                                   >  2:4248101/1‑93 (MQ=255)

       cTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTgtggt                                                                                  <  2:5912177/93‑1 (MQ=255)

        tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTgtggtg                                                                                 >  2:3257219/1‑93 (MQ=255)

        tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTgtggtg                                                                                 >  1:2131176/1‑93 (MQ=255)

        tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTgtgggg                                                                                 >  2:3297235/1‑91 (MQ=255)

         ttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGgtgt                                                                                >  1:1764080/1‑93 (MQ=255)

           aaaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTc                                                                              <  2:1951139/93‑1 (MQ=255)

             aTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCtt                                                                            >  1:2404712/1‑93 (MQ=255)

              tCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGAGGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCttt                                                                           <  1:2544208/93‑1 (MQ=255)

                cACTACAGCAAATTGCCCCTATGGGATTATTGCGGACTGGGGATGATGCCAGATCTGGAGCA‑TCT‑AGG‑‑GGG‑GGGGGTTTGTGGTGTTCTTTCa                                                                         <  2:481367/93‑1 (MQ=255)

                cACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCa                                                                         <  2:1026541/93‑1 (MQ=255)

                cACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCa                                                                         >  2:5031114/1‑93 (MQ=255)

                 aCTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCAc                                                                        >  1:6027434/1‑93 (MQ=255)

                 aCTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTATTTCaa                                                                        >  1:846754/1‑92 (MQ=255)

                 aCTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGGTTGTGGTGTTCTTTCAc                                                                        >  2:1450725/1‑93 (MQ=255)

                  cTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACg                                                                       <  2:5219951/93‑1 (MQ=255)

                  cTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACg                                                                       >  1:3148582/1‑93 (MQ=255)

                   tACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGt                                                                      >  2:2740209/1‑93 (MQ=255)

                      aacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGtata                                                                   <  1:3059372/93‑1 (MQ=255)

                      aacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGtata                                                                   >  2:4761060/1‑93 (MQ=255)

                       acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTataa                                                                  >  1:5340545/1‑93 (MQ=255)

                       acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTataa                                                                  >  1:1173173/1‑93 (MQ=255)

                       acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑C‑‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTataat                                                                 <  2:3438127/93‑1 (MQ=255)

                       acaaATTGACACTATAGGATTATTGCAAGCTTGGGAGGATGCAAGTCCGGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTataa                                                                  <  2:240506/93‑1 (MQ=255)

                        caaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGCGGTGTTCTTTCACGTataat                                                                 <  2:3585219/93‑1 (MQ=255)

                         aaaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTataata                                                                <  1:4624773/93‑1 (MQ=255)

                          aaTTGGCCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCGGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAg                                                               <  2:5483587/93‑1 (MQ=255)

                          aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAg                                                               <  1:3786788/93‑1 (MQ=255)

                          aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGGTTGTGGTGTTCTTTCACGTATAATAg                                                               >  2:4636178/1‑93 (MQ=255)

                           aTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAgg                                                              <  1:1864337/93‑1 (MQ=255)

                            ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGa                                                             <  1:5035699/93‑1 (MQ=255)

                            ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGGGTTCTTTCACGTATAATAGGa                                                             >  1:1299379/1‑93 (MQ=255)

                             tGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGAc                                                            >  1:3890187/1‑93 (MQ=255)

                             tGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTGCTTTCACGTATAATAGGAc                                                            <  2:6144081/93‑1 (MQ=255)

                              gACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATATGAcc                                                           >  2:4349823/1‑93 (MQ=255)

                               aCCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCa                                                          <  1:2712256/93‑1 (MQ=255)

                                gccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCaa                                                         <  1:5911840/92‑1 (MQ=255)

                                 ccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGGTTGCGGCGTTATTTCCCGTATAATAGGACCAAg                                                        >  2:4937241/1‑93 (MQ=255)

                                   tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                                      >  2:3129287/1‑93 (MQ=255)

                                   tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                                      >  1:1917284/1‑93 (MQ=255)

                                   tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                                      >  2:3946269/1‑93 (MQ=255)

                                   tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGCTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                                      >  2:927185/1‑93 (MQ=255)

                                    aTGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGGTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaa                                                     >  2:541435/1‑93 (MQ=255)

                                    aTGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑C‑‑GGG‑GGGGGCTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                                    >  2:1755225/1‑93 (MQ=255)

                                      gggATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAt                                                   <  2:5589483/93‑1 (MQ=255)

                                       ggATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑C‑‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAttt                                                 >  1:2587073/1‑93 (MQ=255)

                                        gATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAttt                                                 >  1:3734143/1‑93 (MQ=255)

                                         attattGCAGACTTGGGATGATCCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTg                                                <  2:3733296/93‑1 (MQ=255)

                                          ttattGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGa                                               <  1:1229032/93‑1 (MQ=255)

                                          ttattGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGa                                               <  2:2282536/93‑1 (MQ=255)

                                           tattGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑C‑‑GGGTGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGAt                                              >  2:3872396/1‑93 (MQ=255)

                                           tatgGCAGTGTTGGGATGATGCAAGAGGTGGAGCA‑TGT‑A‑CGGGGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGa                                               <  1:641304/93‑1 (MQ=255)

                                            attGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATa                                             >  2:1486245/1‑93 (MQ=255)

                                            attGCAGACTTGGGATGATGCCAGATCTGAAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATa                                             <  1:5851170/93‑1 (MQ=255)

                                               gCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTa                                          <  1:909482/93‑1 (MQ=255)

                                               gCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTa                                          <  1:4675748/93‑1 (MQ=255)

                                                gaGACTTGGGATGATGGCAGATCTGGAGCA‑TGTAA‑G‑‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTAc                                         <  1:5264324/92‑1 (MQ=255)

                                                 agtgttGGGATGAGGCCAGATGTGGACCAGTGC‑G‑G‑‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                                        <  1:1169877/89‑1 (MQ=255)

                                                 aGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑C‑‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                                       <  2:50422/93‑1 (MQ=255)

                                                   aCTTGGGATTATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGc                                      <  1:1939643/93‑1 (MQ=255)

                                                   aCTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGc                                      <  1:4235577/93‑1 (MQ=255)

                                                   aCTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑C‑‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCa                                     >  1:5826917/1‑93 (MQ=255)

                                                     ttGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAg                                    >  2:4180404/1‑93 (MQ=255)

                                                      tGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTATGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGa                                   >  1:3018669/1‑93 (MQ=255)

                                                       gggATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAt                                  >  1:5495235/1‑93 (MQ=255)

                                                        ggATGATGCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtg                                 <  2:2517080/93‑1 (MQ=255)

                                                          atgatgGCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgct                               <  1:905862/93‑1 (MQ=255)

                                                             atgCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCg                            <  2:1619359/93‑1 (MQ=255)

                                                             atgCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCg                            >  2:4107656/1‑93 (MQ=255)

                                                              tgCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGt                           >  1:1451079/1‑93 (MQ=255)

                                                               gCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGtt                          >  1:3952365/1‑93 (MQ=255)

                                                               gCCAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGtt                          >  1:4155595/1‑93 (MQ=255)

                                                                ccAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                         >  2:89575/1‑93 (MQ=255)

                                                                ccAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                         <  1:3722142/93‑1 (MQ=255)

                                                                ccAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                         >  2:6247156/1‑93 (MQ=255)

                                                                 cAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                        >  1:5684925/1‑93 (MQ=255)

                                                                 cAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                        <  1:4988427/93‑1 (MQ=255)

                                                                 cAGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                        >  2:2803716/1‑93 (MQ=255)

                                                                  aGATCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGa                       <  2:2860378/93‑1 (MQ=255)

                                                                     tCTGGAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGc                    <  1:4833250/93‑1 (MQ=255)

                                                                        ggAGCA‑TGT‑G‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt                 <  2:126007/93‑1 (MQ=255)

                                                                        ggAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt                 >  2:2047505/1‑93 (MQ=255)

                                                                        ggAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt                 <  2:824479/93‑1 (MQ=255)

                                                                        ggAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt                 <  1:3684297/93‑1 (MQ=255)

                                                                         gAGCA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTg                <  1:996373/93‑1 (MQ=255)

                                                                            cA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaa             >  1:4970688/1‑93 (MQ=255)

                                                                            cA‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaa             >  2:2221483/1‑93 (MQ=255)

                                                                             a‑TGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaag            <  1:4001876/93‑1 (MQ=255)

                                                                               tGT‑A‑CG‑GGG‑GGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga           >  2:2355386/1‑93 (MQ=255)

                                                                                    tccgggg‑gggggTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaaa         <  2:2942308/89‑1 (MQ=255)

                                                                                    tacgggg‑gggggTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaaa         <  2:3812146/89‑1 (MQ=255)

                                                                                    tacaggg‑gggGCTTTGTGGTGTTCTTTCACATATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaaa         <  2:5828047/89‑1 (MQ=255)

                                                                                     tcgggg‑gggggTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAAc        <  2:59447/90‑1 (MQ=255)

                                                                                       gggg‑gggggTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACga      >  2:1114798/2‑92 (MQ=255)

                                                                                         gg‑ggggGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAg    >  2:5891918/1‑93 (MQ=255)

                                                                                          g‑gggggTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGa   >  1:6241374/1‑93 (MQ=255)

                                                                                            ggggGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGaa  <  1:6133707/93‑1 (MQ=255)

                                                                                            ggggGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGaa  >  1:5161862/1‑93 (MQ=255)

CAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGT‑A‑C‑‑GGG‑GGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAA  >  NC_000913/526214‑526392

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