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%	853.9 / ‑3.9	260	coding (699/1527 nt)	ybcK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (126/133); minor base . (1/0); total (127/133)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 266 total aligned reads displayed.

CAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGA  >  NC_000913/568733‑568912

 aGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCt                                                                                        <  2:2976630/93‑1 (MQ=255)

 aGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCt                                                                                        >  2:5224738/1‑93 (MQ=255)

 aGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCt                                                                                        <  1:2840107/93‑1 (MQ=255)

 aGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCt                                                                                        >  2:1674348/1‑93 (MQ=255)

 aGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCt                                                                                        <  2:1903857/93‑1 (MQ=255)

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

   gATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTc                                                                                      <  1:5027975/93‑1 (MQ=255)

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

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

    aTGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCt                                                                                     >  1:3156024/1‑93 (MQ=255)

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

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

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

         aaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAt                                                                                >  1:4576117/1‑93 (MQ=255)

          aGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAtt                                                                               >  1:4560522/1‑93 (MQ=255)

          aGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAtt                                                                               >  1:2586732/1‑93 (MQ=255)

          aGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAtt                                                                               <  2:526583/93‑1 (MQ=255)

              gcTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaa                                                                           <  1:1304418/93‑1 (MQ=255)

               cTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaa                                                                          <  1:1716997/93‑1 (MQ=255)

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

               cTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaa                                                                          <  2:1552795/93‑1 (MQ=255)

                 cATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAat                                                                        <  2:218237/93‑1 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      <  2:1130753/93‑1 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      <  2:1341668/93‑1 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      <  1:6545248/93‑1 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      <  2:4509042/93‑1 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      >  1:2219791/1‑93 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      <  1:1974/93‑1 (MQ=255)

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

                    tGAATGCAATAGCCAAGTATTCAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAattat                                                                     >  2:2573921/1‑93 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   <  2:5660578/93‑1 (MQ=255)

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

                        tGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGc                                                                 <  1:6306393/93‑1 (MQ=255)

                          cAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGa                                                               <  2:3840416/93‑1 (MQ=255)

                          cAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGa                                                               <  1:650830/93‑1 (MQ=255)

                          cAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGa                                                               >  1:4038127/1‑93 (MQ=255)

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

                            aTAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAAt                                                             >  2:4773056/1‑93 (MQ=255)

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

                                 cAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGc                                                        <  1:2456558/93‑1 (MQ=255)

                                   aGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctc                                                      >  1:1658915/1‑93 (MQ=255)

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

                                   aGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctc                                                      >  2:3933990/1‑93 (MQ=255)

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

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

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

                                     tATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTg                                                    <  1:5624457/93‑1 (MQ=255)

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

                                      aTTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGa                                                   <  2:4662958/93‑1 (MQ=255)

                                       tttAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGAt                                                  >  2:5017957/1‑93 (MQ=255)

                                       tttAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGAt                                                  >  2:6398995/1‑93 (MQ=255)

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

                                          aaaTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTG‑AATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGt                                              <  2:1502858/93‑1 (MQ=255)

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

                                             tGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAt                                            >  1:3378522/1‑93 (MQ=255)

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

                                               aTCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAttt                                          >  1:3924371/1‑93 (MQ=255)

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

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

                                                 cATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGc                                        <  2:6188704/93‑1 (MQ=255)

                                                   tGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGt                                      <  2:4290010/93‑1 (MQ=255)

                                                    gCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTa                                     >  2:6559692/1‑93 (MQ=255)

                                                    gCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTa                                     <  1:6232675/93‑1 (MQ=255)

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

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

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

                                                      tGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTAcc                                   <  1:4939774/93‑1 (MQ=255)

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

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

                                                        tAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCtt                                 <  1:1770419/93‑1 (MQ=255)

                                                         aaaGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTc                                <  1:4322050/93‑1 (MQ=255)

                                                          aaGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCa                               >  2:6136718/1‑93 (MQ=255)

                                                          aaGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCa                               >  2:5758121/1‑93 (MQ=255)

                                                           aGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTTCCTTCat                              <  1:1011812/93‑1 (MQ=255)

                                                           aGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCat                              >  2:4873891/1‑93 (MQ=255)

                                                           aGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCat                              >  2:144566/1‑93 (MQ=255)

                                                            gAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCata                             <  2:4141325/93‑1 (MQ=255)

                                                            gAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCata                             <  2:224890/93‑1 (MQ=255)

                                                              aTTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATc                           <  2:5436909/93‑1 (MQ=255)

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

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

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

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

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

                                                                       gAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGATTAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAagag                  <  1:592689/93‑1 (MQ=255)

                                                                       gAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAagag                  >  1:1079233/1‑93 (MQ=255)

                                                                        aaaaGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAgg                 <  2:548048/93‑1 (MQ=255)

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

                                                                        aaaaGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAgg                 <  1:6695014/93‑1 (MQ=255)

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

                                                                          aagaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGa               <  1:636651/93‑1 (MQ=255)

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

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

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

                                                                               agTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAggg          <  1:5672934/93‑1 (MQ=255)

                                                                                gTGCATGGTGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGa         <  2:2846821/93‑1 (MQ=255)

                                                                                 ggCATGGGGGCCTTCTGTAATTTAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGAt        <  1:867836/92‑1 (MQ=255)

                                                                                   cATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATaa      <  1:5331783/93‑1 (MQ=255)

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

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

CAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGA  >  NC_000913/568733‑568912

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