BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	568,823	A→G	G233G (GGA→GGG)	ybcK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	568,823	0	A	G	99.6%	1494.2 / ‑3.8	464	G233G (GGA→GGG)	ybcK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (253/208); total (255/209)
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.87e-01

Only 100 of 476 total aligned reads displayed.

CAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGA  >  NC_000913/568733‑568912

  ggATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCtt                                                                                       >  1:1309761/1‑93 (MQ=255)

   gATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTc                                                                                      >  2:3072923/1‑93 (MQ=255)

    gtGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCt                                                                                     <  2:8264691/92‑1 (MQ=255)

    aTGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCt                                                                                     <  2:7335232/93‑1 (MQ=255)

     tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTg                                                                                    >  1:1230977/1‑93 (MQ=255)

     tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTg                                                                                    >  1:9106252/1‑93 (MQ=255)

     tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTg                                                                                    <  2:7518584/93‑1 (MQ=255)

      ggAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGt                                                                                   >  2:7226878/1‑93 (MQ=255)

        aaaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTaa                                                                                 >  1:9477527/1‑93 (MQ=255)

         aaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAt                                                                                <  2:3995769/93‑1 (MQ=255)

           ggCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTg                                                                              >  2:5782447/1‑93 (MQ=255)

             cgcTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaa                                                                            <  2:1199547/93‑1 (MQ=255)

               cTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaa                                                                          >  1:3853182/1‑93 (MQ=255)

                tCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaaa                                                                         >  1:4063632/1‑93 (MQ=255)

                 cATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGTGCCTTCTGTAATTGAAAAat                                                                        <  1:1715765/93‑1 (MQ=255)

                 cATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAat                                                                        >  1:4724850/1‑93 (MQ=255)

                  aTTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatt                                                                       <  1:7231194/93‑1 (MQ=255)

                  aTTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatt                                                                       >  2:6856399/1‑93 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATATCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      >  2:9017874/1‑93 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGATGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGCGGGCCTTGTGTAATTGAAAAatta                                                                      <  2:1854271/93‑1 (MQ=255)

                    tGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAattat                                                                     <  1:2391205/93‑1 (MQ=255)

                     gAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAattatt                                                                    <  2:9726383/93‑1 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   >  2:548487/1‑93 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAAttatta                                                                   <  2:6675859/93‑1 (MQ=255)

                       aTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAg                                                                  <  1:4970837/93‑1 (MQ=255)

                       aTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAg                                                                  >  2:1100512/1‑93 (MQ=255)

                        tGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGc                                                                 >  2:943390/1‑93 (MQ=255)

                         gCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCg                                                                <  1:6210530/93‑1 (MQ=255)

                         gCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCg                                                                >  1:6522765/1‑93 (MQ=255)

                           aaTAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGaa                                                              >  2:4455859/1‑93 (MQ=255)

                            aTAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAAt                                                             <  1:5136456/93‑1 (MQ=255)

                            aTAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAAt                                                             <  1:8232718/93‑1 (MQ=255)

                             tAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATa                                                            >  2:780699/1‑93 (MQ=255)

                              aGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATaa                                                           <  1:3671585/93‑1 (MQ=255)

                              aGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATaa                                                           <  2:3775744/93‑1 (MQ=255)

                               gCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATaaa                                                          >  2:1593873/1‑93 (MQ=255)

                               gCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATaaa                                                          <  2:9393341/93‑1 (MQ=255)

                                ccAAGTATTTAAATGATCATCCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAg                                                         <  1:6655584/93‑1 (MQ=255)

                                 cAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGc                                                        >  2:23039/1‑93 (MQ=255)

                                  aaGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGct                                                       >  1:7372860/1‑93 (MQ=255)

                                  aaGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGct                                                       >  2:2453586/1‑93 (MQ=255)

                                  aaGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCATCTGTAATTGAAAAATTATTAGCGAATAAAGct                                                       >  1:9676288/1‑93 (MQ=255)

                                   aGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctc                                                      <  2:1492608/93‑1 (MQ=255)

                                    gTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctct                                                     <  1:1588737/93‑1 (MQ=255)

                                    gTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctct                                                     >  1:9302761/1‑93 (MQ=255)

                                    gTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctct                                                     <  2:5695802/93‑1 (MQ=255)

                                     tATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTg                                                    >  2:10021921/1‑93 (MQ=255)

                                      aTTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGa                                                   >  1:7730217/1‑93 (MQ=255)

                                       tttAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGAt                                                  <  1:6730197/93‑1 (MQ=255)

                                        ttAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATa                                                 <  1:5827156/93‑1 (MQ=255)

                                         tAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAg                                                >  1:5166709/1‑93 (MQ=255)

                                         tAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAg                                                >  2:5257598/1‑93 (MQ=255)

                                          aaaTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAgg                                               >  2:6471319/1‑93 (MQ=255)

                                            aTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTa                                             <  1:3304089/93‑1 (MQ=255)

                                            aTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTa                                             <  1:6961281/93‑1 (MQ=255)

                                              gATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAtt                                           >  2:2109260/1‑93 (MQ=255)

                                              gATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAtt                                           <  1:9856720/93‑1 (MQ=255)

                                              gATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAtt                                           <  1:5374897/93‑1 (MQ=255)

                                               aTCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAttt                                          <  2:292855/93‑1 (MQ=255)

                                                tCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTg                                         >  2:5986176/1‑93 (MQ=255)

                                                tCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTg                                         <  2:7753891/93‑1 (MQ=255)

                                                 cATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGc                                        >  1:5719277/1‑93 (MQ=255)

                                                    gCTGTCAAGAGTTTCACAGGAAAAGAATGTACATGAGGGCCTTCTGTAATTGACAAATTATTAGCGAATAAAGCTCTGATAGGAATTTGcatt                                     <  2:8957104/93‑4 (MQ=255)

                                                      tGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTAcc                                   >  1:2011558/1‑93 (MQ=255)

                                                       gTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCt                                  <  2:8583325/93‑1 (MQ=255)

                                                       gTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCt                                  >  2:2984042/1‑93 (MQ=255)

                                                       gTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCt                                  >  1:802260/1‑93 (MQ=255)

                                                        tAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCtt                                 >  1:7013465/1‑93 (MQ=255)

                                                        tAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCtt                                 >  1:1121865/1‑93 (MQ=255)

                                                         aaaGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTc                                >  1:4442025/1‑93 (MQ=255)

                                                               tttCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCg                          <  1:4196257/93‑1 (MQ=255)

                                                                 tctcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTg                        >  2:6564796/1‑93 (MQ=255)

                                                                 tctcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTg                        <  2:2792872/93‑1 (MQ=255)

                                                                  ctcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGc                       >  1:8748768/1‑93 (MQ=255)

                                                                   tcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCa                      <  2:3319228/93‑1 (MQ=255)

                                                                   tcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCa                      >  1:3557028/1‑93 (MQ=255)

                                                                   tcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCa                      >  1:4388545/1‑93 (MQ=255)

                                                                    cAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCaa                     <  1:3965188/93‑1 (MQ=255)

                                                                      tgAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAaga                   <  1:4590223/92‑1 (MQ=255)

                                                                      ggAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAaga                   <  1:2984042/93‑1 (MQ=255)

                                                                      ggAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAaga                   >  1:2740483/1‑93 (MQ=255)

                                                                        aaaaGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAgg                 >  2:4365241/1‑93 (MQ=255)

                                                                         aaagaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAggg                >  1:5897338/1‑93 (MQ=255)

                                                                         aaagaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAggg                >  2:2722874/1‑93 (MQ=255)

                                                                           agaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGaa              >  2:3704258/1‑93 (MQ=255)

                                                                            gaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGaaa             <  2:7003615/93‑1 (MQ=255)

                                                                            gaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGaaa             <  2:3456138/93‑1 (MQ=255)

                                                                             aaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAg            <  1:2179458/93‑1 (MQ=255)

                                                                             aaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAg            <  1:6776729/93‑1 (MQ=255)

                                                                             aaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAg            >  1:4844988/1‑93 (MQ=255)

                                                                              aagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAgg           >  2:4991840/1‑93 (MQ=255)

                                                                                gTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGa         >  2:8851774/1‑93 (MQ=255)

                                                                                gTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGa         <  2:7903550/93‑1 (MQ=255)

                                                                                  gCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATa       >  1:10047508/1‑93 (MQ=255)

                                                                                  gCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATa       <  2:6175444/93‑1 (MQ=255)

                                                                                    aTGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAg     <  2:7645741/93‑1 (MQ=255)

                                                                                    aTGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAg     >  1:4169680/1‑93 (MQ=255)

                                                                                    aTGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAg     <  1:6895719/93‑1 (MQ=255)

                                                                                     tGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGt    <  1:7494018/93‑1 (MQ=255)

                                                                                      ggggGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTg   >  2:7618177/1‑93 (MQ=255)

CAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGA  >  NC_000913/568733‑568912

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